MATLAB simulink 通信系统建模与仿真第五章部分习题答案

MATLAB原理及应用实验报告第五章答案

《MATLAB 原理及应用》实验报告三．练习题1、编写程序绘制如下图形（时间间隔自己设定）图中曲线标注:sin(t)标注在(1.7*pi,-0.3)，cos(t)标注在(1.6*pi,0.8) 2、使用subplot指令在不同区域画出下列曲线y=sinAx;y=cosAx;y=Ae x ,y=x^2。

x 的取值范围是]15,0[，A 为自己的学号。

3、画出衰减振荡曲线t e y t 3sin 3-=及其它的包络线30t e y -=。

t 的取值范围是]4,0[π，图的标题字体大小(fontsize)为10号，字体(fontname)宋体，字体风格（fontangle ）为斜体字体，并选择合适的位置在图上标注自己的学号，姓名，以及作图日期，字体大小、字体风格自己指定。

四．实验程序及结程序为1：t=0:0.1:4*pi;y=sin(t);y1=cos(t);plot(t,y,'r',t,y1,'--b') text(1.7*pi,-0.3,'sint') text(1.6*pi,0.8,'cost')2468101214-1-0.8-0.6-0.4-0.200.20.40.60.812：x=0:0.1:15; y1=sin(1*x); subplot(2,2,1); plot(x,y1); y2=cos(1*x); subplot(2,2,2); plot(x,y2);y3=(1*exp(x)); subplot(2,2,3); plot(x,y3); y4=x.^2; subplot(2,2,4); plot(x,y4);051015-1-0.500.51051015-1-0.500.510510151234x 1060510151002003003 程序为：t=0:0.01:4*pi;y0=exp(-t/3);y1=-exp(-t/3);y=exp(-t/3).*sin(3*t);plot(t,y0,'--b',t,y,'r',t,y1,'--b');title('y=e^(-t/3)sin(3t)及它的包络线','fontsize',10','fontweigh','bold','fontangle','italic');text(9,-0.4,['学号:'],'fontsize',10','fontweigh','bold','fontangle','italic');text(9.5,-0.6,['姓名:某某某'],'fontsize',10','fontweigh','bold','fontangle','italic');text(10,-0.8,['2013.12.01']);。

MATLAB语言与控制系统仿真参考答案第5章

5、6控制系统的时域响应MATLAB 仿真实训5、6、1实训目的 123 4 5学会利用MATLAB绘制系统的单位阶跃响应曲线，掌握读取系统动态性能指标的方法学会利用MATLAB 绘制系统的单位脉冲响应曲线的方法掌握利用MATLAB 绘制系统的零输入响应曲线的方法掌握利用MATLAB 绘制系统的一般输入响应曲线的方法学会通过仿真曲线读取相关信息，并依据有关信息进行系统的时域分析。

5、6、2实训内容 1、编写程序求取下列各系统的单位阶跃响应，完成表5-5并记录相关曲线。

G(s) 16 ~2 s 3.2s 16 G 2(s) 16~2 s 2.4s 16 G a (s)16 s 21.6s 16G 4(s)16 -2s s 16解:>> n1=16; >> d 仁[1,3、2,16]; >> sys1=tf( n1,d1); >> step(sys1)-X r倚 25 3 7百'imeQ -4o.a■rp3±±-i_l r v>> n2=16; >> d2=[1,2、4,16]; >> sys2=tf( n2,d2); >> step(sys2)041.4D Q粘EEdE'luaJ 05 1 15 2 2.3 3 55 4<1.5Tim (lie)11-o02:>> n3=16;>> d3=[1,1、6,16];>> sys3=tf( n3,d3);>> step(sys3)lim* (t«c) >> n4=16; >> d4=[1,1,16];>> sys4=tf( n4,d4);>> step(sys4)表5-5序号nc maxt pt s (5%)计算值实验计算值实验计算值实验值1 0、4 4 1、2538 1、250、85690、8632、 18752、12 0、34 1 3723 1、370、 82330、8282、 91672、813 0、24 1、52661、530、80160、8 4 3750 4、940、125 4 1、6731 1、670、 79160、8037、 00007、33>> z1=0、4; w=4;cmax 1=1+exp(-z1*pi/sqrt(1-z1A2)); tp1=pi/(w*sqrt(1-z1A 2)); ts1=3、5/(z1*w); [cmax1,tp1,ts1] ans =>> z2=0、3;w=4;cmax2=1+exp(-z2*pi/sqrt(1-z2A2)); tp2=pi/(w*sqrt(1-z2A2)); ts2=3、5/(z2*w); [cmax2,tp2,ts2]ELIIds1、 2538 0、 85692、 1875ans =1、3723 0、82332、 9167>> z3=0、2;w=4;cmax3=1+exp(-z3*pi/sqrt(1-z3A2));tp3=pi/(w*sqrt(1-z3A2));ts3=3、5/(z3*w);[cmax3,tp3,ts3] ans =1、5266 0、8016 4、3750 >> z4=0、125;w=4;cmax4=1+exp(-z4*pi/sqrt(1-z4A2));tp4=pi/(w*sqrt(1-z4A2));ts4=3、5/(z4*w);[cmax4,tp4,ts4]ans =1、6731 0、7916 7、0000说明：对于二阶欠阻尼系统（0 1）,若系统的闭环传递函数为则系统单位阶跃响应的输出最大值cmaxt pn 1 2调整时间估算值t s3.5n（以5%为误差带）t s 4.4（以2%为误差带）2、已知二阶系统的闭环传递函数如下，编程求取系统的单位阶跃响应并完成表峰值时间曲线。

MATLAB语言基础与应用(第二版)第5章习题答案

第5章习题与答案5.1用矩阵三角分解方法解方程组123123123214453186920x x x x x x x x x +-=⎧⎪-+=⎨⎪+-=⎩ 解答：>>A=[2 1 -1;4 -1 3;6 9 -1] A =2 1 -1 4 -13 6 9 -1 >>b=[14 18 20]; b =14 18 20 >> [L, U, P]=lu(A) L =1.0000 0 0 0.6667 1.0000 0 0.3333 0.2857 1.0000 U =6.0000 9.0000 -1.0000 0 -7.0000 3.6667 0 0 -1.7143 P =0 0 1 0 1 0 1 0 0 >> y=backsub(L,P*b’) y =20.0000 4.6667 6.0000 >> x=backsub(U,y) x =6.5000 -2.5000 -3.5000 5.2 Cholesky 分解方法解方程组123121332352233127x x x x x x x ++=⎧⎪+=⎨⎪+=⎩ 解答：>> A=[3 2 3;2 2 0;3 0 12] A =3 2 32 2 03 0 12>> b=[5;3;7]b =537>> L=chol(A)L =1.7321 1.1547 1.73210 0.8165 -2.44950 0 1.7321>> y=backsub(L,b)y =-11.6871 15.7986 4.0415>> x=backsub(L',y)x =-6.7475 28.8917 49.93995.3解答:观察数据点图形>> x=0:0.5:2.5x =0 0.5000 1.0000 1.5000 2.0000 2.5000 >> y=[2.0 1.1 0.9 0.6 0.4 0.3]y =2.0000 1.1000 0.9000 0.6000 0.4000 0.3000 >> plot(x,y)图5.1 离散点分布示意图从图5.1观察数据点分布，用二次曲线拟合。

第五章 Simulink系统建模与仿真

第五章 Simulink建模与仿真
本章重点
Simulink基本结构 Simulink模块系统模型及仿真
一、Simulink简介
Simulink 是MATLAB 的工具箱之一，提供交互式动态系统
建模、仿真和分析的图形环境
可以针对控制系统、信号处理及通信系统等进行系统的建模、仿真、分析等工作可以处理的系统包括：线性、非线性系统；离散、连续及混合系统；单任务、多任务离散事件系统。
从模块库中选择合适的功能子模块并移至编辑窗口中，按设计要求设置好各模块的参数，再将这些模块连接成系统 Simulink的仿真过程就是给系统加入合适的输入信号模块和输出检测模块，运行系统，修改参数及观察输出结果等
过程
二、Simulink的基本结构
Simulink窗口的打开
命令窗口：simulink 工具栏图标：
三、Simulink模型创建
7、信号线的标志
信号线注释：双击需要添加注释的信号线，在弹出的文本编辑框中输入信号线的注释内容
信号线上附加说明：（1）粗线表示向量信号：选中菜单Forma t|Wide nonscalar lines 即可以把图中传递向量信号的信号线用粗线标出；（2）显示数据类型及信号维数：选择菜单Format|Port data types 及Format|Signaldimensions，即可在信号线上显示前一个输出的数据类型及输入/输出信号的维数；（3）信号线彩色显示：选择菜单Format|Sample Time Color，SIMULINK 将用不同颜色显示采样频率不同的模块和信号线，默认红色表示最高采样频率，黑色表示连续信号流经的模块及线。
同一窗口内的模块复制：（1）按住鼠标右键，拖动鼠标到目标

MATLAB语言与控制系统仿真-参考答案-第5章

控制系统的时域响应MATLAB 仿真实训实训目的 1. 学会利用MATLAB 绘制系统的单位阶跃响应曲线，掌握读取系统动态性能指标的方法； 2. 学会利用MATLAB 绘制系统的单位脉冲响应曲线的方法； 3. 掌握利用MATLAB 绘制系统的零输入响应曲线的方法； 4. 掌握利用MATLAB 绘制系统的一般输入响应曲线的方法；5.学会通过仿真曲线读取相关信息，并依据有关信息进行系统的时域分析。

实训内容1.编写程序求取下列各系统的单位阶跃响应，完成表5-5并记录相关曲线。

162.316)(21++=s s s G 164.216)(22++=s s s G ;166.116)(23++=s s s G 1616)(24++=s s s G 解：>> n1=16;>> d1=[1,,16]; >> sys1=tf(n1,d1); >> step(sys1)>> n2=16; >> d2=[1,,16]; >> sys2=tf(n2,d2); ~>> step(sys2)>> n3=16;>> d3=[1,,16]; >> sys3=tf(n3,d3); >> step(sys3)>> n4=16;>> d4=[1,1,16]; *>> sys4=tf(n4,d4); >> step(sys4)序号ξnωm axcptst(%5=∆)计算值实验计算值实验计算值实验值1】42[434~44】w=4;cmax1=1+exp(-z1*pi/sqrt(1-z1^2));tp1=pi/(w*sqrt(1-z1^2));ts1=(z1*w);`[cmax1,tp1,ts1]ans =>> z2=;w=4;cmax2=1+exp(-z2*pi/sqrt(1-z2^2));tp2=pi/(w*sqrt(1-z2^2));ts2=(z2*w);[cmax2,tp2,ts2]）ans =>> z3=; w=4;cmax3=1+exp(-z3*pi/sqrt(1-z3^2)); tp3=pi/(w*sqrt(1-z3^2)); ts3=(z3*w); [cmax3,tp3,ts3]》ans =>> z4=; w=4;cmax4=1+exp(-z4*pi/sqrt(1-z4^2)); tp4=pi/(w*sqrt(1-z4^2)); ts4=(z4*w); [cmax4,tp4,ts4] ^ ans =说明：对于二阶欠阻尼系统（10<<ξ），若系统的闭环传递函数为2222)(nn ns s s Φωξωω++= 则系统单位阶跃响应的输出最大值21max 1ξξπ--+=ec峰值时间21ξωπ-=n p t调整时间估算值ns t ξω5.3= （以5%为误差带） ns t ξω4.4=（以2%为误差带）￥2.已知二阶系统的闭环传递函数如下，编程求取系统的单位阶跃响应并完成表5-6，记录相关曲线。

MATlab 仿真部分习题答案

[4.1]控制系统结构如图4.1所示（1）利用MATLAB对以上单位负反馈控制系统建立传递函数；（2）将第一问中求得的传递函数模型转化为零极点增益形式和状态空间形式。

解:（1）num=[2 2];den=[1 2 1];[num1,den1]=cloop(num,den);sys=tf(num1,den1)程序运行结果如下：Transfer function:2 s + 2-------------s^2 + 4 s + 3（2）[z,p,k]=tf2zp(num1,den1);g_zpk=zpk(z,p,k);[A,B,C,D]=zp2ss(z,p,k);g_ss=ss(A,B,C,D)程序运行结果如下：z = -1 p = -3 -1 k = 2a = x1 x2x1 -4 -1.732x2 1.732 0b = u1x1 1x2 0c = x1 x2y1 2 1.155d = u1y1 0K[5.1]设单位负反馈控制系统的开环传递函数为G(S)=+s(s^2+7s17)（1）试绘制k=10、100时闭环系统的阶跃响应曲线，并计算稳态误差、上升时间、超调量和过渡过程时间；（2）绘制k=1000时闭环系统的阶跃响应曲线，与k=10、100时所得的结果相比较，分析增益系数与系统稳定性的关系；解：（1）k=10时，K=100时，K=10时，利用MA TLAB工作区输入程序：num=[10];den=[1,7,17,0];[z,p,k]=tf2zp( num,den);运行得z,p,k的值p= -3.5+2.1794*i -3.5-2.1794*i k=10G=zpk([ ],[-3.5+2.1794*i,-3.5-2.1794*i],10)；c=dcgain(G);[y,t]=step(G);plot(t,y)[Y,K]=max(y);timetopeak=t(k);percentovershoot=100*(y-c)/cn=1;while y(n)<c n=n+1;end risetime=t(n) i=length(t);while(y(i)>0.98*c&y(i)<1.02*c) i=i-1;End setllingtime=t(i) 运行程序结果为：稳定值c=0.5882，响应时间setllingtime=1.0096s,上升时间risetime=1.1989s,最大峰值时间timetopeak=1.4356,超调量percentovershoot=0.555% 同理得k=100时，稳定值c=5.882，响应时间setllingtime=1.0096s, 上升时间risetime=1.1989s,最大峰值时间timetopeak=1.4356,超调量percentovershoot=0.555%（3）k=1000时，由其响应曲线可知，增益系数越大，其稳定性越差。

MATLAB教程课后实验报告题目及解答[第一至第五章]

MATLAB教程实验报告实验项目名称实验一 Matlab基本操作学生姓名汪德旺专业班级 09数教（1）班学号 0301090131实验成绩日期一. 实验目的和要求1、了解MATLAB 的开发环境。

2、熟悉Matlab的基本操作。

3、掌握建立矩阵的方法。

4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。

5、填写实验报告，实验报告文件取名为report1.doc。

6、于邮件附件形式将实验报告文件report1.doc 发到邮箱*******************，邮件主题为班级学号姓名，如：09数教1班15号张三。

二、实验内容1、先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。

（1）0 122sin851ze =+（2）2212 1ln(0.4552i z x x+⎡⎤==⎢⎥-⎣⎦其中（3）0.30.330.3sin(0.3)ln,22a ae e az a--+=++a=-3.0,-2.9,-2.8,…, 2.8, 2.9,3.0（4）2242,011,12,0:0.5:2.521,23t tz t tt t t⎧≤<⎪=-≤<=⎨⎪-+≤<⎩其中t2.已知：1234413134787,2033657327 A B--⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：（1）A+6*B和A-B+I(其中I为单位矩阵) （2）A*B和A.*B（3）A^3和A.^3（4）A/B和B\A（5）[A,B]和[A([1,3],:);B^2]3、设有矩阵A 和B123453016678910A=,B=17-691112131415023-41617181920970212223242541311⎡⎤⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦（1）求它们的乘积C 。

（2）将矩阵C 的右下角3*2子矩阵赋给D 。

（3）查看MATLAB 工作空间的使用情况。

第五讲基于MATLAB-Simulink的建模与仿真

MATLAB的功能包括：数值分析，数值和符号计算，工程和科学绘图，通讯和控制系统的设计与仿真，数字图像与信号处理，财务与金融工程等。
MATLAB软件简介？
MATLAB软件的典型应用领域：
❖科学研究； ❖工程技术应用研究 ❖CAI（Computer Aided Instruct） ❖数学实验（Mathematical Experiment） ❖数学建模（Mathematical Modeling）
模型 Transfer-Fcn：线性传递函数模型 Zero-Pole：以零极点表示的传递
函数模型 Memory：存储上一时刻的状态值 Transport Delay：输入信号延时一个给定时间再输出 Variable Transport Delay：输入信号延时一个可变时间再输出
✓ 离散模块（Discrete）
For循环不能用For循环内重新赋值循环变
量n来终止。
在For循环中循环控制量的范围可以是任
何有效的MATLAB矩阵。比如
data=[11 9 45 6; 7 16 -1 5];
for n=data
x=n(1)-n(2)
end 这时程序的输出有四个数值，分别是矩阵
data的两列相减的结果
x = 4 x = -7
x = 46 x = 1
For循环可按需要嵌套，即For循环体内的命令组中可以出现另一个For循环体，这体现了 For循环体也是命令组。比如 for n=1:5
for m=5:-1:1
A(n,m)=n^2+m^2; End
end
MATLAB软件简介？
While-end循环以不定的次数求一组语句的值。 Whil-end 循环的一般形式是： while expression（控制表达式） {commands} end 只要在控制表达式（expression）里的所有元素为真，就执行While和end语句之间的命令串（{commands}）。

5章习题答案matlab

张卫华 MATLAB课堂
例13：
分析以下程序，并运行观察。
clf;x=3*pi*(-1:0.05:1);y=x; [X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps; Z=sin(R)./R; h=surf(X,Y,Z);colormap(jet); axis off n=12;mmm=moviein(n); for i=1:n rotate(h,[0 0 1],25); mmm(:,i)=getframe; end movie(mmm,5,10)
张卫华 MATLAB课堂
例5程序：
fplot('cos(tan(pi*x))',[ 0,1],1e-4)
张卫华 MATLAB课堂
例6：
绘制r=sin(t)cos(t)的极坐标图，并标记数据点。
张卫华 MATLAB课堂例6 Nhomakorabea序：t=0:pi/50:2*pi; r=sin(t).*cos(t); polar(t,r,'-*');
title('y1=0.2e^{-0.5x}cos(4\pix) 和y2=2e^{0.5x}cos(4\pix)比较') text(2,2*exp(-1),'\fontsize{20}\bf \leftarrow x_{1}=2,y_{1}=2e^{-1}')
张卫华 MATLAB课堂
例5：
绘制f(x)=cos(tan(πx))的曲线
张卫华 MATLAB课堂
例11：
绘制柱形，剪切掉x、y小于零的部分
张卫华 MATLAB课堂
例11程序：
t=linspace(0,2*pi,100); r=1-exp(-t/2).*cos(4*t); [X,Y,Z]=cylinder(r,60); ii=find(X<0&Y<0); Z(ii)=NaN; surf(X,Y,Z);colormap(spring), shading interp light('position',[-3,-1,3],'style','local')

MATLAB基础习题第五章习题答案

（1）如果当前时间为 7-8 点则提醒用户，该吃早饭了；（2）如果当前时间为 9-11 点则提醒用户，该学习了；（3）如果当前时间为 12-14 点则提醒用户，该午休了；（4）如果当前时间为 14-17 点则提醒用户，该锻炼了；答： %%程序为：remind.fig、remind.m

6.求解六元线性方程组：（1）方程组的系数矩阵由用户通过键盘输入；（2）得到系数矩阵后给出方程的解；（3）程序要具有友好性答：
%% %该模块实现的功能是：求解六元线性方程组，方程组的系数矩阵由用户通过键盘输入；得到系数矩阵后给出方程的解； clc; clear all; close all; %% %方程输入模块 x=inputdlg({'第一个方程系数','第二个','第三个','第四个','第五个','第六个'}); A=cell2mat(x); %% y=inputdlg({'第一个方程等号右边数','第二个','第三个','第四个','第五个','第六个'});
you=input('请做出你的选择，石头（1），剪刀（2），布（3）： '); end disp('您的选择是： '); disp(b(2*you-1:2*you)); compute=ceil(3*rand(1,1)); disp('电脑的选择是：'); disp(b(2*compute-1:2*compute)); %% %%输赢判断模块 %如果电脑与选手出的一样，则显示平手，否则根据石头剪刀布的规则来判断输赢 if you==compute disp('平手'); end switch(you-compute) case{1,-2} disp('您输了'); case{-1,2} disp('您赢了'); end %% %是否继续判别模块 n=input('是否继续玩该游戏？否（0），是（1）'); while n~=0&n~=1 disp('您输入的不是正确数字，请正确输入'); n=input('是否继续玩该游戏？否（0），是（1）'); end end >> 请选择，石头（1），剪刀（2），布（3）： 1 您的选择是：石头电脑的选择是：、剪您赢了是否继续玩该游戏？否（0），是（1）1 请选择，石头（1），剪刀（2），布（3）： 2 您的选择是：、剪电脑的选择是：石头您输了是否继续玩该游戏？否（0），是（1）0 >> 5.编写一个日程提醒程序实现如下功能：

MATLAB课后习题集附标准答案

第2章 MATLAB概论1、与其他计算机语言相比较，MATLAB语言突出的特点是什么？答：起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富，可扩展性强.2、MATLAB系统由那些部分组成？答：开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装.第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可.矚慫润厲钐瘗睞枥庑赖。

4、MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？聞創沟燴鐺險爱氇谴净。

答：在MATLAB操作桌面上有五个窗口，在每个窗口的右下角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口称为独立的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面称为独立窗口，在独立窗口的view菜单中选择Dock，菜单项就可以将独立的窗口重新防止的桌面上.残骛楼諍锩瀨濟溆塹籟。

5、如何启动M文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动.在命令窗口中键入edit命令时也可以启动M文件编辑/调试器.酽锕极額閉镇桧猪訣锥。

6、存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可.彈贸摄尔霁毙攬砖卤庑。

7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中.謀荞抟箧飆鐸怼类蒋薔。

Matlab第五章答案

第一题（1）a=[1 9 8;7 2 5;3 -2 7] %产生矩阵det(a) %检验是否可逆ans=-442,非0，可逆div(a) %求逆矩阵（2）b=[1 0 -7 5;0 -26 7 2;7 4 3 5;8 -3 2 15]det(b)div(b)第二题(1)A=[1 2 3;2 2 5;3 5 1];B=[11 12 31];X=B/A(2) A=[3 1 0 5;0 6 7 3;0 4 3 0;2 -1 2 6];B=[2 4 7 8];X=B/A第三题（1）t=[1 2 3 4 5 6 7 8 9 10]';y=[4.842 4.362 3.754 3.368 3.169 3.083 3.034 3.016 3.012 3.005]';A=[ones(size(t)) exp(-t)];C=A\yT=[0:.1:10]';Y=[ones(size(T)) exp(-T)]*C;plot(T,Y,'-',t,y,'o')title( '采用y(t)≈c1+c2e–t的拟合' )xlabel('\itt'), ylabel('\ity')（2）t=[0 .3 .8 1.1 1.6 2.3]';y=[.82 .72 .63 .60 .55 .50]';A=[ones(size(t)) t.*exp(-t)];C=A\y;T=[0:.1:2.5]';Y=[ones(size(T)) T.*exp(-T)]*C;plot(T,Y,'-',t,y,'o')title('采用y(t)≈d1+d2te–t拟合')xlabel('\itt'), ylabel('\ity')第四题A=[11.59 12.81 15.66; 15.2 4.18 13.61; 10.597.59 9.22];[L,U]=lu(A)[Q R]=qr(A)B=[16.00 4.41 -10.37 -21.61; 0.88 -20.04 12.86 8.56; -1.43 10.71 18.81 -5.99; -12.48 24.35-23.9 10.34];[C,D]=lu(B)[E F]=qr(B)第五题（1）A=[5 -5 -6;3 -2 5;2 -1 -4];x0=[1;-4;5];X=[];for t=0:.01:1X=[X expm(t*A)*x0];endplot3(X(1,:),X(2,:),X(3,:),'-o')grid on(2)A=[1 2 -3 1;3 0 1 -2;1 -2 0 5;2 3 0 1];x0=[1;-1;2;1];X=[];for t=0:.01:1X=[X expm(t*A)*x0];endplot3(X(1,:),X(2,:),X(3,:),'-o')grid on第六题（1）A=[11.59 12.81 15.66; 15.2 4.18 13.61;10.59 7.59 9.22];lambda=eig(A)[V,D]=eig(A)（2）B=[16.00 4.41 -10.37 -21.61; 0.88 -20.04 12.86 8.56; -1.43 10.71 18.81 -5.99; -12.48 24.35 -23.9 10.34];lambda=eig(B)[V,D]=eig(B)第七题（1）x=[1 2 3 4 5 6 7 8 9 10];y=[15.0 39.5 66.0 85.5 89.0 67.5 12.0 -86.4 -236.9 -448.4];p=polyfit(x,y,2);x2=1:.1:10;y2=polyval(p,x2);figure(1)plot(x,y,'o',x2,y2)grid ontitle('二阶多项式曲线拟合')(2)x=[1 2 3 4 5 6 7 8 9 10];y=[15.0 39.5 66.0 85.5 89.0 67.5 12.0 -86.4 -236.9 -448.4];p=polyfit(x,y,3);x2=1:.1:10;y2=polyval(p,x2);figure(1)plot(x,y,'o',x2,y2)grid ontitle('三阶多项式曲线拟合')第八题p1=[1,-2-3,4,2];p2=[1,-7,5,31,-30];p3=[1,-1,-25,25];p4=[-2,3,1,5,8,0];[L1,U1]=lu(p1)r1=roots(p1)[L2,U2]=lu(p2)r2=roots(p2)[L3,U3]=lu(p3)r3=roots(p3)[L4,U4]=lu(p4)r4=roots(p4)第九题p1=[1,-2-3,4,2];p2=[1,-7,5,31,-30];p3=[1,-1,-25,25];p4=[-2,3,1,5,8];p1_x=polyval(p1,[-1.5,2.1,3.5]) p2_x=polyval(p2,[-1.5,2.1,3.5]) p3_x=polyval(p3,[-1.5,2.1,3.5]) p4_x=polyval(p4,[-1.5,2.1,3.5])第十题a=[2,3,-4];b=[4,-2,5];c=[3,0,-2,5,6];d1=conv(a,b)[d2,r2]=deconv(c,a)[d3,r3]=deconv(c,b)第十一题a=[2,3,-4];b=[4,-2,5];c=[3,0,-2,5,6];dao1=polyder(a,b)[dao2,r2]=polyder(c,a)[dao3,r3]=polyder(c,b)第十二题x=-5:.25:5;y=10*exp(-x);xi=-5:5;y1=interp1(x,y,xi,'nearest');y2=interp1(x,y,xi,'linear');y3=interp1(x,y,xi,'spline');y4=interp1(x,y,xi,'cubic'); figure(1);subplot(2,2,1)plot(x,y,'-',xi,y1,'o');title('最邻近内插');grid on;xlabel('x');ylabel('y');subplot(2,2,2)plot(x,y,'-',xi,y2,'o');title('线性内插');grid on;xlabel('x');ylabel('y');subplot(2,2,3)plot(x,y,'-',xi,y3,'o');title('三次样条内插');grid on;xlabel('x');ylabel('y');subplot(2,2,4)plot(x,y,'-',xi,y4,'o');title('三次曲线内插');grid on;xlabel('x');ylabel('y');第十三题x=rand(1,50);y=randn(1,50);minx=min(x)miny=min(y)maxx=max(x)maxy=max(y)avx=mean(x)avy=mean(y)Ex=(std(x)).^2Ey=(std(y)).^2第十四题t=[0 .2 .4 .6 .8 1.0 2.0 5.0 ]';y=[1.0 1.51 1.88 2.13 2.29 2.40 2.60 24.00]'; X1=[ones(size(t)) t t.^2];a=X1\y;X2=[ones(size(t)) exp(-t) t.*exp(-t)];b=X2\y;T=[0:.1:6]';Y1=[ones(size(T)) T T.^2]*a;Y2=[ones(size(T)) exp(-T) T.*exp(-T)]*b; figure(1)subplot(1,2,1)plot(T,Y1,'-',t,y,'o'),grid ontitle('多项式回归')subplot(1,2,2)plot(T,Y2,'-',t,y,'o'),grid ontitle('指数函数回归')第十五题t=0:1/119:1;x=3*sin(2*pi*20*t)+10*sin(2*pi*200*t+pi/4)+10*randn(size(t)); y=fft(x);m=abs(y);f=(0:length(y) -1)'*119/length(y);figure(1)subplot(2,1,1),plot(t,x),grid ontitle('被噪声污染的信号')ylabel('Input \itx'),xlabel('Time ')subplot(2,1,2),plot(f,m)ylabel('Abs. Magnitude'),grid onxlabel('Frequency (Hertz)')第十六题w=input('w=');t=0:1/119:1;x1=sin(w.*t)+randn(size(t));x2=cos(w.*t)+randn(size(t));x3=sin(w.*t)+randn(size(t));a=corrcoef(x1,x2)b=corrcoef(x1,x3)若没有正弦分量w=input('w=');t=0:1/119:1;x1=randn(size(t));x2=randn(size(t));x3=randn(size(t));a=corrcoef(x1,x2)b=corrcoef(x1,x3)第十七题z1=quad('exp(-2*t)',0,2)z2=quad('exp(2*t)',0,2)z3=quad('exp(t.^2-3*t+.5)',-1,1)第十八题function y=five(x)y=exp(-x)-1.5*exp(2*cos(2*x));%主函数x0=input('x0='); %执行时，按要求输入[-1,1]z=fzero('five',x0)第十九题function f=five(x,y)f=exp(-x.*y)-2*x.*y;%主函数z=dblquad('five',0,1,-1,1)第二十题function dy=five(t,y)dy=[0.5-y(1);y(1)-4*y(2)];%主函数X0=[1; -0.5];tspan=[0,25];[T,X]=ode45('five',tspan,X0);figure(1)subplot(2,1,1),plot(T,X(:,1),'r'),title('x_{1}'),grid onsubplot(2,1,2),plot(T,X(:,2),'k'),title('x_{2}'),grid onfigure(2)plot(X(:,1),X(:,2)),title('系统轨迹'),grid onxlabel('x_{1}'),ylabel('x_{2}')。

MATLAb与数学实验第五章习题解答

01 10 00 0 -2
(1,1)
1
(3,2)
1
(2,3)
2
(1,4)
-1
(3,5)
3
A2 =
1 0 0 -1 0 00200
01003 (3) A1 =
(1,1)
1
(5,1)
2
(4,2)
3
(3,3)
1
(2,4)
3
(1,5)
2
(5,5)
1
A2 =
10002 00030 00100 03000 20001 18．创建一个 4 阶稀疏矩阵，使副对角线上元素为 1 答 A=sparse(1:4,1:4,1)
B=
-2 0 0 010 004
C=
-2 1 4
14 47 7 10
D=
-2 0 0 -2 1 4 010147 0 0 4 4 7 10
F=
-2 1 4 147 4 7 10 -2 0 0 010 004
(4) C=[-2 1 4;1 4 7;4 7 10] C(:,1)=[]
C=
-2 1 4 147 4 7 10
A=
(1,1)
1
(2,2)
1
(3,3)
1
(4,4)
1
19．创建如下稀疏矩阵，查看其信息，并将其还原成全元素矩阵
1 0 2 0 0 1 0 1 0 1 0 0
0 1 0 2 0
0 2 0 2 0
2
0

（1） 3 0 1 0 2 (2) 0 0 3 0 3 0 3
工种
天数在木工家的工作天数在电工家的工作天数在油漆工家的工作天数
木工
2 4 4

Matlab与通信仿真实验报告册答案

Matlab 与通信仿真34.0000 22.0000 62.00001. 在Command Window 里面计算_7 8 914 5 6，计算：C = AxB,D =A + B,A C,C/B ； 12 3>> a=[1 2 3;4 5 6;7 8 9]; >> b=[7 8 9;4 5 6;1 2 3]; >> c=a*b c =18 24 30 5469 84 90114138>> d=a+b d =8 10 12 8 10 12 81012>> e=a\cWarning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 1.54 佃76e-018.(1) (3 5 8)- 5 10 ;>>(3+5+8)/5*10 ans=32>>sin(3*pi)/sqrt(9/5) ans=2.7384e-161 2 3(3)A= 4 5 6 ,B 78 9JMatlab 与通信仿真-50.0000-23.0000 -100.0000 28.0000 16.000056.0000>> f=c/bWarning: Matrix is singular to working precision.NaN NaN NaNNaN -Inf Inf NaNNaNNaN3 1.2 4(4) A= 7.5 6.6 3.1，求A ：A, A ；(求矩阵的行列式)5.4 3.46.1 _>> a=[3 1.2 4;7.5 6.6 3.1;5.4 3.4 6.1]; >> a' ans =2.1555 0.4555 -1.6449 -2.1040-0.2393 1.5013 -0.7354-0.26980.7833>> det(a) ans =13.78803.0000 1.20004.0000>> inv (c )7.5000 6.6000 3.10005.4000 3.40006.1000Matlab 与通信仿真1 2 0 0 0(5) Z =『+2i 站如I 输入复数矩阵;§+6i 7 + 8i>> z=[1+2*i 3+4*i;5+6*i 7+8*i] 1.0000 + 2.0000i 3.0000 + 4.0000i 5.0000 + 6.0000i7.0000 + 8.0000i1 2川23 HI2.建立.m 文件，用for 循环语句生成5>5的矩阵A ：* ..■fa ri H *RF弓 6 III>> a=zeros(5); for m=1:5for n=1:5 a(m,n)=m+n-1 end end1 0 0 0 0a = 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0直翻转得到矩阵 B 和C 。