连续系统的时域、频域分析报告

1、专业资料学生实验报告实验课程:信号与系统EDA实验地点:东1教414学 院 :专 业 :学 号 :姓 名 :word 完美格式专业资料实验题目连续系统的时域、频域分析实验时间实验台号操作成绩报告成绩一、实验目的1学习利用 Matlab 软件实现系统时域输入输出分析法，求解微分方程；2观察和掌握各种常用信号卷积后的结果；3利用 Matlab 实现系统的冲激响应和阶跃响应分析;4利用 Matlab 实现信号的频谱分析二、实验性质验证性实验三、预习内容系统时域输入 - 输出分析方法；信号卷积；系统的阶跃响应和冲激响应的概念；信号的时域和频域的关系。四、实验内容1系统时域输入输出分析法某系统的输入输出

2、描述方程为：dy (t )1 y (t)x (t)Matlab 软件实现系dt2，给出运行结果，理解统时域输入 - 输出分析方法，说明函数lsim()用法。程序：t=0:10;%确定信号时间范围x=ones(1,length(t); %定义输入信号形式b=1;%方程描述a=1 0.5;s=lsim(b,a,x,t);%方程求解plot(t,s,'y-')%系统输出信号波形绘制lsim 函数的用法：对于线性时不变系统，任意输入，输出时间响应。系统方程为状态方程时，可以得到状态轨迹。word 完美格式专业资料2信号卷积，根据 PPT中的实验 2.2 和 2.3 内容完成课堂练习，

3、写出程序及运行结果。用 Matlab 实现卷积运算f (t) * h(t ) ，其中 f (t )2 (t )(t2), h(t )e t (t ) ，h2(t )h( 2t ) ；对比说明信号f (t ) 分别输入系统 h(t )和 h2(t ) 时的输出有什么区别并分析原因。>> p=0.01;nf=0:p:4;f=2*(heaviside(nf)-heaviside(nf-2);nh=0:p:6;h=exp(-nh).*(nh>0);y=conv(f,h);t=0:length(y)-1;subplot(3,1,1),stairs(nf,f);title('f(

4、t)');axis(0 6 0 2.1);subplot(3,1,2),plot(nh,h);title('h(t)');axis(0 6 0 1.1);subplot(3,1,3),plot(0.01*t,y); title('y(t)=f(t)*h(t)');word 完美格式专业资料>> p=0.01;nf=0:p:4; f=2*(heaviside(nf)-heaviside(nf-2);nh=0:p:6;h=exp(-2*nh).*(2*nh>0);y=conv(f,h);t=0:length(y)-1;subplot(3,1

5、,1),stairs(nf,f);title('f(t)');axis(0 6 0 2.1);subplot(3,1,2),plot(nh,h);title('h(t)');axis(0 6 0 1.1);subplot(3,1,3),plot(0.01*t,y); title('y(t)=f(t)*h(t)');word 完美格式专业资料区别： h（t ）横轴缩短一半y（ t ）纵轴横轴缩短一半原因： t 扩大 2 倍 横轴缩短其卷积缩小到原来的 4倍 故纵轴缩小2 倍3系统的冲激响应和阶跃响应分析uc'' ( t)R u

6、9; (t)1uc ( t)1(t )4 ,L 1H,C1/ 3F（2）R 2 ,L 1H,C1F（3）已知二阶系统方程LLCLC，（1）RR 1 ,L 1H,C 1F（4）R 0, L1H ,C1 F ，根据不同情况下的实验结果分析系统参数变化时系统输出有什么变化规律。程序：R=input('电阻 R='); %以交互方式输入电阻R 的值L=input('电感 L='); %以交互方式输入电阻L 的值C=input('电容 C='); %以交互方式输入电阻C 的值b=1/(L*C);a=1 R/L 1/(L*C);impulse(b,a);wo

7、rd 完美格式专业资料word 完美格式专业资料分析：电阻增大 峰值降低 曲线坡度变缓 电容变大 峰值变小 坡度变缓 电阻为 0 发生错误4. 信号的频谱分析：根据PPT中程序 2.6 和 2.7 完成课堂作业，写出程序及运行结果。word 完美格式专业资料f (t )sin( t )sin( 2t )Matlabt( )用分析信号以及的频谱。根据结果分析信号时域变f t2t化与信号频谱变化间关系。关系；时域与频域成反比? m32.m：?syms t w f ft;%定义符号变量?f=sin(t)/t;%信号?ft=f*exp(-j*w*t);%计算被积函数?F=int(ft,t,-2,2);

8、%计算傅立叶变换 F(w)?F=simple(F);F%化简?subplot(2,1,1),ezplot(f,-2 2); %绘制信号?axis(-3 3 0 1.1);title('sin(t)/t信号 ');?subplot(2,1,2),ezplot(abs(F),-8:0.01:8);%绘制信号的频谱?title('sin(t)/t信号的频谱 ');word 完美格式专业资料? m32.m：?syms t w f ft;%定义符号变量?f=sin(2*t)/(2*t);%信号?ft=f*exp(-j*w*t);%计算被积函数?F=int(ft,t,-2,2);%计算傅立叶变换 F(w)?F=simple(F);F%化简?subplot(2,1,1),ezplot(f,-2 2);%绘制信号?axis(-3 3 0 1.1);title(&#