哈工大威海信号系统实验报告完整版
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
《信号与系统》实验报告
实验一 典型连续时间信号描述及运算 实验报告要求:
(1)仿照单边指数信号的示例程序,按要求完成三种典型连续信号,即:正弦信号、衰减正弦信号、钟型信号的波形绘制。(要求:要附上程序代码,以下均如此,不再说明) (2)根据《信号与系统》教材第一章的习题1.1(1,3,5,8)函数形式绘制波形。 (3)完成三种奇异信号,即:符号函数、阶跃信号、单位冲激信号的波形绘制。 (4)完成实验一中信号的运算:三、6 实验内容中的 (1)(2)(3)(4)。 (5)求解信号的直流/交流分量,按第四部分的要求完成。 正文: (1)
<1>正弦信号:
代码:>> 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:1:250;
>> f1=400*exp(-1.*t.*t./10000); >> f1=400*exp(-1.*t.*t./22500); >> f1=400*exp(-1.*t.*t./62500);
姓 名:
学 号:
同组人:
无
指导教师:
成 绩:
>> plot(t,f1,'-',t,f2,'--',t,f3,'-.')
(2)习题1,3,5,8
<1>
代码:t=0:1:10;
f=t;
plot(t,f)
<3>
代码:t=1:1:10;
f=t;
plot(t,f)
<5>
代码:t=0:1:10;
f=2-exp(-1.*t.);
plot(t,f)
<8>
代码:t=1:0.1:2;
f=exp(-1.*t.)*cos(10*pi*t);
plot(t,f)
(3)三种奇异函数
<1>符号函数
代码: t=-5:0.05:5;
f=sign(t);
plot(t,f)
<2>阶跃信号
代码:>>
t=-5:0.1:5;
>>
f=u(t);
>> plot(t,f)
<3>单位冲激信号
代码:function chongji(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
<1>
syms t
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);
function f=u(t) f=(t>0);
<2>
4、
function f=u(t)
f=(t>0)
syms t
f1=sym('(-t+4)*(u(t)-u(t-4))'); subplot(1,3,1);
ezplot(f1);
f2=sym('sin(2*pi*t)');
subplot(1,3,2);
ezplot(f2);
f6=f1.*f2;
subplot(1,3,3);
ezplot(f6);
5、
function f=u(t)
f=(t>0)
syms t
f1=sym('(-t+4)*(u(t)-u(t-4))'); f2=sym('sin(2*pi*t)');
subplot(1,3,1);
ezplot(f2);
f6=f1.*f2;
y6=subs(f6,t,t-2);
subplot(1,3,2);
ezplot(y6);
f7=y6+f2;
subplot(1,3,3);
ezplot(f7);
四、
t=0:0.1:500;
f=100.*abs(sin(2.*pi.*t./50)); plot(t,f,t,fD,t,fA)
调用子程序:
function fD=fDC(f)
fD=mean(f);
function fA=fAC(f,fD)
fA=f-fD;
(5)求解信号的交直流分量
代码:function fD=fDC(f)
fD=mean(f);
function fA=fAC(f,fD)
fA=f-fD;
t=0:0.1:500;
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)+3 e(t) ,e(t)=u(t) 。
要求:1) 用解析法求系统的零状态响应r(t );(手写)
2)利用MATLAB绘出系统零状态响应的时域仿真波形,并验证1)的结果是否正确,附程序代码和波形;
3)利用MATLAB绘出系统的冲激响应和阶跃响应波形,附程序代码和波形。
正文:
(1)求解下面两个信号的卷积积分。
<1>在实验报告中推导出这两个信号卷积积分运算表达式;(手写)
<2>利用MATLAB 进行求解验证,附程序代码和波形。
代码:p=0.01;
k1=-1/2:p:1;
f1=1;
k2=0:p:2;
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:p:(k3*p+k0);
subplot(2,2,1)
plot(k1,f1)
xlabel('t')
ylabel('f1(t)')
subplot(2,2,2)
plot(k2,f2)
ylabel('f2(t)')
subplot(2,2,3)
plot(k,f);
h=get(gca,'position');
h(3)=2.5*h(3);