哈工大(威海)信号系统实验报告完整版Word文档下载推荐.docx
《哈工大(威海)信号系统实验报告完整版Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《哈工大(威海)信号系统实验报告完整版Word文档下载推荐.docx（27页珍藏版）》请在冰豆网上搜索。
代码：
>
t=-250:
1:
250;
f1=150*sin（2*pi*t/100）;
f2=150*sin（2*pi*t/200）;
f3=150*sin（2*pi*t/200+pi/5）;
plot（t,f1,'
-'
t,f2,'
--'
t,f3,'
-.'
）
2>
衰减正弦信号
3>
t=-250:
f1=400*exp（-1.*t.*t./10000）;
f1=400*exp（-1.*t.*t./22500）;
f1=400*exp（-1.*t.*t./62500）;
（2）习题1，3，5，8
代码：
t=0:
10;
f=t;
plot（t,f）
t=1:
5>
f=2-exp（-1.*t.）;
8>
0.1:
2;
f=exp（-1.*t.）*cos（10*pi*t）;
（3）三种奇异函数
符号函数
代码:
t=-5:
0.05:
5;
f=sign（t）;
plot（t,f）
阶跃信号
代码:
>
t=-5:
f=u（t）;
plot（t,f）
单位冲激信号
functionchongji（t1,t2,t0）
dt=0.01;
t=t1:
dt:
t2;
n=length（t）;
x=zeros（1,n）;
x（1,（-t0-t1）/dt+1）=1/dt;
stairs（t,x）;
axis（[t1,t2,0,1.2/dt]）title（'
δ（t）'
（4）实验三1234
symst
f1=sym（'
（-t+4）*（u（t）-u（t-4））'
）;
subplot（1,2,1）;
ezplot（f1）;
y1=subs（f1,t,-t）;
f3=f1+y1;
subplot（1,2,2）;
ezplot（f3）;
functionf=u（t）f=（t>
0）;
4、
functionf=u（t）
f=（t>
0）
symst
subplot（1,3,1）;
f2=sym（'
sin（2*pi*t）'
subplot（1,3,2）;
ezplot（f2）;
f6=f1.*f2;
subplot（1,3,3）;
ezplot（f6）;
5、
y6=subs（f6,t,t-2）;
ezplot（y6）;
f7=y6+f2;
ezplot（f7）;
四、
500;
f=100.*abs（sin（2.*pi.*t./50））;
plot（t,f,t,fD,t,fA）
调用子程序：
functionfD=fDC（f）
fD=mean（f）;
functionfA=fAC（f,fD）
fA=f-fD;
（5）求解信号的交直流分量
fD=mean（f）;
fA=f-fD;
f（t）=100|sin（2*PI*t/50）|;
plot（t,fD,t,fA）
实验二线性系统时域分析
（1）求解下面两个信号的卷积积分。
要求：
1）在实验报告中推导出这两个信号卷积积分运算表达式；
（手写）
2）利用MATLAB进行求解验证，附程序代码和波形。
（2）已知描述系统的微分方程和激励信号如下
r¢
¢
（t）+3r¢
（t）+2r（t）=e¢
（t）+3e（t），e（t）=u（t）。
1）用解析法求系统的零状态响应r（t）；
2）利用MATLAB绘出系统零状态响应的时域仿真波形，并验证1）的结果是否正确，附程序代码和波形；
3）利用MATLAB绘出系统的冲激响应和阶跃响应波形，附程序代码和波形。
在实验报告中推导出这两个信号卷积积分运算表达式；
利用MATLAB进行求解验证，附程序代码和波形。
p=0.01;
k1=-1/2:
p:
1;
f1=1;
k2=0:
f2=0.5*k2;
[f,k]=sconv（f1,f2,k1,k2,p）
function[f,k]=sconv（f1,f2,k1,k2,p）
f=conv（f1,f2）;
f=f*p;
k0=k1
（1）+k2
（1）;
k3=length（f1）+length（f2）-2;
k=k0:
（k3*p+k0）;
subplot（2,2,1）
plot（k1,f1）
xlabel（'
t'
）
ylabel（'
f1（t）'
subplot（2,2,2）
plot（k2,f2）
f2（t）'
subplot（2,2,3）
plot（k,f）;
h=get（gca,'
position'
h（3）=2.5*h（3）;
set（gca,'
h）
title（'
f（t）=f1（t）*f2（t）'
xlabel（'
f（t）'
用解析法求系统的零状态响应r（t）；
利用MATLAB绘出系统零状态响应的时域仿真波形，并验证<
的结果是否正确，附程序代码和波形；
a=[1,3,2];
b=[0,1,3];
impulse（a,b）
利用MATLAB绘出系统的冲激响应和阶跃响应波形，附程序代码和波形。
a=[1,3,2];
b=[0,1,3];
step（b,a）
代码；
p=0.01;
t=0:
x=exp（-t）;
lsim（b,a,x,t）
实验三：
实验报告内容：
对所给音频信号，进行时域压缩和扩展，画出时域波形与幅度谱，使其满足以下要求。
（1）将music1.wav的音调变低a倍（0.8<
a<
1），及变高a倍（1<
1.3），比较变换前后的10%带宽（最大值10%处的带宽）的变化情况。
（2）将music2.wav的音调变化到与儿童和男声相似的声音，分析变换前后的10%带宽的变化情况（变宽或变窄多少）。
①a=[0.91.01.2];
a1=30;
b=11025;
[x0b]=wavread（'
F:
\信号与系统\music1.wav'
form=1:
3;
a2=a1*a（m）;
x=resample（x0,a1,a2）;
wavplay（x,b）;
1/b:
（length（x）-1）/b;
subplot（3,3,3*m-2）;
plot（t,x）;
axis（[06-infinf]）;
title（['
a='
num2str（a（m））,'
的时域图'
]）;
N=（length（resample（x0,a1,a1*max（a））））;
f=（-N/2:
N/2-1）/N*b;
X=fftshift（fft（x,N））;
AMP=abs（X）;
AMP0=AMP/max（AMP）;
subplot（3,3,3*m-1）;
plot（f,AMP0）;
的频域图'
subplot（3,3,3*m）;
axis（[-infinf00.1]）
的10%带宽图'
end
②
a=[10.71.4];
str=['
原'
'
男'
‘童'
];
\信号实验_2018\music2.wav'
a2=a1*a（m）;
x=resample（x0,a1,a2）;
wavplay（x,b）;
t=0:
subplot（3,3,3*m-2）;
plot（t,x）;
axis（[03-infinf]）
升级会员 