matlab对非正弦信号频谱分析,急,y=(40/pi)*sin(100*pi*x)+(40/(3*pi))*sin(300*pi*x)+(40/(5*pi))*sin(500*pi*x)函数已绘出,取的是x=[0:0.00001:0.06](是否需要改?)怎么频谱分析?
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/16 01:17:48
matlab对非正弦信号频谱分析,急,y=(40/pi)*sin(100*pi*x)+(40/(3*pi))*sin(300*pi*x)+(40/(5*pi))*sin(500*pi*x)函数已绘出,取的是x=[0:0.00001:0.06](是否需要改?)怎么频谱分析?
matlab对非正弦信号频谱分析,急,
y=(40/pi)*sin(100*pi*x)+(40/(3*pi))*sin(300*pi*x)+(40/(5*pi))*sin(500*pi*x)
函数已绘出,取的是x=[0:0.00001:0.06](是否需要改?)
怎么频谱分析?
matlab对非正弦信号频谱分析,急,y=(40/pi)*sin(100*pi*x)+(40/(3*pi))*sin(300*pi*x)+(40/(5*pi))*sin(500*pi*x)函数已绘出,取的是x=[0:0.00001:0.06](是否需要改?)怎么频谱分析?
Fs=2048;%采样频率
dx=1/Fs;
N=2048;
df=Fs/N;
x=(0:N-1)*dx;
%从y的表达式可看出信号最高频率为500/2=250Hz,选采样频率至少为2.56倍
y=(40/pi)*sin(100*pi*x)+(40/(3*pi))*sin(300*pi*x)+(40/(5*pi))*sin(500*pi*x);
figure(1)
plot(x(1:200),y(1:200));
figure(2)
Y=abs(fft(y));
f=(0:N/2)*df;
Y(1)=Y(1)/N;
Y(2:N/2)=Y(2:N/2)/(N/2);
Y(N/2+1)=Y(N/2+1)/N;
plot(f,Y(1:N/2+1))