计算机仿真实验报告

上传人:无*** 文档编号:135247566 上传时间:2022-08-15 格式:DOC 页数:9 大小:313.50KB
返回 下载 相关 举报
计算机仿真实验报告_第1页
第1页 / 共9页
计算机仿真实验报告_第2页
第2页 / 共9页
计算机仿真实验报告_第3页
第3页 / 共9页
点击查看更多>>
资源描述
山东工商学院计算机仿真及应用实验报告实验七 MATLAB的基本应用(二)及Simulink仿真(验证性实验) 学 院: 专业班级: 实验时间: 学 号: 姓 名: 一、实验目的1、掌握连续信号的仿真和傅里叶分析方法2、掌握连续系统的分析方法(时域分析法,拉氏变换法和傅里叶分析法);3、掌握离散信号的仿真和分析运算方法4、掌握离散系统的分析方法(时域分析法);5、掌握符号运算方法;6、掌握Simulink仿真工具;二、实验原理1、连续信号的仿真和分析法,参考教材第6.1节,重点:单位冲激信号的仿真方法;单位阶跃信号的仿真方法;复指数信号的仿真方法2、连续系统的分析方法,参考教材第6.1节,重点:例6.2,LTI系统的零输入响应的求解方法;例6.3,LTI系统的冲激响应的求解方法例6.5,LTI系统的零状态响应的求解方法例6.6,系统中有重极点时的计算3、系统的频域分析方法,参考教材第6.2节,重点:例6.7,方波分解为多次正弦波之和例6.8:全波整流电压的频谱例6.10:调幅信号通过带通滤波器例6.12:用傅里叶变换计算滤波器的响应和输出4、离散信号的仿真和分析法,参考教材第6.3节,7.1节,重点:单位脉冲序列impseq,单位阶跃序列stepseq例7.1: 序列的相加和相乘例7.2:序列的合成与截取例7.3:序列的移位和周期延拓运算三、实验内容(包括内容,程序,结果)以自我编程练习实验为主,熟悉各种方法和设计,结合课堂讲授,实验练习程序代码。1、 根据教材第6.1节的内容,练习连续信号和系统的时域分析和拉氏变换方法。q602clear,clca=input(输入分母系数向量a=a1,a2,.= );n=length(a)-1;Y0=input(输入初始条件向量Y0=y0,Dy0,D2y0,.= );p=roots(a);V=rot90(vander(p);c=VY0;dt=input(dt= );tf=input(tf= );t=0:dt:tf;y=zeros(1,length(t);for k=1:n y=y+c(k)*exp(p(k)*t); endplot(t,y),gridhold on输入分母系数向量a=a1,a2,.= 3 5 7 1输入初始条件向量Y0=y0,Dy0,D2y0,.= 1 0 0dt= 0.2tf= 8Warning: Imaginary parts of complex X and/or Y arguments ignored. In q602 at 9输入分母系数向量a=a1,a2,.= 3 5 7 1输入初始条件向量Y0=y0,Dy0,D2y0,.= 0 1 0dt= 0.2tf= 8Warning: Imaginary parts of complex X and/or Y arguments ignored. In q602 at 9输入分母系数向量a=a1,a2,.= 3 5 7 1输入初始条件向量Y0=y0,Dy0,D2y0,.= 0 0 1dt= 0.2tf= 8Warning: Imaginary parts of complex X and/or Y arguments ignored. In q602 at 9q603clear,clca=input(多项式分母系数向量a= );b=input(多项式分子系数向量b= );r,p=residue(b,a),disp(解析式h(t)=r(i)*exp(p(i)*t)disp(给出时间数组t=0:dt:tf)dt=input(dt= );tf=input(tf= );t=0:dt:tf;h=zeros(1,length(t);for i=1:length(a)-1 h=h+r(i)*exp(p(i)*t); endplot(t,h),grid多项式分母系数向量a= poly(0 -1+2i -1-2i -2 -5)多项式分子系数向量b= 8 3 1r = 0.6200 0.1300 - 0.3900i 0.1300 + 0.3900i -0.9000 0.0200 p = -5.0000 -1.0000 + 2.0000i -1.0000 - 2.0000i -2.0000 0 解析式h(t)=r(i)*exp(p(i)*t)给出时间数组t=0:dt:tfdt= 0.2tf= 82、 根据教材第6.2节练习傅里叶分析方法。q607t =0:.01:2*pi;y =sin(t);plot(t,y),figure(gcf),pausey =sin(t) + sin(3*t)/3;plot(t,y), pausey =sin(t) + sin(3*t)/3 + sin(5*t)/5 + sin(7*t)/7 + sin(9*t)/9;plot(t,y)z= zeros(10, max(size(t);x = zeros(size(t);for k=1:2:19 x= x + sin(k*t)/k;y(k+1)/2, : )=x;endpause,figure(1),plot(t,y(1:9, : ),gridline(0,pi+0.5,pi/4,pi/4)text(pi+0.5,pi/4,pi/4)halft=ceil(length(t)/2);pause,figure(2),mesh(t(1:halft),1:10,y(:, 1:halft)q609clear, format compactL=0.4;C=10e-6;R=200;Um=100;w1=100*pi;N=input(需分析的谐波次数2N= );n=1:N/1;w=eps,2*n*w1;Us = 4*Um/pi*0.5,-1./(4*n.2-1);z1=j*w*L;z2=1./(j*w*C);z3=R;z23=z2.*z3./(z2+z3);UR=Us.*z23./(z1+z3)UR = Columns 1 through 6 63.6620 - 0.0000i 3.6951 +16.0354i 0.8430 + 0.7970i 0.2077 + 0.1185i 0.0711 + 0.0294i 0.0302 + 0.0099i Columns 7 through 11 0.0149 + 0.0040i 0.0081 + 0.0019i 0.0048 + 0.0010i 0.0030 + 0.0005i 0.0020 + 0.0003i 谐波次数 谐波幅度 谐波相移 0 63.6620 -0.0000 2.0000 16.4556 77.0238 4.0000 1.1601 43.3940 6.0000 0.2391 29.7121 8.0000 0.0769 22.5035 10.0000 0.0318 18.0861 12.0000 0.0154 15.1100 14.0000 0.0083 12.9713 16.0000 0.0049 11.3613 18.0000 0.0031 10.1058 20.0000 0.0020 9.0997disp( 谐波次数 谐波幅度 谐波相移 );disp(2*0,n,abs(UR),angle(UR)*180/pi)3、 根据教材第6.3节的内容,练习离散信号和系统的时域分析方法。q613clear,n0=0;nf=10;ns=3;n1=n0:nf;x1=zeros(1,ns-n0),1,zeros(1,nf-ns);n2=n0:nf;x2=(n1-ns)=0;n3=n0:nf;x3=exp(-0.2+0.5j)*n3);subplot(2,2,1),stem(n1,x1);title(单位脉冲序列)subplot(2,2,2),stem(n2,x2);title(单位阶跃序列)subplot(2,2,3),stem(n3,real(x3);line(0,10,0,0)title(复指数序列),ylabel(实部)subplot(2,2,4),stem(n3,imag(x3);line(0,10,0,0)ylabel(虚部)4、 根据教材第9.1节的内容,以及表9.1,练习符号运算。5、 根据教材第9.2节练习Simulink方法。四、实验总结与体会用matlab解决各种信号、电路问题很方便,省去了大量的计算,提高了准确性。但是编程需要对matlab有很高的运用能力,必须多加练习。9 / 9文档可自由编辑打印
展开阅读全文
相关资源
正为您匹配相似的精品文档
相关搜索

最新文档


当前位置:首页 > 压缩资料 > 基础医学


copyright@ 2023-2025  zhuangpeitu.com 装配图网版权所有   联系电话:18123376007

备案号:ICP2024067431-1 川公网安备51140202000466号


本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!