信号处理上机实验结果Word格式.docx
《信号处理上机实验结果Word格式.docx》由会员分享，可在线阅读，更多相关《信号处理上机实验结果Word格式.docx（16页珍藏版）》请在冰豆网上搜索。
3;
ft=A*sin（w0*t+phi）;
axis（[03-1.11.1]）;
（3）抽样信号
x=linspace（-20,20）;
y=sinc（x/pi）;
plot（x,y）
axis（[-2121-0.51.1]）;
（4）矩形脉冲
0.001:
4;
T=1;
ft=rectpuls（t-2*T,2*T）;
axis（[-15-0.11.1]）;
（5）阶跃信号
functionf=Heaviside（t）
f=（t>
0）
------------
t=-1:
ft=Heaviside（t）;
axis（[-13-0.1
（6）复指数信号的时域波形
0.1:
60;
f=exp（-0.1*t）.*sin（2/3*t）;
plot（t,f）;
ylabel（'
f（t）'
）;
xlabel（'
Time（sec）'
axis（[060-11]）;
（7）加入随机噪声的正弦波
t=（0:
50）;
y=sin（2*pi*50*t）;
s=y+randn（size（t））;
Subplot（2,1,1）;
plot（t（1:
100）,y（1:
100））;
Subplot（2,1,2）;
100）,s（1:
（8）周期矩形波
t=（0:
0.0001:
5）;
y=A*square（2*pi*t,20）;
plot（t,y）;
axis（[05-1.51.5]）
2信号的基本运算
设
求出
symst;
f=sym（'
（t/2+1）*（heaviside（t+2）-heaviside（t-2））'
subplot（2,3,1）,ezplot（f,[-3,3]）;
title（'
y1=subs（f,t,t+2）;
subplot（2,3,2）,ezplot（y1,[-5,1]）;
f（t+2）'
y2=subs（f,t,t-2）;
subplot（2,3,3）,ezplot（y2,[-1,5]）;
f（t-2）'
y3=subs（f,t,-t）;
subplot（2,3,4）,ezplot（y3,[-3,3]）;
f（-2t）'
y4=subs（f,t,2*t）;
实验二
例1求系统y”（t）+2y’（t）+100y（t）=10f（t）的零状态响应，已知f（t）=（sin2πt）δ（t）。
ts=0;
te=5;
dt=0.01;
sys=tf（[10],[12100]）;
t=ts:
dt:
te;
f=10*sin（2*pi*t）;
y=lsim（sys,f,t）;
）
y（t）'
）
例2求系统y”（t）+2y’（t）+100y（t）=10f（t）的零状态响应，已知f（t）=（t）。
y=impulse（sys,t）;
h（t）'
实验三
[例1]试用MATLAB画出图示周期三角波信号的频谱
解：
周期信号的频谱为
N=8;
n1=-N:
-1;
%计算n=-N到-1的Fourier系数
c1=-4*j*sin（n1*pi/2）/pi^2./n1.^2;
c0=0;
%计算n=0时的Fourier系数
n2=1:
N;
%计算n=1到N的Fourier系数
c2=-4*j*sin（n2*pi/2）/pi^2./n2.^2;
cn=[c1c0c2];
n=-N:
subplot（2,1,1）;
stem（n,abs（cn））;
Cn的幅度'
subplot（2,1,2）;
stem（n,angle（cn））;
Cn的相位'
\omega/\omega0'
[例2]求周期矩形脉冲的Fourier级数表示式。
并用MATLAB求出由前N项Fourier级数系数得出的信号近似波形。
取A=1,T=2,t=1,w0=p
•t=-2:
2;
•N=input（'
Numberofharmonics='
•c0=0.5;
•xN=c0*ones（1,length（t））;
%dccomponent
•forn=1:
2:
N
•xN=xN+cos（pi*n*t）*sinc（n/2）;
•end
•plot（t,xN）;
•title（['
N='
num2str（N）]）;
N=5N=15
N=50N=500
[例]试用数值方法近似计算三角波信号的频谱
解:
图示三角波可表示为
三角波信号频谱的理论值为
F（jw）=Sa2（w/2）
w=linspace（-6*pi,6*pi,512）;
N=length（w）;
F=zeros（1,N）;
fork=1:
N
F（k）=quad（'
sf1'
-1,1,[],[],w（k））;
end
figure
（1）;
>
plot（w,real（F））;
F=zeros（1,N）;
title（'
'
xlabel（'
\omega'
ylabel（'
升级会员 