matlab中求下边带频谱函数出错无图像t=1000e-6; % 信号的时间长度,单位sfs=60e6; % 采样频率,单位Hzts=1/fs;fc=17.5e6; % 线性调频信号的起始频率f0=20e6; % 本振信号的中心频率B=5e6; % 调频带宽ft=0:1/fs*200:t

来源:学生作业帮助网 编辑:作业帮 时间:2024/04/28 04:05:56
matlab中求下边带频谱函数出错无图像t=1000e-6; % 信号的时间长度,单位sfs=60e6; % 采样频率,单位Hzts=1/fs;fc=17.5e6; % 线性调频信号的起始频率f0=20e6; % 本振信号的中心频率B=5e6; % 调频带宽ft=0:1/fs*200:t

matlab中求下边带频谱函数出错无图像t=1000e-6; % 信号的时间长度,单位sfs=60e6; % 采样频率,单位Hzts=1/fs;fc=17.5e6; % 线性调频信号的起始频率f0=20e6; % 本振信号的中心频率B=5e6; % 调频带宽ft=0:1/fs*200:t
matlab中求下边带频谱函数出错无图像
t=1000e-6; % 信号的时间长度,单位s
fs=60e6; % 采样频率,单位Hz
ts=1/fs;
fc=17.5e6; % 线性调频信号的起始频率
f0=20e6; % 本振信号的中心频率
B=5e6; % 调频带宽
ft=0:1/fs*200:t-1/fs; % 时间序列步长
N=length(ft);
k=B/fs*2*pi/max(ft); % modulate 函数产生
y=modulate(ft,fc,fs,'fm',k);
y_fft_result=fft(y);
figure,subplot(2,1,1),plot(ft,y),xlabel('t(单位:秒)'),ylabel('y(单位:伏)'),title('线性调频信号y(t)');
subplot(2,1,2),plot((0:fs/N:fs/2-fs/N),10*log(abs(y_fft_result(1:N/2)))),xlabel('频率f(单位:Hz)'),title('线性调频信号y(t)的频谱');
figure(3);
c1=cos(2*pi*f0*ft);
c2=sin(2*pi*f0*ft);
subplot(3,2,3);
u1=y(1:N).*c1(1:N)+imag(hilbert(y(1:N))).*c2(1:N);
plot(ft,u1);
axis([0 1e-3 -1 1]);
title('下边带调制信号');
xlabel('t/s');
ylabel('幅度/v');
grid on;
y2=fft(u1);
mag2=(abs(y2));
subplot(3,2,4);
plot((0:fs/N:fs/2-fs/N),mag2);
title('下边带频域信号');
xlabel('f/KHz');
ylabel('幅度/v');grid on;
报错信息为:
Error using ==> plot
Vectors must be the same lengths.
Error in ==>
plot((0:fs/N:fs/2-fs/N),mag2);
本人是新手请见谅.

matlab中求下边带频谱函数出错无图像t=1000e-6; % 信号的时间长度,单位sfs=60e6; % 采样频率,单位Hzts=1/fs;fc=17.5e6; % 线性调频信号的起始频率f0=20e6; % 本振信号的中心频率B=5e6; % 调频带宽ft=0:1/fs*200:t
clear
t=1000e-6; % 信号的时间长度, 单位s
fs=60e6; % 采样频率,单位Hz
ts=1/fs;
fc=17.5e6; % 线性调频信号的起始频率
f0=20e6; % 本振信号的中心频率
B=5e6; % 调频带宽
ft=0:1/fs*200:t-1/fs; % 时间序列步长
N=length(ft);
k=B/fs*2*pi/max(ft); % modulate 函数产生
y=modulate(ft,fc,fs,'fm',k);
y_fft_result=fft(y);
figure,subplot(2,1,1),plot(ft,y),xlabel('t(单位:秒)'),ylabel('y(单位:伏)'),title('线性调频信号y(t)');
subplot(2,1,2),plot((0:fs/N:fs/2-fs/N),10*log(abs(y_fft_result(1:N/2)))),xlabel('频率f(单位:Hz)'),title('线性调频信号y(t)的频谱');
figure(3);
c1=cos(2*pi*f0*ft);
c2=sin(2*pi*f0*ft);
subplot(3,2,3);
u1=y(1:N).*c1(1:N)+imag(hilbert(y(1:N))).*c2(1:N);
plot(ft,u1);
axis([0 1e-3 -1 1]);
title('下边带调制信号');
xlabel('t/s');
ylabel('幅度/v');
grid on;
y2=fft(u1);
mag2=(abs(y2));
subplot(3,2,4);
gg = 0:fs/N:fs/2-fs/N;
plot(gg,mag2(1:150));
title('下边带频域信号');
xlabel('f/KHz');
ylabel('幅度/v');grid on;