发布时间 : 星期日 文章信号与系统matlab实验傅里叶分析及应用报告答案更新完毕开始阅读
w
f = sym('Heaviside(t+1)-Heaviside(t-1)'); Fw = simplify(fourier(f)); F = Fw.*Fw;
ezplot(abs(F)),grid on title('Fw*Fw');
triangle = sym('(t+2)*Heaviside(t+2)-2*t*Heaviside(t)+(t-2)*Heaviside(t-2)'); Fttriangle = fourier (triangle); F = simplify (Fttriangle); ezplot(abs(F),[-6,6]),grid on; title('triangle FT');
5、设有两个不同频率的余弦信号,频率分别为 仏=100Hz , f^ 3800Hz ;现在使用 抽样频率fs = 4000Hz对这三个信号进行抽样,使用 MATLAB命令画出各抽样信号的波形 和频谱,并分析其频率混叠现象[建议:抽样信号的频谱图横坐标范围不小于 -10000Hz~10000Hz或-20000*pi~20000*pi rad/s]。
F1=100HzM,
实验代码如下: Ts=14000;
5
dt=0.0001;
t1=-0.006:dt:0.006; ft=cos(2*pi*100*t1); subplot(221);
plot(t1,ft),grid on; axis([-0.006 0.006 -1.6 1.6]); xlabel('Time(sec)'),ylabel('f(t)') title('余弦信号波形');N=5000; k=-N:N;
W=2*pi*k/((2*N+1)*dt);
Fw=dt*ft*exp(-j*t1'*W); subplot(222); plot(W,abs(Fw)); grid on;
axis([-20000 20000 0 0.006]); xlabel('\\omega'),ylabel('F(w)'); title(' 余弦信号的频谱t2=-0.006:Ts:0.006; fst=cos(2.*pi.*100*t2); subplot(223)
plot(t1,ft,':'),hold on stem(t2,fst),grid on axis([-0.006 0.006 -1.5 1.5]) xlabel('Time (sec)'),ylabel('fs(t)') title(' 抽样后的信号 '),hold off Fsw=Ts*fst*exp(-j*t2'*W); subplot(224) plot(W,abs(Fsw)),grid on axis([-20000 20000 0 0.006]) xlabel('\\omega'),ylabel('Fsw') title(' 抽样信号的频谱 ')
6
');
余弦信号液形
6
買1屮余弦信号的频谱
4
Ll_
2
-6 0 5 -2-1012
Time(sec) x 让 抽样
后的信号
北1代抽样信号的频谱
i ■ i i i ■ i 1A ___________ ■ ___ __ 1 4 r I i ■ ■丄 ~ —八 丁 八-=— 1 > 1 1 A | 0 Time (sec)
5 x 10'3
u -1 0 ILB 1 2
F2=3800Hz时,
实验代码如下:
Ts=14000; dt=0.0001;
t仁-0.006:dt:0.006; ft=cos(2*pi*3800*t1); subplot(221);
plot(t1,ft),grid on; axis([-0.006 0.006 -1.6 1.6]); xlabel('Time(sec)'),ylabel('f(t)') title('余弦信号波形');N=5000; k=-N:N; W=2*pi*k/((2*N+1)*dt);
Fw=dt*ft*exp(-j*t1'*W); subplot(222); plot(W,abs(Fw)); grid on;
axis([-20000 20000 0 0.006]); xlabel('\\omega'),ylabel('F(w)'); title('余弦信号的频谱'); t2=-0.006:Ts:0.006; fst=cos(2.*pi.*100*t2);
7
subplot(223)
plot(t1,ft,':'),hold on stem(t2,fst),grid on
axis([-0.006 0.006 -1.5 1.5]) xlabel('Time (sec)'),ylabel('fs (t)') title('抽样后的信号'),hold off Fsw=Ts*fst*exp(-j*t2'*W);
subplot(224) plot(W,abs(Fsw)),grid on axis([-20000 20000 0 0.006]) xlabel('\\omega'),ylabel('Fsw') title('抽样信号的频谱')
余弦信号液形
Time(sec)
抽样启的信号
6
1屮抽样信号的频谱
\\ I \\
I
I
||
LJL
2
-10 12
①
x 104
(六)结合抽样定理,利用MATLAB编程实现Sa(t)信号经过冲激脉冲抽样后得到的抽样信 号fs t及其频谱[建议:冲激脉冲的周期分别取 4*pi/3 s、pi s、2*pi/3 s三种情况对比],并 利用fs t构建Sa(t)信号。
周期取4*pi/ 3 s:
8