杭电信号与系统实验
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第十一章：
1.
clear all;
close all;
t=-10:0.01:10;
ft=exp(abs(-0.4*t));
plot(t,ft)
2
clear all;
close all;
h=0.001
t=-3:h:3;
%ft的函数表达式%
ft=0.5*exp(-2*t).*stepfun(t,0);
%ft1的函数表达式%
ft1=0.5*exp(-2*(t-1)).*stepfun(t-1,0); %ft的导数表达式%
y1=diff(ft)*1/h;
%产生ft的波形%
subplot(3,1,1)
plot(t,ft);
title('f(t）波形');
%产生ft1的波形%
subplot(3,1,2)
plot(t,ft1);
title('f(t-1)的波形');
%产生ft的微分波形%
subplot(3,1,3)
plot(t(1:length(t)-1),y1);
title('f(t）的微分波形');
第十二章：
1). clear all;
close all;
f=fourier(cos(2*t).*sin(3*t)); ezplot(f);
2). clear all;
close all;
r=0.02;
t=-5:r:5;
N=200;
w=2*pi;
k=-N:N;
w=k*w/N;
f1=sinc(t/pi);
F=r*f1*exp(-j*t'*w);
F1=abs(F);
P1=angle(F);
subplot(3,1,1);
plot(t,f1);
grid on;
xlabel('t');
ylabel('f(t)');
title('f(t)');
subplot(3,1,2)
plot(w,F1);
xlabel('w');
grid on;
ylabel('F(jw)');
subplot(3,1,3);
plot(w,P1*180/pi);
grid;
xlabel('w');
ylabel('相位(度)');
3). clear all;
close all;
syms t w;
ifourier((((sin(w/4))/(w/4))^2),t);ans 4).
clear all;
close all;
r=0.02;
t=-5:r:5;
w=2*pi;
k=-N:N;
w=k*w/N;
f1=sin(2*pi*3*t+eps)./(pi*3*t+eps); F=r*f1*exp(-j*t'*w);
F1=abs(F);
P1=angle(F);
subplot(3,1,1);
plot(t,f1);
grid on;
xlabel('t');
ylabel('f(t)');
title('f(t)');
subplot(3,1,2)
plot(w,F1);
xlabel('w');
grid on;
ylabel('F(jw)');
subplot(3,1,3);
plot(w,P1*180/pi);
grid;
xlabel('w');
ylabel('相位(度)');
clear all;
close all;
r=0.02;
t=-5:r:5;
N=200;
w=2*pi;
k=-N:N;
w=k*w/N;
f1=sin(2*pi*(t-2)+eps)./(pi*(t-2)+eps); F=r*f1*exp(-j*t'*w);
F1=abs(F);
P1=angle(F);
subplot(3,1,1);
plot(t,f1);
grid on;
xlabel('t');
ylabel('f(t)');
title('f(t)');
subplot(3,1,2)
plot(w,F1);
xlabel('w');
grid on;
ylabel('F(jw)');
subplot(3,1,3);
plot(w,P1*180/pi);
grid;
xlabel('w');
ylabel('相位(度)');
clear all;
close all;
r=0.02;
t=-5:r:5;
N=200;
w=2*pi;
k=-N:N;
w=k*w/N;
f1=sin(2*pi*(t/3)+eps)./(pi*(t/3)+eps); F=r*f1*exp(-j*t'*w);
F1=abs(F);
P1=angle(F);
subplot(3,1,1);
plot(t,f1);
grid on;
xlabel('t');
ylabel('f(t)');
title('f(t)');
subplot(3,1,2)
plot(w,F1);
xlabel('w');
grid on;
ylabel('F(jw)');
subplot(3,1,3);
plot(w,P1*180/pi);
grid;
xlabel('w');