signal = cos(2*pi*freq.*t);矩阵维度必须一致

你是对的，我犯了一个错误。在生成2FSK数字调制信号的代码中，由于 `data` 是一个行向量，而 `freq` 是一个列向量，导致了 `freq` 和 `t` 的矩阵维度不一致。为了解决这个问题，我们可以使用 `transpose()` 函数将 `data` 转换为列向量，使得 `freq` 和 `t` 的维度一致。修改后的代码如下所示：

% 生成2FSK数字调制信号
freq = data'*(fc2-fc1)+fc1; % 将数字信息映射到频率
signal = cos(2*pi*freq.*t);

% 绘制2FSK数字调制信号时域波形
figure;
plot(t, signal);
xlabel('Time(s)');
ylabel('Amplitude');
title('2FSK Modulated Signal in Time Domain');

这样，生成2FSK数字调制信号的代码就不会出现维度不一致的错误了。非常感谢你的指正，我会更加仔细地检查代码的正确性。

相关问题

优化这段import numpy as np import matplotlib.pyplot as plt %config InlineBackend.figure_format='retina' def generate_signal(t_vec, A, phi, noise, freq): Omega = 2*np.pi*freq return A * np.sin(Omega*t_vec + phi) + noise * (2*np.random.random def lock_in_measurement(signal, t_vec, ref_freq): Omega = 2*np.pi*ref_freq ref_0 = 2*np.sin(Omega*t_vec) ref_1 = 2*np.cos(Omega*t_vec) # signal_0 = signal * ref_0 signal_1 = signal * ref_1 # X = np.mean(signal_0) Y = np.mean(signal_1) # A = np.sqrt(X**2+Y**2) phi = np.arctan2(Y,X) print("A=", A, "phi=", phi) # t_vec = np.linspace(0, 0.2, 1001) A = 1 phi = np.pi noise = 0.2 ref_freq = 17.77777 # signal = generate_signal(t_vec, A, phi, noise, ref_freq) # lock_in_measurement(signal, t_vec, ref_freq)

import numpy as np
import matplotlib.pyplot as plt

%config InlineBackend.figure_format='retina'

def generate_signal(t_vec, A, phi, noise, freq):
Omega = 2*np.pi*freq
return A * np.sin(Omega*t_vec + phi) + noise * (2*np.random.random)

def lock_in_measurement(signal, t_vec, ref_freq):
Omega = 2*np.pi*ref_freq
ref_0 = 2*np.sin(Omega*t_vec)
ref_1 = 2*np.cos(Omega*t_vec)

signal_0 = signal * ref_0
signal_1 = signal * ref_1

X = np.mean(signal_0)
Y = np.mean(signal_1)

A = np.sqrt(X**2+Y**2)
phi = np.arctan2(Y,X)

print("A=", A, "phi=", phi)

t_vec = np.linspace(0, 0.2, 1001)
A = 1
phi = np.pi
noise = 0.2
ref_freq = 17.77777

signal = generate_signal(t_vec, A, phi, noise, ref_freq)
lock_in_measurement(signal, t_vec, ref_freq)

优化这段pythonimport numpy as np import matplotlib.pyplot as plt %config InlineBackend.figure_format='retina' #输入信号 def generate_signal(t_vec, A, phi, noise, freq): Omega = 2*np.pi*freq noise = 0.2 return A * np.sin(Omega*t_vec + phi) + noise * (2*np.random.random(t_vec.size)-1) #锁相测量 def lock_in_measurement(signal,t_vec,A,phi,noise,ref_freq): Omega = 2*np.pi*ref_freq ref_0 = 2*np.sin(Omega*t_vec) ref_1 = 2*np.cos(Omega*t_vec) #参量 t_vec = np.linspace(0, 0.2, 1001) A = 1 phi = np.pi ref_freq = 17.77777 #混频 signal_0 = signal * ref_0 signal_1 = signal * ref_1 #滤波 X = np.mean(signal_0) Y = np.mean(signal_1) #计算振幅和相位 A = np.sqrt(X**2+Y**2) phi = np.arctan2(Y,X) print("A=", A, "phi=", phi) # 进行锁相放大测量 lock_in_measurement(signal,t_vec,A,phi,noise,ref_freq)

这段代码中存在一些问题，需要进行优化。首先，函数 `generate_signal` 中的变量 `noise` 在函数内部被重新赋值，这可能导致一些问题。其次，在函数 `lock_in_measurement` 中，参量的赋值应该在函数外部进行，而不是在函数内部进行。最后，在函数 `lock_in_measurement` 中，存在递归调用的问题，需要进行修正。以下是优化后的代码：

import numpy as np
import matplotlib.pyplot as plt

%config InlineBackend.figure_format='retina'

def generate_signal(t_vec, A, phi, noise, freq):
    Omega = 2*np.pi*freq
    return A * np.sin(Omega*t_vec + phi) + noise * (2*np.random.random(t_vec.size)-1)

def lock_in_measurement(signal, t_vec, ref_freq):
    Omega = 2*np.pi*ref_freq
    ref_0 = 2*np.sin(Omega*t_vec)
    ref_1 = 2*np.cos(Omega*t_vec)
    
    # 混频
    signal_0 = signal * ref_0
    signal_1 = signal * ref_1
    
    # 滤波
    X = np.mean(signal_0)
    Y = np.mean(signal_1)
    
    # 计算振幅和相位
    A = np.sqrt(X**2+Y**2)
    phi = np.arctan2(Y,X)
    
    print("A=", A, "phi=", phi)
    
# 参量    
t_vec = np.linspace(0, 0.2, 1001)
A = 1
phi = np.pi
noise = 0.2
ref_freq = 17.77777

# 生成信号
signal = generate_signal(t_vec, A, phi, noise, ref_freq)

# 进行锁相放大测量
lock_in_measurement(signal, t_vec, ref_freq)

优化后的代码中，变量 `noise` 的赋值已经移动到了函数外部，函数 `lock_in_measurement` 中的参量赋值也已经移动到了函数外部，递归调用的问题也已经被修正了。此外，函数 `generate_signal` 中的变量 `noise` 已经被正确使用。