资源描述：
《信号与线性系统实验报告材料.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、课程实验报告课程名称：信号与线性系统专业班级：计算机科学与技术11级10班学号：U201114445姓名：王涛指导教师：伟平报告日期：2013/7/1计算机科学与技术学院实验一连续时间系统的时域分析%1_1_1.mT=0.01;tx=0:T:2;x=rectpuls((tx-1),2).*exp(-tx);th=0:T:2;h=rectpuls((th-1),2).*2;t=(0+0):T:(2+2);y=conv(x,h);figuresubplot(3,1,1);plot(tx,x)ylabel('输入激励');subplot(3,1,2);plot(t
2、h,h)ylabel('单位冲激响应');subplot(3,1,3);plot(t,y)ylabel('输出响应');%1_1_2.mT=0.01;tx=-4:T:4;x=rectpuls((tx-0),8).*(1-abs(tx)/4);th=0:T:4;h=(th>=0);t=(-4):T:8;y=conv(x,h);figuresubplot(3,1,1);plot(tx,x)ylabel('输入激励');subplot(3,1,2);plot(th,h)ylabel('单位冲激响应');subplot(3,1,3);plot(t,y)ylabel(
3、'输出响应');%1_2_1.mR=10;%电阻值C=4;%电容值U0=2;%电容电压初始值S=[num2str(C),'*',num2str(R),'*Dy+','y=heaviside(t)'];%得到微分方程的字符串表达式init=['y(0)=',num2str(U0)];%得到初始条件的字符串表达式y=dsolve(S,init,'t');%求解微分方程，得到符号解t=2:0.01:2*pi;x=ones(1,length(t));figuresubplot(1,2,1);plot(t,x)title('激励电压');subplot(1,2,2);
4、%2号子图显示输出响应ysezplot(y,[0,2*pi])title('输出响应')%1_2_2.mR=0.5;%电阻值L=1;%电容值I0=2;%电感电流初始值S=[num2str(L),'*Dy+',num2str(R),'*y=heaviside(t)'];%得到微分方程的字符串表达式init=['y(0)=',num2str(I0)];%得到初始条件的字符串表达式y=dsolve(S,init,'t');%求解微分方程，得到符号解t=0:0.01:2*pi;x=ones(1,length(t));%得到时间围t的阶跃信号u(t)的离散抽样序列xf
5、iguresubplot(1,2,1);%多子图显示，将图形框分为1x2个子图，1号子图显示阶跃信号xplot(t,x)title('激励电压');subplot(1,2,2);%2号子图显示输出响应ysezplot(y,[0,2*pi])title('输出响应')实验分析：本次实验应用连续系统时域分析的两种方法，卷积计算和微分方程的求解。实验二信号的傅里叶分析实验%2_1_1.mN=10;%希望看到的谐波次数symstTk;%定义时间t、周期Ta、下标k三个符号y=subs(sym('Heaviside(t+T/20)-Heaviside(t-T/20)'
6、),'T',TT);%锯齿脉冲基础波形，tao/T=1/20%修改上面的sym中的参数，将T/20改为T/40皆可%得到占空比为1/20的。A0=int(y,t,-TT/2,TT/2)/TT;%直流分量A0%利用符号法得到各频率分量的复数Fourier系数Ak的符号表达式Ak=int(y*exp(-2*i*pi*k*t/TT),t,-TT/2,TT/2)/TT;fk=sym(Ak)*sym(exp(2*i*k*pi*t/TT));%计算k在区间[-N,N]的Fourier系数，存入向量a（长度为2*N+1）form=-N:-1%a(m+N+1)=numeri
7、c(subs(Ak,k,m));a(m+N+1)=double(subs(Ak,k,m));enda(N+1)=double(sym(A0));form=1:Na(m+N+1)=double(subs(Ak,k,m));end%for%利用symsum命令进行符号求和，得到综合信号ff=symsum(fk,k,-N,-1)+A0+symsum(fk,k,1,N);%图形显示结果figuren=-N:N;as=abs(a)*2;%由Fourier系数得到频谱幅度subplot(3,1,1)ezplot(y,[-TT,TT])ylabel('原函数');subp
8、lot(3,1,2)ezplot(f,[-TT,TT
