实验三 信号与系统的频域分析

信号与系统实验报告信号与系统实验报告4-4-PAGE4姓名：姓名：实验名称：实验三连续系统的频率响应分析学号：实验时间：评语：一.实验目的深入理解系统频率响应的物理意义Matlab分析系统频率响应的方法理解系统对信号的作用关系二.实验内容RC电路如图：R+y(t)f(t) C-对不同RC值，freqs画出系统的幅度响应曲线H()。f(t)=cos(100t)+RCf(t)f(t)y(t)在0-0.2s范围内的波形。（3）50Hz的交流信号经过全波整流后可表示为f(t)=10t)试取不同的RC值，计算并画出f(t)通过系统的响应y(t)。计算f(t)的直流分量。f(t)和f(t)见下图：1 2f(tf(t)11011f2(t)t012（1）取t:0.05:2.5,计算信号f(t)=f(t)+f

(t)*cos(50t)的值，并画出波形;1 2一可实现的系统的频率响应函数为H(j)

104(j)426.131(j)3341.42(j)22613.1(j)10000用freqs画出H(j)的幅度响应和相位响应曲线。和f(t)*cos(50t)y(t)和y(t),并根据理论知识解1 2释所得结果。三.仿真分析1 RCH(j)

1/RCj1/RC%(1)RC1=0.001;RC2=0.002;RC3=0.005;a1=1/RC1;a2=1/RC2;a3=1/RC3;A1=[1a1];A2=[1a2];A3=[1a3];B1=a1;B2=a2;B3=a3;w1=linspace(0,3000,200);w2=linspace(0,3000,400);w3=linspace(0,3000,500);H1=freqs(B1,A1,w1);H2=freqs(B2,A2,w2);

xlabel('\omega(title('RC=0.002')subplot(3,1,3);plot(w3,abs(H3));xlabel('\omega(rad/s)')title('RC=0.005')RC=0.00110.500 500 1000 1500 2000 2500 3000(rad/s)RC=0.002H3=freqs(B3,A3,w3);figure(1)subplot(3,1,1);plot(w1,abs(H1));

10.500 500 1000 1500 2000 2500 3000(rad/s)xlabel('\omega(rad/s)')title('RC=0.001')

10.5

RC=0.005subplot(3,1,2);plot(w2,abs(H2));%(2)t=0:0.0001:0.2;

00 500 1000 1500 2000 2500 3000(rad/s)f(t)=cos(100*t)+cos(3000*t)f=cos(100*t)+cos(3000*t);Hw1=a3/(j*100+a3);Hw2=a3/(j*3000+a3);y=abs(Hw1)*cos(100*t)...+abs(Hw2)*cos(3000*t);figure(2)subplot(2,1,1);plot(t,f);xlabel('t (s)')ylabel('f(t)')title('f(t)=cos(100*t)+cos(3000*t)')subplot(2,1,2);plot(t,y);xlabel('t (s)')ylabel('y(t)')

210-1-20 0.02 0.04 0.06 0.08 0.1t(s)10.50-0.5-10 0.02 0.04 0.06 0.08 0.1t(s)

0.12 0.14 0.16 0.18 0.20.12 0.14 0.16 0.18 0.2%(3)clearRC31=0.01;RC32=0.02;RC33=0.05;a31=1/RC31;a32=1/RC32;a33=1/RC33;A31=[1a31];A32=[1a32];A33=[1a33];B31=a31;B32=a32;B33=a33;sys=[tf(B31,A31);tf(B32,A32);tf(B33,A33)];t3=0:0.0005:0.25;f3=abs(10*sin(100*pi*t3));y3=lsim(sys,f3,t3);

figure(1)plot(t3,f3);xlabel('t title('f(t)=|10*sin(100*pi*t3)|')figure(2)plot(t3,y3);legend('RC=0.01时系统的输出',...'RC=0.02时系统的输出','RC=0.05时系统的输出',4)f(t)=|10*sin(100*pi*t3)|10 7968576 45432100 0.05

t(s)

0.15 0.2 0.25

32100 0.05 0.1

RC=0.01时系统的输出RC=0.02时系统的输出RC=0.05时系统的输出0.2 0.25直流分量的计算：fDC

(t)=

12%(1)2%(1)t=0:0.01:2.5;f1t=rectpuls(t-0.5,1);f2t=tripuls(t-1,2,0);ft=f1t+f2t.*cos(50*t);subplot(3,1,1);plot(t,f1t);axis([-0.52.6-0.51.2])ylabel('f1(t)');subplot(3,1,2);plot(t,f2t);axis([-0.52.6-0.51.2])ylabel('f2(t)')subplot(3,1,3);plot(t,ft);ylabel('f(t)')axis([-0.52.6-0.51.2])10.5f0-0.5-0.500.511.522.51f0.50-0.5-0.500.511.522.510.5f0-0.5-0.500.511.522.5

10

0.01sin(100t)dt0

206.37%(2)%(2)a=[126.131341.422613.110000];b=[1000];w=linspace(0,15,500);w1=w*pi;H=freqs(b,a,w1);figure(2)subplot(2,1,1);plot(w,abs(H));xlabel('\omega(\pirad/s)')ylabel('|H(j\omega)|')subplot(2,1,2);plot(w,angle(H));xlabel('\omega(\pirad/s)')ylabel('\phi(\omega)')0.20.15)0.1|0.050051015(rad/s)