南华大学MATLAB实验报告1
matlab实验报告(实验1)
开课学院及实验室：机电学院 2012年11 月16日
（2）在E盘新建一个文件夹并把其加入搜索路径。
截屏相关过程结果并保存。
（3）逐步按照以下步骤进行并记录各个步骤的结果：
1）在指令窗中输入clear指令；
5）再次在指令窗中输入clear指令后，输入
8）把cow.m暂时移动到另一个无关的文件夹，如“我的文档”。
（2）保存全部变量为数据文件Mydata.mat；在对应文件夹中找到该文件，截屏后保存结果；
（4）删除全部内存变量；
Editor/Debugger和脚本编写初步
把以下指令编写为M脚本文件，并运行：
t=0:pi/50:4*pi;
y0=exp(-t/3); y=exp(-t/3).*sin(3*t); plot(t,y,'-r',t,y0,':b',t,-y0,':b') grid
什么？此方法也是查询函数注释的有效方法，请牢记。
MATLAB实验报告(打印版)
MATLAB实验报告班别：09电气3班姓名：李嘉明学号：200924122301实验一/二 MATLAB的基础操作3．自定义一个变量，然后分别用8种不同的数字显示格式显示查看。
>> a=[1.1,1.2,1.3;1.4,1.5,1.6;1.7,1.8,1.9]a =1.1000 1.2000 1.30001.4000 1.5000 1.60001.7000 1.8000 1.9000>> format long>> aa =1.100000000000000 1.200000000000000 1.3000000000000001.400000000000000 1.500000000000000 1.6000000000000001.700000000000000 1.800000000000000 1.900000000000000>> format short e>> aa =1.1000e+000 1.2000e+000 1.3000e+0001.4000e+000 1.5000e+000 1.6000e+0001.7000e+000 1.8000e+000 1.9000e+000>> format long e>> aa =1.100000000000000e+000 1.200000000000000e+000 1.300000000000000e+0001.400000000000000e+000 1.500000000000000e+000 1.600000000000000e+0001.700000000000000e+000 1.800000000000000e+000 1.900000000000000e+000 >> format hex>> aa =3ff199999999999a 3ff3333333333333 3ff4cccccccccccd3ff6666666666666 3ff8000000000000 3ff999999999999a3ffb333333333333 3ffccccccccccccd 3ffe666666666666>> format +>> aa =+++++++++>> format rat>> aa =11/10 6/5 13/107/5 3/2 8/517/10 9/5 19/10>> format short>> aa =1.1000 1.2000 1.30001.4000 1.5000 1.60001.7000 1.8000 1.90005、己举例比较给矩阵的全行赋值，取出矩阵中某一个或几个元素，然后组成一个新的矩阵。
MATLAB 语言及其应用 实验报告
《MATLAB 语言及其应用》实验报告目录实验一Matlab 使用方法和程序设计........................实验二控制系统的模型及其转换.............................实验三控制系统的时域、频域和根轨迹分析...........实验四动态仿真集成环境-Simulink.........................实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法；2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容：1、帮助命令使用help命令，查找 sqrt（开方）函数的使用方法；答：help sqrt2、矩阵运算（1）矩阵的乘法已知A=[1 2;3 4]; B=[5 5;7 8];求A^2*B解：>> A=[1 2;3 4];B=[5 5;7 8];A^2*Bans =105 115229 251（2）矩阵除法已知 A=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];A\B,A/B解：>> A=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];>> A\Bans =1.0e+016 *-0.4504 1.8014 -1.35110.9007 -3.6029 2.7022-0.4504 1.8014 -1.3511>> A/Bans =1.0000 1.0000 1.00004.0000 2.5000 2.00007.0000 4.0000 3.0000（3）矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i];求A.', A'解：>> A=[5+i,2-i,1;6*i,4,9-i];>> A.'ans =5.0000 + 1.0000i 0 +6.0000i2.0000 - 1.0000i 4.00001.0000 9.0000 - 1.0000i>> A'ans =5.0000 - 1.0000i 0 -6.0000i2.0000 + 1.0000i 4.00001.0000 9.0000 + 1.0000i（4）使用冒号表达式选出指定元素已知： A=[1 2 3;4 5 6;7 8 9];求A中第3列前2个元素；A中所有列第2，3行的元素；>> A=[1 2 3;4 5 6;7 8 9];>> B=A(1:2,3)B =36>> B1=A(2:3,:)B1 =4 5 67 8 9用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列>> A=magic(4)A =16 2 3 135 11 10 89 7 6 124 14 15 1>> B=A(:,1:3)B =16 2 35 11 109 7 64 14 153、多项式（1）求多项式4xxp的根2=x)(3-->> p=[1 0 -2 -4];>> r=roots(p)r =2.0000-1.0000 + 1.0000i-1.0000 - 1.0000i（2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征多项式;>> A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ;>> p=poly(A)p =1.0000 -6.9000 -77.2600 -86.1300 604.5500把矩阵A作为未知数代入到多项式中；>> A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ;>> syms x;P=x^3-2*x-4;>> P1=subs(P,x,A)P1 =-4.6720 17.0000 111.0000 -5.0710111.0000 -2.4870 111.0000 200.000017.0000 707.0000 -4.0000 -5.0000-5.0000 0 17.0000 52.00004、基本绘图命令（1）绘制余弦曲线 y=cos(t)，t∈[0，2π]（2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π]解：（1）程序为：t=[0:0.05:2*pi];y=cos(t);plot(t,y)（2）程序如下：x=[0:0.05:2*pi];y1=cos(x-0.25);y2=sin(x-0.5);plot(x,y1,x,y2)绘制曲线如下：5、基本绘图控制绘制[0，4π]区间上的x1=10sint曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线（3）标注控制：坐标轴名称、标题、相应文本；程序为：>> t=[0:0.05:4*pi];x1=10*sin(t);plot(t,x1,'r-.+')grid on;axis([0,15,-10,10]);title('曲线x1=10sint');xlabel('T轴');ylabel('X1轴');set(gca,'xminortick','on');set(gca,'yminortick','on');6、基本程序设计（1）编写命令文件：计算1+2+…+n<2000 时的最大n值；>> s=1;n=1;while(s<2000-n),n=n+1;s=s+n;end,[s,n]ans =1953 62（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。
MATLAB实验报告
实验一、Matlab 软件的安装与基本使用一、实验目的和要求：学会MA TLAB 软件的安装方法，熟悉MA TLAB 的工作环境及其基本功能的设置，以及基本指令的使用，了解MA TLAB 数值数组及其矩阵操作的使用、复数的输入方法。
二、实验内容和步骤：1、完成MA TLAB 软件的安装。
2、安装完成后启动MA TLAB 软件，熟悉MA TLAB 的工作环境。
打开MA TLAB 7，默认打开的窗口包括： (1)命令窗口(Command Window)； (2)命令历史窗口(Command History)； (3)工作间管理窗口(Workspace)；(4)当前路径窗口(Current Directory)此外，还有捷径键、编译窗口、图形窗口和帮助窗口等其他种类的窗口，结合课本实例，熟悉每个窗口的功能和基本使用。
3、（1）求23)]47(212[÷-⨯+的运算结果，注意命令后有无分号的作用。
命令及结果如下：t=[12+2*(7-4)]/3^2t = 2（2）求20012x x v t at =++的值，其中t=5s 、010x m =、015/v m s =、29.81/sec a m =-。
命令及结果如下：syms t x0 v0 a; t=5,x0=10,v0=15,a=-9.18; x=x0+v0*t+1/2*a*t^2 t = 5 x0 = 10 v0 = 15 x = -29.7500(3) 设3x =, 4v =，用MA TLAB 命令求下列各式的值：1）22log()x v + 命令及结果如下：syms x v x=3,v=4; log(x^2+v^2)x = 3 ans =3.21892)2(2)x v - 命令及结果如下：syms x vx=3,v=4;f=(sqrt(x-3))/(x-2*v)^2 x = 3 f = 0 3) |sin 2|vx e命令及结果如下： syms x v x=3,v=4;t=abs(2*x)*exp(v) x = 3t = 327.58894、用三种不同的方法输入矩阵1234246836912a ⎛⎫ ⎪= ⎪ ⎪⎝⎭，运行并观察实验结果，矩阵的三种基本输入方法如下：a=[1,2,3,4;3,4,6,8;3,6,9,12] a =1 2 3 4 3 4 6 8 3 6 9 12>> a=[1,2,3,43,4,6,8 3,6,9,12]a =1 2 3 4 3 4 6 8 3 6 9 12>> >>利用M 文件创建 1 2 3 4 3 4 6 8 3 6 9125、输入复数矩阵⎥⎦⎤⎢⎣⎡-+++=i ii i24356223B ，并求其实部、虚部、模和幅角（分别用弧度和角度表示），了解MA TLAB 向量化运算的特点。
Matlab实验报告
实验结果及分析实验1：程序如下x=1:10y=2*x;plot(x,y)仿真结果：实验结果分析：仿真结果是条很规则的直线，X轴和Y轴一一对应，清楚明了，而序又特别简单。
所以用Maltab 软件很方便地画出规则的直线，方便研究。
实验结果及分析1、A=2、A=1A=实验结果及分析实验三 Matlab在信号与系统中的应用实验名称实验1、掌握信号与系统课程中基本知识的Matlab编程、仿真方法目的实验原理实验1程序：b=[1];a=[1 1];p=;t=0:p:5;x=exp(-3*t);subplot(1,2,1);impulse(b,a,0:p:5);title('冲激响应');subplot(1,2,2);step(b,a,0:p:5);title('阶跃响应');实验内容<设计性实验>1、用MATLAB在时域中编程求解y′(t)+y(t)=f(t), f(t)= exp(-3t)ε(t)的冲激响应、阶跃响应。
在simulink仿真环境下，设计系统框图，分析系统的冲激响应、阶跃响应。
<设计性实验>（选做）2、用MATLAB在时域中编程求解y′(t)+y(t)=f(t), f(t)=(1+exp(-3t))ε(t)的冲激响应、阶跃响应，要求用conv编程实现系统响应。
在simulink仿真环境下，设计系统框图，分析系统的冲激响应、阶跃响应。
实验结果及分析实验1仿真结果：simulink仿真环境下冲激响应阶跃响应实验名称实验四 Matlab在数字信号处理中的应用实验结果及分析实验1仿真结果:6khz12kHZ。
MATLAB实验报告
MATLAB实验报告MATLAB实验报告姓名：专业：学号：实验⼀MATLAB环境的熟悉与基本运算⼀、实验⽬的：1．熟悉MATLAB开发环境2．掌握矩阵、变量、表达式的各种基本运算⼆、实验基本知识：1.熟悉MATLAB环境:MATLAB桌⾯和命令窗⼝、命令历史窗⼝、帮助信息浏览器、⼯作空间浏览器⽂件和搜索路径浏览器。
2.掌握MATLAB常⽤命令3.MATLAB变量与运算符变量命名规则如下：（1）变量名可以由英语字母、数字和下划线组成（2）变量名应以英⽂字母开头（3）长度不⼤于31个（4）区分⼤⼩写MATLAB中设置了⼀些特殊的变量与常量，列于下表。
MATLAB运算符，通过下⾯⼏个表来说明MATLAB的各种常⽤运算符表2MATLAB算术运算符表3MATLAB关系运算符表4MATLAB逻辑运算符表5MATLAB特殊运算4.MATLAB的⼀维、⼆维数组的寻访表6⼦数组访问与赋值常⽤的相关指令格式5.MATLAB的基本运算表7两种运算指令形式和实质内涵的异同表6.MATLAB的常⽤函数表8标准数组⽣成函数表9数组操作函数三、实验内容1、学习安装MATLAB软件。
2、学习使⽤help命令，例如在命令窗⼝输⼊helpeye，然后根据帮助说明，学习使⽤指令eye（其它不会⽤的指令，依照此⽅法类推）3、学习使⽤clc、clear，观察commandwindow、commandhistory和workspace等窗⼝的变化结果。
4、初步程序的编写练习，新建M-file，保存（⾃⼰设定⽂件名，例如exerc1、exerc2、exerc3……），学习使⽤MATLAB的基本运算符、数组寻访指令、标准数组⽣成函数和数组操作函数。
注意：每⼀次M-file的修改后，都要存盘。
四、实验结果练习A：（1）helprand，然后随机⽣成⼀个2×6的数组，观察commandwindow、commandhistory和workspace等窗⼝的变化结果。
matlab实验报告
（2）将上述信号左平移或右平移 得到两个新信号y2和y3，利用hold命令在同一图中显示三条曲线并通过线型和颜色加以区分。
（3）利用subplot命令将y，y2和y3以子图的形式显示。
3.
实验结果及分析
1．（1）
t1 =
0.0780
实验结果及分析
（2）
t2 =
0.0780
实验结果及分析
（3）
t3 =
0
实验结果及分析
2.
实验
名称
实验四Matlab在通信系统中的应用
实验目的
1、练习通信原理中各种调制方法及简单通信系统的Matlab编程实现；
2、学习Matlab与Simulink的混合编程。
（2）y变成虚数后，重复上述运算。
（3）创建矩阵 ，并进行如下操作，取出矩阵A的前两行元素，生成一个新矩阵B。生成一个3*3的单位矩阵B，与矩阵A进行加、减、乘、除、相等、不相等运算。
（4）随机生成2个复数x1和x2，计算x3=x1÷x2，x3的实部real_x3、虚部image_x3、模abs_x3、幅角angle_x3，并把x1，x2，x3，real_x3，image_x3，abs_x3，angle_x3保存至文件complex_学号.mat。利用save命令保存数据为mat格式文件后，用clear命令清除工作空间中的变量，再利用load命令从文件读入数据。
2、流程控制语句相关操作。
（1）使用for语句及while语句求解1至100整数中奇数的和。
（2）求连续自然数的和，当和大于等于1000时，求最后一个自然数以及自然数的和。
matlab实验报告总结
matlab实验报告总结1.求一份matlab的试验报告计算方法试验报告3【实验目的】检查各种数值计算方法的长期行为【内容】给定方程组x'(t)=ay(t),y'(t)=bx(t), x(0)=0, y(0)=b的解是x-y 平面上的一个椭圆，利用你已经知道的算法，取足够小的步长，计算上述方程的轨道，看看那种算法能够保持椭圆轨道不变。
（计算的时间步长要足够多）【实验设计】用一下四种方法来计算：1. Euler法2. 梯形法3. 4阶RK法4. 多步法Adams公式【实验过程】1. Euler法具体的代码如下：clear;a=2;b=1;A=[0 a; -b0];U=[];u(:,1)=[0;b];n=1000000;h=6*pi/n;fori=1:n delta(i)=((u(1,i)/a)^2+(u(2,i)/b)^2)^0.5; u(:,i+1)=u(:,i)+h*A*u(:,i);endt=1:n+1;subplot(1, 2,1);plot(1:n,delta);gridon;subplot(1,2,2);plot(u(1,:),u(2,:));gridon;max(abs(delta-ones(1,length(delta)）））；结果如下：2. 梯形法具体的代码如下：clear;a=2;b=1;A=[0 a; -b 0];U=[];u(:,1)=[0;b];n=300;h=6*pi/n;for i=1:n delta(i)=((u(1,i)/a)^2+(u(2,i)/b)^2)^0.5;v1=u(:,i)+h*A*u(:,i);v2=u(:,i)+h*A*(u(:,i)+v1)/2;1u(:,i+1)=u(:,i)+h*A*(u(:,i)+v2)/2;endt=1:n+1;sub plot(1,2,1);plot(1:n,delta);gridon;subplot(1,2,2)；结果如下 3. 4阶RK法clear;a=2;b=1;A=[0 a; -b 0];U=[];u(:,1)=[0;b];n=70;h=6*pi/n;for i=1:n delta(i)=((u(1,i)/a)^2+(u(2,i)/b)^2)^0.5;k1=A*u(:,i); k2=A*(u(:,i)+h/2*k2); k3=A*(u(:,i)+h*k3); k4=A*(u(:,i)+h*k3); u(:,i+1)=u(:,i)+h/6*(k1+2*k2+2*k3+k4);endt=1:n+1 ;subplot(1,2,1);plot(1:n,delta);gridon;subplot(1,2,2)；结果如下：4. 多步法Adams公式clear;a=2;b=1;A=[0 a; -b 0];U=[];u(:,1)=[0;b];n=200;h=6*pi/n;u(:;2)=u(u,1)+h*A*u(:,1);u(:;3)=u(u,2)+h/2*A*(3*u(:,2)-u(:,1));u(:;4)=u(u,3)+h/12*A*(23*u(:,3)-16*u(:,2)+5*u(:, 1)); delta(1)=((u(1,1)/a)^2+(u(2,1)/b^2)^0.5 delta(2)=((u(1,2)/a)^2+(u(2,2)/b^2)^0.5delta(3)=((u(1,3)/a)^2+(u(2,3)/b^2)^0.5for i=4:n delta(i)=((u(1,i)/a)^2+(u(2,i)/b)^2)^0.5;u(:,i+1)=u(:,i)+h/24*A*(55*u(:,i)-59*u(:,i-1)+37 *u(:,i-1)+37*u(:,i-2)-9*u(:,i-3));endt=1:n+1;sub plot(1,2,1);plot(1:n,delta);gridon;subplot(1,2,2)；结果如下：【实验分析】通过这几种方法对比，发现最为稳定的是多步法Adams公式和4阶RK法，其次是梯形法，而欧拉法最为不稳定。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
核科学技术学院实验报告实验项目名称MATLAB数值计算所属课程名称MATLAB及应用实验类型上机实验实验日期指导教师班级学号姓名成绩一、实验名称MATLAB数值计算二、实验目的（1）掌握MATLAB变量的使用（2）掌握MATLAB数组的创建（3）掌握MATLAB数组和矩阵的运算（4）熟悉MATLAB多项式的运用三、实验原理1. 矩阵分析矩阵转置：单引号（’）矩阵的旋转：rot90(A,k)，功能是将矩阵A逆时针旋转90度的k倍，缺省值是1矩阵的左右翻转：fliplr(A)矩阵的上下翻转：flipud(A)矩阵的逆：inv(A)，与A^(-1)等价矩阵的行列式：det(A)矩阵的秩：rank(A)矩阵的迹：trace(A)将矩阵化为最简式：rref(A)矩阵的特征值与特征向量：(1) E = eig(A)；矩阵A的所有特征值构成向量E；(2) [V,D]=eig(A)；A的所有特征值构成对角阵D，A的特征向量构成V的列向量；2. 多项式多项式的建立：若多的项的全部根构成的向量为X，则以X为根的多项式为poly(X)多项式的根：roots(p)计算以向量p为系数的多项式的根，包括重根，复根多项式求值：polyval(p,x)，p是多项式的系数，x可以是一个数也可以是一个矩阵多项式求拟合次数：polyfit(x,y,n)，x可以是一个数也可以是一个矩阵，y 是x对应的数或矩阵多项式的四则运算：(1)P1+P2；(2)P1-P2；(3)conv(P1,P2), (4)deconv(P1,P2)四、实验容1. 已知矩阵11 12 13 1421 22 23 24A= 31 32 33 3441 42 43 44(1) 如何输出A的第1列？(2) 如何输出A的第2行？(3) 如何输出A的第2列和第3列？(4) 要得到22233233⎛⎫⎪⎝⎭，怎么办？(5) 如何输出A的第1列和第3列？(6) 如何将21和31按一列输出？(7) 如何将A所有元素按从最左列至最右列新排列为一列输出？(8) 如何将A原阵输出？(9) 如何输出2阶全1阵？(10) 如何输出2阶单位阵？(11) 试用分块矩阵的方法生成阵11 12 13 14 1 121 22 23 24 1 131 32 33 34 1 041 42 43 44 0 1(12) 求A的主对角线(13) 求除去A阵第1列后新阵的主对角线(14) 求除去A阵第1行后新阵的主对角线(15) 求除去A 阵从左到右前2列后新阵的主对角线(16) 要分别单独输出A 阵元素11、22、33、44，怎么办？(17) 如何将11、31、12、32、13、33、14、34按顺序一行输出？2. 已知 2354127354217623A -⎡⎤⎢⎥--⎢⎥=⎢⎥-⎢⎥--⎣⎦.求（1）A 的逆；（2）A 的行列式；（3）A 的迹； （4）A 的所有特征向量和特征值。
3. A=magic(3)，B=7*rand(3)，计算数组A 、B 乘积，计算A&B ，A|B ，~A ，A==B ，A>B 如何生成一个三阶幻方、3阶随机阵、1行2列全零阵、3行4列全一阵、3阶单位阵？4. 生成一个4阶Hilbert 矩阵H ，（1）求H 的转置；（2）将H 旋转90度；（3）对H 实行左右翻转；（4）对H 实行上下翻转。
5 输入如下矩阵A0 π/3A ＝ π/6 π/2(1) 求矩阵B1, B1中每一元素为对应矩阵Ａ中每一元素的正弦函数(2) 求矩阵B2, B2中每一元素为对应矩阵Ａ中每一元素的余弦函数(3) 求 B 12+B 22(4) 求矩阵Ａ的特征值与特征矢量：称特征矢量为Ｍ，而特征值矩阵为Ｌ(5) 求Msin(L)M -1(6) 使用funm 命令求矩阵A 的正弦函数（结果应该与（5）同）(7) 求cosA(8) 证明 sin 2A+cos 2A=I I 是全1阵6 某专业有三名研究生,本学期选修了四门课程,若这些研究生的,学号,性别,出生年月,课程名称,考试成绩可任意假定,(1)分别用结构型变量和细胞型变量表示以上信息;(2)举例说明查阅以上任何一条信息的方法;(3)求每一个研究生的平均成绩.7. 已知多项式32()3279p x x x x =-++，试求：（1）p(x)的根；（2）由其根生成一个多项式q(x)并与p(x)比较；（3）计算p(1.5)，p(-2)，p(5)的值。
8 已知矩阵A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4],求(1)A的特征多项式(2)特征多项式中未知数为20 时的值(3)特征多项式的根(4)特征多项式的导数9. 在实验中测得如下10组数据:X 1 2 3 4 5 6 10 12 15 16 Y 7 15 19 30 38 37 23 69 39 28 （1）求最多能拟合多项式的次数是多少？并求出各项系数。
（2）将数据点和拟合曲线在同一图中绘出。
五、实验过程及结果（含源代码）1、在MA TLAB中新建.m文件，输入实验代码如下并保存：clc,clearA=[11 12 13 14;21 22 23 24;31 32 33 34;41 42 43 44];A1=A(:,1)A2=A(2,:)A3=A(:,2:3)A4=A(2:3,2:3)A5=A(:,1:2:3)A6=A(2:3)A7=A(:)A8=A(:,:)A9=ones(2,2)A10=eye(2)A11=[A,[ones(2,2);eye(2)]]A12=diag(A)A13=diag(A,1)A14=diag(A,-1)A15=diag(A,2)点击“run”运行得到结果如下：A1 =11213141A2 =21 22 23 24A3 =12 1322 2332 3342 43A4 =22 2332 3311 1321 2331 3341 43 A6 =21 31 A7 =112131411222324213233343142444A8 =11 12 13 1421 22 23 2431 32 33 3441 42 43 44A9 =1 11 1A10 =1 00 1A11 =11 12 13 14 1 121 22 23 24 1 131 32 33 34 1 041 42 43 44 0 1 A12 =112233A13 =122334A14 =213243A15 =13242、在命令窗口输入代码如下：A=[2 3 -5 4;-1 2 7 -3;5 -4 2 1;7 6 -2 -3];inv(A)ans =0.0198 0.0010 0.1069 0.06100.1199 0.1116 -0.0636 0.02710.0709 0.1616 0.0667 -0.04480.2388 0.1178 0.0777 -0.1069 det(A)-1918trace(A)ans =3E=eig(A)E =-5.65454.68881.9828 + 8.2711i1.9828 - 8.2711i[V,D]=eig(A)V =0.3850 0.4655 -0.0889 + 0.4794i -0.0889 - 0.4794i-0.1196 0.4083 -0.1447 - 0.5624i -0.1447 + 0.5624i-0.1972 0.4872 0.5903 0.5903-0.8936 0.6157 -0.1444 + 0.2359i -0.1444 - 0.2359iD =-5.6545 0 0 00 4.6888 0 00 0 1.9828 + 8.2711i 00 0 01.9828 - 8.2711i3、在命令窗口键入实验代码如下：>> A=magic(3),B=7*rand(3)A =8 1 63 5 74 9 2B =6.6509 3.4019 3.19531.6180 6.2391 0.12954.24795.3347 5.7499>> A*Bans =80.3126 65.4622 60.190857.7778 78.7438 50.482449.6611 80.4287 25.4465 >> A&Bans =1 1 11 1 11 1 1>> A|Bans =1 1 11 1 11 1 1>> ~Aans =0 0 00 0 00 0 0>> A==Bans =0 0 00 0 00 0 0>> A>Bans =1 0 11 0 10 1 0>> A=magic(3)A =8 1 63 5 74 9 2>> B=rand(3)B =0.9501 0.4860 0.45650.2311 0.8913 0.01850.6068 0.7621 0.8214 >> C=zeros(1,2)C =0 0>> D=ones(3,4)D =1 1 1 11 1 1 11 1 1 1>> E=eye(3)E =1 0 00 1 00 0 14、在命令窗口键入实验代码及运行结果如下：>> H=hilb(4)H =1.0000 0.5000 0.3333 0.25000.5000 0.3333 0.2500 0.20000.3333 0.2500 0.2000 0.16670.2500 0.2000 0.1667 0.1429 >> H'ans =1.0000 0.5000 0.3333 0.25000.5000 0.3333 0.2500 0.20000.3333 0.2500 0.2000 0.16670.2500 0.2000 0.1667 0.1429 >> rot90(H)ans =0.2500 0.2000 0.1667 0.14290.3333 0.2500 0.2000 0.16670.5000 0.3333 0.2500 0.20001.0000 0.5000 0.3333 0.2500 >> fliplr(H)ans =0.2500 0.3333 0.5000 1.00000.2000 0.2500 0.3333 0.50000.1667 0.2000 0.2500 0.33330.1429 0.1667 0.2000 0.2500>> flipud(H)ans =0.2500 0.2000 0.1667 0.14290.3333 0.2500 0.2000 0.16670.5000 0.3333 0.2500 0.20001.0000 0.5000 0.3333 0.2500 5、在命令窗口键入实验代码及运行结果如下：>> A=[0 pi/3;pi/6 pi/2]A =0 1.04720.5236 1.5708>> B1=sin(A)B1 =0 0.86600.5000 1.0000>> B2=cos(A)B2 =1.0000 0.50000.8660 0.0000>> B1^2+B2^2ans =1.8660 1.3660>> [M,L]=eig(A)M =-0.9628 -0.48960.2703 -0.8719L =-0.2940 00 1.8648>> M*sin(L)*inv(M)ans =-0.1200 0.60480.3024 0.7873>> funm(A,sin)ans =-0.1200 0.60480.3024 0.7873>> funm(A,cos)ans =0.7873 -0.6048-0.3024 -0.1200>> (funm(A,sin))^2+(funm(A,cos))^2 ans =-0.0000 1.00006、实验代码及运行结果如下：(1) 分别用细胞型变量和结构型变量表示题示信息如下：% 细胞型变量A1={'明',20130001,'男',19890214,'MATLAB',99,'ENGLISH',88,'Matrix',96,'Math',94};A2={'远',20130006,'男',19900425,'MATLAB',98,'ENGLISH',80,'Matrix',90,'Math',96};A3={'王云',20130009,'女',19901001,'MATLAB',97,'ENGLISH',89,'Matrix',88,'Math',90};% 结构型变量B1=struct('Name','明','StuNumber',20130001,'Sex','男','Birth',19890214,...'MATLAB',99,'ENGLISH',88,'Matrix',96,'Math',94);B2=struct('Name','远','StuNumber', 20130006,'Sex','男','Birth', 19900425,...'MATLAB',98,'ENGLISH',80,'Matrix',90,'Math',96);B3=struct('Name','王云','StuNumber', 20130009,'Sex','女','Birth', 19901001,...'MATLAB',97,'ENGLISH',90,'Matrix',88,'Math',90);C=[B1 B2 B3];(2) 查询信息，如查询明的学号，实验代码如下：name='明'information='StuNumber';for i=1:3if strcmp(C(i).Name,name)==1disp(information);getfield(C(i),information)break;endend运行结果如下：name =明StuNumberans =20130001(3)求每位研究生的平均成绩：for i=1:3Average(i)=(C(i).MATLAB+C(i).ENGLISH+C(i).Matrix+C(i).Math)/4; endAverage结果如下：Average =94.2500 91.0000 91.25007、在命令窗口键入实验代码及运行结果如下：>> p=[3 -2 7 9];>> x=roots(p)x =0.7513 + 1.7392i0.7513 - 1.7392i-0.8358>> q=poly(x)q =1.0000 -0.66672.33333.0000>> b=[1.5 -2 5];>> polyval(p,b)ans =25.1250 -37.0000 369.00008、在命令窗口键入实验代码及运行结果如下：>> A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4];>> p=poly(A);>> poly2sym(p)ans =x^4 - (69*x^3)/10 - (3863*x^2)/50 - (8613*x)/100 + 12091/20>> polyval(p,20)ans =7.2778e+04>> roots(p)ans =13.0527-4.1671 + 1.9663i-4.1671 - 1.9663i2.1815>> Dp=polyder(p);>> poly2sym(Dp)ans =4*x^3 - (207*x^2)/10 - (3863*x)/25 - 8613/1009、因题目给出的实验数据是10组，可以看做是10个多项式方程组，仅当多项式的次数小于方程组的个数时，方程组为超定方程，此时才可以进行最小二乘多项式拟合，因此，最多拟合次数为9次，实验代码及运行结果如下：X=[1 2 3 4 5 6 10 12 15 16];Y=[7 15 19 30 38 37 23 69 39 28];a=polyfit(X,Y,9);x2=1:0.1:16;y2=polyval(a,x2);plot(X,Y,'r*',x2,y2,'b-')legend('原数据点','拟合曲线')拟合多项式系数如下：a =Columns 1 through 70.0000 -0.0005 0.0150 -0.2693 2.9713 -20.5126 85.9563Columns 8 through 10-204.4409 250.3796 -107.0989实验图形如下所示：（此处按题号顺序写上所写的程序语句容以及matlab中出现的程序运行结果）六、实验总结。