信号分析与处理MATLAB仿真实验报告Word文件下载.docx
《信号分析与处理MATLAB仿真实验报告Word文件下载.docx》由会员分享，可在线阅读，更多相关《信号分析与处理MATLAB仿真实验报告Word文件下载.docx（12页珍藏版）》请在冰点文库上搜索。
LTI系统的单位脉冲响应为h（n）=（0.9）
ε（n）,输入序列为x（n）=ε（n）-ε（n-10）,求系统的输出y（n）。
程序
function[y,ny]=conv_m（x,nx,h,nh）;
y=conv（x（5:
45）,h）;
ny=nh;
y=y（1:
41）;
subplot（3,1,1）;
stem（nx,x,'
filled'
）;
axis（[-4,40,0,1]）;
title（'
x[n]'
subplot（3,1,2）;
stem（nh,h,'
h[n]'
subplot（3,1,3）;
stem（ny,y,'
axis（[-4,40,0,8]）;
y[n]'
实验结果：
结果分析：
改变参数以后的程序：
axis（[-4,40,0,2]）;
输出结果：
实验8连续时间信号的卷积
学会用MATLAB实现对连续时间信号的卷积，掌握利用卷积的方法求系统零状态响应，并与离散系统比较。
使用MATLAB中求卷积函数的conv（），并对结果分析总结。
已知某连续系统的单位冲激响应h（t）=e-5t*ε（t）,输入信号x（t）=ε（t）-ε（t-1）,求解系统的零状态响应y（t）
p=0.01;
t=0:
p:
2;
x=heaviside（t）-heaviside（t-1）;
x
（1）=1;
x（101）=1;
h=exp（-5*t）;
hl=h*p;
y=conv（x,hl）;
length（t））;
plot（t,x）;
axis（[-1,2,0,1.2]）;
x（t）'
plot（t,h）;
axis（[0,2,0,1.2]）;
h（t）'
plot（t,y）;
axis（[0,2,0,0.22]）;
y（t）'
实验结果
结果分析
改变参数后程序
axis（[-1,2,0,1.0]）;
axis（[0,2,0,1.0]）;
axis（[0,2,0,0.5]）;
输出图形
实验18周期信号的频谱分析
学会用MATLAB绘制周期信号的频谱
使用MATLAB中的符号积分函数int（）
对于周期矩形脉冲信号，若A=1，T=5s,г=T/4=1.25s,试求其傅里叶系数Xk，并绘制出频谱图。
symstk;
T=5;
tao=1.25;
A=1;
x0=int（A,t,-tao/2,tao/2）/T
x=A*exp（-j*k*2*pi/T*t）;
xk=int（x,t,-tao/2,tao/2）/T;
xk=simple（xk）
t=[-2*T:
0.01:
2*T];
x1=rectpuls（t,tao）;
subplot（1,2,1）;
plot（t,x1）
holdon
x1=rectpuls（t-5,tao）;
x1=rectpuls（t+5,tao）;
holdoff
周期矩形脉冲（tao=T/4）'
）
xlabel（'
t'
axis（[-8,8,0,1.2]）
k=[-20:
-1,eps,1:
20];
xk=subs（xk,k,'
k'
subplot（1,2,2）;
stem（k,xk,'
line（[-20,20],[0,0]）
周期矩形脉冲的频谱'
ylabel（'
Xk'
减小脉冲的宽度，令tao=T/4=2.5
则输出图形为
实验21非周期信号的傅里叶变换
学会用MATLAB求解非周期信号的傅里叶变换
使用MATLAB中直接求傅里叶变换域逆变换的函数fourier（）与ifourier（）
已知矩形脉冲信号x（t）=ε（t+1）-ε（t-1），将其乘以载波信号cos（10∏t）,得到已调信号y（t）=x（t）cos（10∏t）,试用MATLAB画出信号x（t）,y（t）的波形及其频谱，并观察频谱搬移情况。
程序：
R=0.005;
t=-1.2:
R:
1.2;
x=Heaviside（t+1）-Heaviside（t-1）;
x（41）=1;
x（441）=1;
y=x.*cos（10*pi*t）;
subplot（2,2,1）;
plot（t,x）
axis（[-2,2,0,1.2]）
subplot（2,2,3）;
plot（t,y）
y（t）=x（t）*cos（10*pi*t）'
W1=40;
N=1000;
k=-N:
N;
W=k*W1/N;
X=x*exp（-j*t'
*W）*R;
X=real（X）;
Y=y*exp（-j*t'
Y=real（Y）;
subplot（2,2,2）;
plot（W,X）
w'
X（w）'
subplot（2,2,4）;
plot（W,Y）
Y（w）'
输出结果
修改采样点数后的输出结果W1=25
实验22离散时间信号的频谱分析
学会用MATLAB求解离散时间福利叶变换
已知序列x（n）=（-0.9）
-5≤n≤5,求其离散时间傅里叶变换X（ejΩ）
n=-5:
5;
x=（-0.9）.^n;
k=-200:
200;
w=（pi/100）*k;
X=x*（exp（-j*pi/100））.^（n'
*k）;
magX=abs（X）;
angX=angle（X）;
subplot（2,1,1）;
plot（w/pi,magX）;
axis（[-2,2,0,15]）
\Omega/（\pi）'
幅度'
subplot（2,1,2）;
plot（w/pi,angX）/pi;
axis（[-2,2,-4,4]）
相位/（\pi）'
改变参数后的程序
axis（[-2,1,0,5]）
axis（[-2,1,-4,4]）
升级会员 