信号的频谱分析实验报告Word文档格式.docx
《信号的频谱分析实验报告Word文档格式.docx》由会员分享,可在线阅读,更多相关《信号的频谱分析实验报告Word文档格式.docx(21页珍藏版)》请在冰豆网上搜索。
clc;
clear;
t=0:
0.0025:
0.5-0.0025;
f=400*t;
w0=100*pi;
y=sin(w0*t);
a=fft(y);
b=abs(a)/200;
d=angle(a)*180/pi;
subplot(311);
plot(t,y);
title('
y=sin(wt)'
);
xlabel('
t'
ylabel('
y(t)'
subplot(312);
stem(f,b);
振幅'
f'
subplot(313);
stem(f,d);
相位'
混叠
0.0115:
0.46-0.0115;
f=(t/0.0115)*2;
b=abs(a)/40;
泄漏
0.5-0.0075;
f=800*t;
b=abs(a)/198;
(2)cos(100*pi*t);
y=cos(w0*t);
y=cos(wt)'
gridon;
holdon;
2.利用FFT,分析并对比方波以及半波对称的正负方波的频谱,改变采样间隔与截断长度,分析混叠与泄漏对信号频谱的影响。
方波正常
0.05:
5-0.05;
f=40*t;
w0=2*pi;
x=square(w0*t)+1;
stem(t,x);
p=fft(x);
q=abs(p)/100;
stem(f,q);
幅频特性曲线'
s=angle(p)*180/pi;
stem(f,s);
相频特性曲线'
方波混叠
closeall
clear
clc
0.011:
0.22-0.011;
f=(t/0.011)*2;
w0=100*pi
x=square(w0*t)
subplot(311)
stem(t,x)
p=fft(x)
q=abs(p)/20
subplot(312)
stem(f,q)
s=angle(p)*180/pi
subplot(313)
stem(t,s)
方波泄漏
5-0.25;
p1=fftshift(p);
q=abs(p1)/96;
半波相对称正常
0.1:
10-0.1;
f=20*t;
x=square(w0*t);
a=0:
9.9;
半波相对称混叠
0.4:
10-0.4;
f=5*t;
q=abs(p)/25;
半波相对称泄漏
10-0.5;
q=abs(p)/95;
3.利用FFT,分析并画出信号的频谱,改变采样间隔与截断长度,分析混叠与泄漏对信号频谱的影响。
e-tu(t)不同Fs的影响
t2=0:
0.2:
20-0.2;
y2=exp(-1*t2);
y22=fft(y2);
e=abs(y22)/10;
y1=exp(-1*t);
y=fft(y1);
a=abs(y)/10;
c=angle(y)*pi/180;
h=angle(y22)*pi/180;
w=2*pi*t;
x=1./(1+j*w);
b=abs(x);
d=angle(x)*pi/180;
figure
(1)
plot(f,a);
plot(f,e,'
g'
plot(f,b,'
r'
幅频特性比较'
figure
(2);
plot(f,c);
plot(f,h,'
plot(f,d,'
相频特性比较'
截断长度L对泄漏的影响
a=abs(y)/500;
b=abs(y)/100;
截断长度对泄漏的影响'
4.利用不同窗函数对内容3.中的信号进行加窗处理,分析对信号频谱的影响;
clearall;
0.01:
5-0.01;
f=200*t;