Matlab作业第1-4题

合集下载

MATLAB实验报告(1-4)

信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。

2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。

4.学会运用MATLAB进行连续信号时移、反折和尺度变换。

5.学会运用MATLAB进行连续时间微分、积分运算。

6.学会运用MATLAB进行连续信号相加、相乘运算。

7.学会运用MATLAB进行连续信号的奇偶分解。

二、实验任务将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。

三、实验内容1.MATLAB软件基本运算入门。

1). MATLAB软件的数值计算：算数运算向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。

2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn 为结束值。

矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开；矩阵的不同行之间必须用分号”;”或者ENTER分开。

2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。

3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。

举例：计算一个函数并绘制出在对应区间上对应的值。

2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y)2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’)6).输出：grid on举例1：举例2：3.matlab程序流程控制1).for循环：for循环变量=初值:增量:终值循环体End2).while循环结构：while 逻辑表达式循环体End3).If分支：(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号：y=k*exp(a*t)举例：2).正弦信号：y=k*sin(w*t+phi)3).复指数信号：举例：4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号：y=tripuls(t,width,skew)(skew的取值在-1~+1之间，若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算：语句格式：d iff(function,’variable’,n)Function:需要进行求导运算的函数，variable:求导运算的独立变量，n：求导阶数2).连续时间信号的积分运算：语句格式：int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件，并且能够输入简单的语句并输出相应的结果和波形图，也在一定程度上巩固了c语言的一些语法。

matlab编程题整理

作业一1输出x,y两个中值较大的一个值x=input(‘x’);y=input(‘y’);if x>yxelseyend2输入x，计算y的值。

计算函数的值y=x+1,x<0,y=2x-1,x≧0x=input错误！未指定书签。

(‘x);if x<0y=x+1elsey=2*x-1end3输入一学生成绩，评定其等级，方法是：90~100分为“优秀”，80~89分为“良好”，70~79分为“中等”，60~69分为“及格”，60分为“不合格”x=input(‘x’)if x>100|x<0y=’输入错误’elseif x>=90y=’优秀’elseif x>=80y=’良好’elseif x>=70y=’中等’elseif x>=60y=’及格’elsey=’不合格’emd4某超市节日期间举办购物打折的促销活动，优惠办法是：每位顾客当天一次性购物在100元以上者，按九五折优惠；在200元以上者，按九折优惠；在300元以上者，按八五折优惠；在500元以上者，按八折优惠。

x=input(‘x’);if x>=500y=x*0.8elseif x>=300y=x*0.85elseif x>=200y=x*0.9elseif x>=100y=x*0.95else y=xend 5编程计算：s=1+2+3+…+100sum=0;for i=1:100sum=sum+i;endsum引申1!-2!+3!-4!+5!- (99)sum=0;for i=1:99pdr=1;for k=1:ipdr=pdr*k;endsum=sum+pdr*(-1)^(i-1);endsum引申1*2*3*4*……*100sum=1；for i=1:100sum=sum*iendsum6计算1~100的奇数和sum=0;for i=1:2:100sum=sum+i;endsum7百元买百鸡问题。

Matlab优化设计作业答案

长江大学机械工程学院机械优化设计大作业班级2012年5月31-361、⎩⎨⎧=+-=++⋅+-++=22422min 321321321232221x x x x x x t s x x x x x x f解： function f=fun1(x)f=x(1)^2+2*x(2)^2+x(3)^2-2*x(1)*x(2)+x(3)clearclcx0=[1;1;1];A=[];b=[];Aeq=[1 1 1;2 -1 1];beq=[4;2];lb=[];ub=[];[x,fval,exitflag]=fmincon('fun1',x0,A,b,Aeq,beq,lb,ub)x =1.90911.95450.1364fval =3.9773exitflag =12、221212min ()2130f x x x x s t x =+-+⋅-≤解： function f=fun2(x)f=x(1)^2+x(2)^2-2*x(1)+1clearclcx0=[1;1];A=[0 -1];b=-3;Aeq=[];beq=[];lb=[];ub=[];[x,fval,exitflag]=fmincon('fun2',x0,A,b,Aeq,beq,lb,ub)x =13fval =9exitflag =13、 2212min ()(2)f X x x =-+112222312..()0()0()10s t g X x g X x g X x x =-≤=-≤=-+-≤解：function f=fun3(x)f=(x(1)-2)^2+x(2)^2function [g,ceq]=mycon3(x)g=[-x(1);-x(2);-x(1)^2+x(2)^2-1];ceq=[];clearclcx0=[0;0];A=[];b=[];Aeq=[];beq=[];lb=[];ub=[];[x,fval]=fmincon('fun3',x0,A,b,Aeq,beq,lb,ub,'mycon3')x =2.0000fval =2.2204e-0164、⎩⎨⎧≤--≤+--⋅++++=01005.1)12424(min 21212122122211x x x x x x t s x x x x x e f x解： function f=fun4(x)f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1);function [g,ceq]=mycon4(x)g(1)=1.5+x(1)*x(2)-x(1)-x(2);g(2)=-x(1)*x(2)-10;ceq=[];clearclcx0=[0;0];A=[];b=[];Aeq=[];beq=[];lb=[];ub=[];[x,fval]=fmincon('fun4',x0,A,b,Aeq,beq,lb,ub,'mycon4')x =1.1825-1.7398fval =3.06085 喜糖问题：需要购买甲乙两种喜糖，喜糖甲10元/斤，喜糖乙20元/斤。

MATLAB作业4-1

MATLAB作业4-1MATLAB 作业41、⽤2sin(103)y t =+在（0，3）区间内⽣成⼀组较稀疏的数据，并⽤⼀维数据插值的⽅法对给出的数据进⾏曲线拟合，并将结果与理论曲线相⽐较。

解：2、⽤242231(,)sin()3x y f x y e xy x y x y--=++原型函数⽣成⼀组⽹络数据或随机数据，分别拟合出曲⾯，并和原曲⾯进⾏⽐较。

解：⽹格数据：随机数据：x∈-区间内的光滑函数曲线，⽐较各种插值算法3、假设已知⼀组数据，试⽤插值⽅法绘制出(2,4.9)>> x=[-2,-1.7,-1.4,-1.1,-0.8,-0.5,-0.2,0.1,0.4,0.7,1,1.3,...1.6,1.9,2.2,2.5,2.8,3.1,3.4,3.7,4,4.3,4.6,4.9];y=[0.10289,0.11741,0.13158,0.14483,0.15656,0.16622,0.17332,...0.1775,0.17853,0.17635,0.17109,0.16302,0.15255,0.1402,...0.12655,0.11219,0.09768,0.08353,0.07019,0.05786,0.04687,...0.03729,0.02914,0.02236];解：x y在(0.1,0.1)~(1.1,1.1)区域内的点进⾏插值，并⽤三维曲⾯4、假设已知实测数据由下表给出，试对(,)>> [x,y]=meshgrid(0.1:0.1:1.1);z=[0.83041,0.82727,0.82406,0.82098,0.81824,0.8161,0.81481,0.81463,0.81579,0.81853,0.82304 ;0.83172,0.83249,0.83584,0.84201,0.85125,0.86376,0.87975,0.89935,0.92263,0.94959,0.9801; 0.83587,0.84345,0.85631,0.87466,0.89867,0.9284,0.96377,1.0045,1.0502,1.1,1.1529;0.84286,0.86013,0.88537,0.91865,0.95985,1.0086,1.0642,1.1253,1.1904,1.257,1.3222;0.85268,0.88251,0.92286,0.97346,1.0336,1.1019,1.1764,1.254,1.3308,1.4017,1.4605;0.86532,0.91049,0.96847,1.0383,1.118,1.2046,1.2937,1.3793,1.4539,1.5086,1.5335;0.88078,0.94396,1.0217,1.1118,1.2102,1.311,1.4063,1.4859,1.5377,1.5484,1.5052;0.89904,0.98276,1.082,1.1922,1.3061,1.4138,1.5021,1.5555,1.5573,1.4915,1.346;0.92006,1.0266,1.1482,1.2768,1.4005,1.5034,1.5661,1.5678,1.4889,1.3156,1.0454;0.94381,1.0752,1.2191,1.3624,1.4866,1.5684,1.5821,1.5032,1.315,1.0155,0.62477;0.97023,1.1279,1.2929,1.4448,1.5564,1.5964,1.5341,1.3473,1.0321,0.61268,0.14763];解：5、习题3和4给出的数据分别为⼀元数据和⼆元数据，试⽤分段三次样条函数和B样条函数对其进⾏拟合。

MATLAB作业

MATLAB作业⼀、必答题：1. MATLAB系统由那些部分组成？答：MATLAB系统主要由开发环境、MATLAB语⾔、MATLAB数学函数库、图形功能和应⽤程序接⼝五个部分组成。

2. 如何启动M⽂件编辑/调试器？答：在操作界⾯上选择“建⽴新⽂件”或“打开⽂件”操作时，M⽂件编辑/调试器将被启动。

在命令窗⼝中键⼊“edit”命令也可以启动M⽂件编辑/调试器。

3. 存储在⼯作空间中的数组能编辑吗？如何操作？答：存储在⼯作空间的数组可以通过数组编辑器进⾏编辑：在⼯作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输⼊修改内容即可。

4. 在MATLAB中有⼏种获得帮助的途径？答：在MATLAB中有多种获得帮助的途径：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器；（2）help命令：在命令窗⼝键⼊“help” 命令可以列出帮助主题，键⼊“help 函数名”可以得到指定函数的在线帮助信息；（3）lookfor命令：在命令窗⼝键⼊“lookfor 关键词”可以搜索出⼀系列与给定关键词相关的命令和函数（4）模糊查询：输⼊命令的前⼏个字母，然后按Tab键，就可以列出所有以这⼏个字母开始的命令和函数。

5. 有⼏种建⽴矩阵的⽅法？各有什么优点？答：（1）以直接列出元素的形式输⼊；（2）通过语句和函数产⽣；（3）.在m⽂件中创建矩阵；（4）从外部的数据⽂件中装⼊。

6. 命令⽂件与函数⽂件的主要区别是什么？答：命令⽂件: M⽂件中最简单的⼀种，不需输出输⼊参数，⽤M ⽂件可以控制⼯作空间的所有数据。

运⾏过程中产⽣的变量都是全局变量。

运⾏⼀个命令⽂件等价于从命令窗⼝中顺序运⾏⽂件⾥的命令，程序不需要预先定义，只要依次将命令编辑在命令⽂件中即可。

函数⽂件：如果M⽂件的第⼀个可执⾏⾏以function开始，便是函数⽂件，每⼀个函数⽂件定义⼀个函数。

数学建模作业题+答案

数学建模MATLAB 语言及应用上机作业11. 在matlab 中建立一个矩阵135792468101234501234A ⎡⎤⎢⎥⎢⎥=⎢⎥-----⎢⎥⎣⎦答案：A = [1,3,5,7,9;2,4,6,8,10;-1,-2,-3,-4,-5;0,1,2,3,4]2. 试着利用matlab 求解出下列方程的解（线性代数22页例14）123412423412342583692254760x x x x x x x x x x x x x x +-+=⎧⎪--=⎪⎨-+=-⎪⎪+-+=⎩ 答案：A=[2 ,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6]; B=[8;9;-5;0]; X=A\B 或A=[2,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6] b=[8,9,-5,0]' X=inv(A)*b3. 生成一个5阶服从标准正态分布的随机方阵，并计算出其行列式的值，逆矩阵以及转置矩阵。

答案:A=randn(5) det(A) inv(A) A'4. 利用matlab 求解出110430002A -⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦的特征值和特征向量。

答案：A=[-1,1,0;-4,3,0;0,0,2] [V,D]=eig(A)5.画出衰减振荡曲线3sin3t y et -=在[0,4]π上的图像。

要求，画线颜色调整为黑色，画布底面为白色。

（在实际中，很多打印机时黑白的，因此大多数作图要考虑黑白打印机的效果。

）给出恰当的x ，y 坐标轴标题，图像x 轴的最大值为4π。

6. 生成一个0-1分布的具有10个元素的随机向量，试着编写程序挑选出向量中大于0.5的元素。

数学建模和Matlab 上机作业2（2016-9-20）跟老师做（不用整合进作业中）：上机演示讲解：函数，递归的两个例子的写法。

附：1. Fibonacci Sequence （斐波那契数列）在数学上，费波那西数列是以递归的方法来定义： F1= 1；F2= 1；F （n ）=F （n-1）+F （n-2） 2. 阶乘举例：数学描述：n!=1×2×……×n ；计算机描述：n!=n*(n-1)!自己做（需要整合进作业中，提交到系统中）：1. 写一个m 文件完成分值百分制到5分制的转换（即输入一个百分制，转换后输出一个5级对应的得分，联系条件控制语句）。

MATLAB平时作业(图文版)

MATLAB 平时作业第一章习题16. 以下两种说法对吗？（1）“MATLAB 的数值表达精度与其指令窗中的数据显示精度相同。

”答：此种说法错误。

MATLAB 提供了控制数据显示格式的控制指format ，该指令并不改变MATLAB 内存中变量的精度，只是改变其显示精度。

（2）“MATLAB 指令窗中显示的数值有效位数不超过7位。

”答：此种说法错误。

当变量小于1000时，使用format 或format short 后，或者默认情况下，变量的显示精度最多不超过7位，但显示精度不等于变量的精度。

7. 想要在MATLAB 中产生二维数组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321S ，下面哪些指令能实现目的？ S=[1,2,3;4,5,6;7,8;9]S=[1 2 3;4 5 6;7 8 9]S=[1，2，3；4，5，6；7，8，9] %整个指令在中文状态下输入答：操作如图：第1、2条指令可以实现，第3条指令不可实现。

第三章习题31．在MATLAB 中，先运行指令A=magic(3), B=[1,2,1;3,4,3;5,6,7], C=reshape(1:6,3,2)生成阵列33⨯A ，23⨯B ，23⨯C ，然后根据运行结果回答以下问题：运行结果如图：（1）计算A*B, B*A ，这两个乘积相同吗？计算结果如图：答：不同。

（2）计算A\B, B/A ，左除、右除结果相同吗？计算结果如图：答：不同。

（3）计算B( : ,[1,2]).*C和C.*B( : , [1,2])，这两个乘积相同吗？计算结果如图答：相同。

（4）计算A\A和A.\A，这两个计算结果相同吗？计算结果如图：答：相同。

（5）计算A\eye(3)和inv(A)，这两个计算结果相同吗？计算结果如图：答：不同。

（提示：根据对计算结果的目测回答问题）2．在MATLAB中，先运行A=[1, 2; 3, 4],b=0.5,C=[4, 2; 1, 0.5], 然后根据计算结果回答以下问题：创建数据步骤略（1）计算A^b和A.^b, 这两个计算结果相同吗？答：不同。

matlab章课后作业

matlab章课后作业MATLAB基础教程1~8章作业Matlab第⼀章1.阐述Matlab的功能Matlab作为⼀种⾼级计算软件，是进⾏算法开发、数据可视化、数据分析以及数值计算的交互式应⽤开发环境，已被⼴泛应⽤于不同领域。

Matlab的基本功能包括：数学计算功能、图形化显⽰功能、M 语⾔编程功能、编译功能、图形⽤户界⾯开发功能、Simulink建模仿真功能、⾃动代码⽣成功能。

Matlab第⼆章1．创建double的变量，并进⾏计算。

(1)a=87，b=190，计算a+b、a-b、a*b。

(2)创建uint8 类型的变量，数值与(1)中相同，进⾏相同的计算。

>> a=87,b=190a =87b =190>> a+bans =277>> a-bans =-103>> a*bans =16530>> c=uint8(87), d=uint8(190)c =87d =190>> c+dans =255>> c-dans =>> c*dans =2552.计算(1)sin(60)(2)e^3(3)cos(3π/4)>> sind(60)ans =0.8660>> exp(3)ans =20.0855>> cos(3*pi/4)ans =-0.70713.设u=2，v=3，计算：(1)(2)(3)>> u=2;>> v=3;>> 4*u*v/log(v)ans =21.8457>> (exp(u)+v)^2/(v^2-u) ans = 15.4189>> sqrt(u-3*v)/(u*v) ans =0 + 0.4410i 4.计算如下表达式：(1)(2)>> (3-5*i)*(4+2*i)ans =22.0000 -14.0000i >> sin(2-8*i) 1.3553e+003 +6.2026e+002i5.判断下⾯语句的运算结果。

MATLAB作业4参考答案

MATLAB作业4参考答案MATLAB作业四参考答案2yt,，sin(103)1、用在(0，3)区间内生成一组较稀疏的数据，并用一维数据插值的方法对给出的数据进行曲线拟合，并将结果与理论曲线相比较。

【求解】类似于上面的例子～可以用几乎一致的语句得出样本数据和插值效果。

>> t=0:0.2:3;y=sin(10*t.^2+3); plot(t,y,'o')ezplot('sin(10*t^2+3)',[0,3]); hold on x1=0:0.001:3;y1=interp1(t,y,x1,'spline'); plot(x1,y1)由于曲线本身变换太大～所以在目前选定的样本点下是不可能得出理想插值效果的～因为样本数据提供的信息量不够。

为了得到好的插值效果～必须增大样本数据的信息量～对本例来说～必须在快变化区域减小样本点的步长。

>> hold offt=[0:0.1:1,1.1:0.04:3]; y=sin(10*t.^2+3); plot(t,y,'o')ezplot('sin(10*t^2+3)',[0,3]); hold on x1=0:0.001:3;y1=interp1(t,y,x1,'spline'); plot(x1,y1)241,,xy22fxyexyxy,，(,)sin()2、用原型函数生成一组网络数据或随机数据，分别拟合出曲3xy，3面，并和原曲面进行比较。

【求解】由下面的语句可以直接生成一组网格数据～用下面语句还可以还绘制出给定样本点是三维表面图。

>> [x,y]=meshgrid(0.2:0.2:2);z=exp(-x.^2-y.^4).*sin(x.*y.^2+x.^2.*y)./(3*x.^3+y); surf(x,y,z) 选择新的密集网格～则可以通过二元插值得出插值曲面。

Matlab 作业1

喻晓磊 200731510103 电子科学与技术习题一：编写M文件，计算以下问题：口口口口×口＝口口口口以上9个口代表1～9这9个数字，不得遗漏或重复要求：给出解算思路和M文件代码注意算法的效率解题理念：1，设以上方框中的数字分别是1~9之间的数字a、b、c、d、e、f、g、h、i。

2，我们可以推出，个位数e肯定不是1或者5 ，否则两个四位数的个位必然相等，即d=i；再者，其不可能等于9，由不重复的1~9组成的四位数，最大为9876，最小为1234，故e的最大可能取值为9876/1234=8.003……，即最大取8 。

3，我们看到，e的最小取值为2，那么a肯定不能大于或等于5 ，否则乘法运算后的结果为5位数，所以a 只能取1~4。

4，在此基础上，我们采取“穷举法”，逐个试验，从 a 开始，每一位与前面出现位的数字不相等，然后检测其是否满足算式。

程序代码如下：global aglobal bglobal cglobal dglobal eglobal fglobal gglobal hglobal ifor a=1:4for b=1:9if b~=a;for c=1:9if (c~=b)&&(c~=a);for d=1:9if (d~=b)&&(d~=a)&&(d~=c);for e=2:8 e~=5;if (e~=a)&&(e~=b)&&(e~=c)&&(e~=d);for f=1:9if (f~=a)&&(f~=b)&&(f~=c)&&(f~=d)&&(f~=e);for g=1:9if(g~=a)&&(g~=b)&&(g~=c)&&(g~=d)&&(g~=e)&&(g~=f);for h=1:9if(h~=a)&&(h~=b)&&(h~=c)&&(h~=d)&&(h~=e)&&(h~=f)&&(h~=g);for i=1:9if(i~=a)&&(i~=b)&&(i~=c)&&(i~=d)&&(i~=e)&&(i~=f)&&(i~=g)&&(i~=h);x=1000*a+100*b+10*c+d;y=e*x;z=1000*f+100*g+10*h+i;if y==z;r=xs=et=yend;end;******end;end;得到的结果如下：我们看到，结果有两组：1738 * 4= 69521963 * 4= 7852这两组结果都显示到了MATLAB软件的命令窗（COMMAND WINDOW）中，我们是使用了一个小技巧：将x、e、y的值赋给r、s、t时并没有在句末加分号，所以计算的过程也显示出来了，但是看工作区（work space）中，r、s、t 的值却只对应后面一组，这是因为，当循环进行到算出第二组结果时，便替换掉了第一组的值。

MATLAB课程作业练习题

内蒙古工业大学电力学院《MATLAB在电力系统中的应用》练习题汇总编者：顾秀芳单位：电力学院电力工程系时间：2012年8月实验一：《MATLAB操作环境及基本操作》实验指导书1、实验目的：了解MATLAB的运行环境，界面功能；熟悉MATLAB编程的基本方法；掌握MATLAB的基本画图功能；熟悉Simulink模块的基本操作及子系统的封装方法。

2、实验内容：操作MATLAB界面的主要菜单并了解其中的功能；练习矩阵的运算；编写或读懂并运行MATLAB程序；练习并实践MATLAB的画图功能的实现方法。

练习Simulink模块的基本操作；熟悉并练习子系统的封装方法；3、实验要求：通过预习，了解MATLAB的界面的组成，各界面实现的主要功能；了解MATLAB实现运算、逻辑判断的基本语法、函数的含义等；掌握MATLAB基本绘图功能，并实现简单的绘图；熟悉Simulink的及子系统封装的主要功能及基本操作。

4、实验步骤：（1）熟悉MATLAB界面的基本功能；（2）熟悉编程基础，并读通指定程序的含义，运行结果，并分析；（3）熟悉MATLAB的基本绘图语句，读通指定并读通指定程序的含义，运行结果并分析；实现给定的绘图任务；（4）熟悉Simulink模块的操作，并完成给定模型的搭建；（5）熟悉Simulink子系统模型搭建的步骤，并完成给定子系统的搭建及参数设置。

6 、实验题目：1）MATLAB环境的了解包括：窗口菜单的主要功能；各窗口实现的功能，包括：命令窗口，workspace/history,command history三个窗口的应用。

2）MATLAB编程基础内容包括：矩阵的形成及运算（矩阵相乘、点乘）；逻辑运算；函数的应用。

练习内容一：矩阵的形成及操作；1：形成矩阵2:从矩阵中提取子矩阵：3：矩阵A中元素A41的值；4：修改矩阵A的元素，使A41=3；5：矩阵A中最后2行和最后3行交汇的子矩阵的值；练习内容二：矩阵的运算1：矩阵定义为：求出下列矩阵运算的结果，并解释部分操作无法实现的原因：（1） R=A+B;（2） R=A*D;（3） R=A.*D;（4） R=A*C;（5） R=A.*C;（6） R=A\B;（7） R=A.\B;（8） R=A.^B.1.网络如图，图中标注为支路导纳yij.试通过程序完成Y导纳矩阵的形成，j含对地支路练习内容三：编程语句的理解1:运行书P20中例题，比较标量循环和矩阵运算的差异；比较程序：1.一个窗口画一个图t=0:0.1:2*pi;y1=sin(t);plot(t,y1,'r:')axis([0,2*pi,-2,2]);xlabel('huduzhi');ylabel('hanshuzhi');title('正弦函数和余弦函数');grid ongtext({'y1=sin(t)'});2.一个窗口绘两个图(1)一条plot语句画t=0:0.1:2*pi;y1=sin(t); y2=cos(t);plot(t,y1,'r:', t,y2,'b-.')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y1,y2');title('正弦函数和余弦函数');grid ongtext({'y1=sin(t)'; 'y2=cos(t)'});(1)分多条plot语句画t=0:0.1:2*pi;y1=sin(t); y2=cos(t); y3=tan(t) plot(t,y1,'r-')holdplot(t,y2,'b--')plot(t,y3,'k-.')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y1,y2,y3');grid ongtext({'y1=sin(t)'; 'y2=sin(t)';'y3=sin(t)'}); legend('y1','y2','y3')3、一个窗口分割成三个小窗口t=0:0.1:2*pi;y1=sin(t);subplot(3,1,1)plot(t,y1,'r:')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y1');title('正弦函数');grid ongtext('y1=sin(t)');subplot(3,1,2)plot(t,cos(t),'r:')axis([0,2*pi,-2,2]);xlabel('t');ylabel('cost');title('余弦函数');grid ongtext('y2=sin(t)');subplot(3,1,3)plot(t,tan(t),'r:')axis([0,2*pi,-2,2]);xlabel('t');ylabel('tant');title('正切函数');grid ongtext('y3=sin(t)');作业：t=0:0.1:2*pi;%y1=sin(t); y2=cos(t); y3=tan(t);y4=cot(t) subplot(2,2,1)plot(t,y1,'r-')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y1'); subplot(2,2,2)plot(t,y2,'b--*')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y2'); subplot(2,2,3)plot(t,y3,'k:')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y3'); subplot(2,2,4)plot(t,y4,'m-.',t,y1.*y4,'m-.+') axis([0,2*pi,-2,2]);xlabel('t');ylabel('y4');grid ongtext({'y4=cot(t)'; 'y1*y4'});legend('y4','y1*y4')2、编制程序或其它途径实现下列绘图功能结合上面矩阵的运算及绘图命令读书中P149\152\154程序第2章：了解Simulink的基本操作；熟悉simulink及SPS库中模块，练习搭建：1）P45图2-22a)2）P102图5-3，课后完成P5-29， P111图5-13，3）搭建下下题表达的模型3-5、考虑如图所示感应电机的等效电路，输入的交流电压源为220V、50Hz，其它参数为。

MATLAB第一章作业答案

第一章M A T L A B概况与基本操作1.选择题：(1)最初的MATLAB核心程序是采用A语言编写的。

(2)即将于2011年9月发布的MATLAB新版本的编号为D。

2011Ra 2011Rb R2011a R2011b(3)在默认设置中，MATLAB中的注释语句显示的颜色是D。

A.黑色B.蓝色C.红色D.绿色(4)如果要以科学计数法显示15位有效数字，使用的命令是B。

long long e long g long d(5)在命令窗口新建变量a、b，如果只查看变量a的详细信息，使用的命令为B。

a a(6)如果要清除工作空间的所有变量，使用的命令为 C 。

all C.两者都可 D.两者都不可(7)在创建变量时，如果不想立即在命令窗口中输出结果，可以在命令后加上D。

A.冒号B.逗号C.空格D.分号(8)如果要重新执行以前输入的命令，可以使用B键。

A.下箭头↓B.上箭头↑C.左箭头←D.右箭头→(9)如果要查询函数inv的相关信息，并显示在命令窗口，应使用命令A。

D.三者均可(10)如果要启动Notebook文档，下列D操作是可行的。

A.在命令窗口输入notebook命令B.在命令窗口输入notebook filename命令C.在Word中启动M-book文档D.三者均可2.填空题：(1)MATLAB是MATrix和LABoratory两个单词前三个字母的组合，意为“矩阵实验室”，它的创始人是Cleve Moler和Jack Little。

(2)在MATLAB的默认设置中，关键字显示的字体为蓝色，命令、表达式、计算结果显示的字体为黑色，字符串显示的字体为紫色，注释显示的字体为绿色，错误信息显示的字体为红色。

(3)在命令窗口中，输出结果显示为各行之间添加空行的命令为format compact，各行之间不添加空行的命令为format compact。

备注：本题布置给大家时有一点小错误，现在予以更正。

Matlab上机作业部分参考答案

3. 设A为数组，B为一个行数大于3的数组，请给出（1）删除A的第4、8、12三列的命令；（2）删除B的倒数第3行的命令；（3）求符号极限（4）求的3阶导数
lim tan( mx ) 的命令集； x 0 nx x3 y arctan ln(1 e 2 x ) 的命令集； x2
1 0.8 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -0.8
0
2
4
6
8
10
12
14
某校60名学生的一次考试成绩如下: 93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55 1)计算均值、标准差、极差、偏度、峰度，画出直方图； 2)检验分布的正态性； 3)若检验符合正态分布，估计正态分布的参数并检验参数。解答： x=[93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55]; mean(x) std(x) range(x) skewness(x) kurtosis(x) hist(x) h=normplot(x) [muhat,sigmahat,muci,sigmaci]=normfit(x) [H,sig,ci]=ttest(x,80.1)

上海电院Matlab大作业

1、编写Matlab函数命令M文件，完成下列功能，要求程序语句尽可能简洁：（1）函数输入参数为正整数n，要求3<n<10，如果输入的参数不符合3<n<10，则给出提示。

（2）生成n阶随机矩阵A；（3）生成n阶矩阵B，B的每个元素是对应位置上A矩阵元素的自然对数；（4）求矩阵B的所有对角线元素之和m。

（5）返回值为矩阵B和m。

（6）要求：进行上机编程，调试完成后将程序书写在大作业报告中、并加以注释，将调试结果截图打印粘贴在大作业报告中。

解：①根据题目要求编写一个函数文件，设函数名为heer ，则函数文件名为heer.m②因为题目要求最终的返回值为B和m，故在定义函数时可设置为 function [B,m]=heer(n)③根据要求，矩阵的阶数n是有限定范围（3<n<10）的，故可采用if语句实现④在这个大前提下，矩阵A是系统产生的，且矩阵B在A的基础上进行对数运算，可以直接用赋值语句实现：A=rand(n,n)B=log(A)⑤因为题目还要求矩阵B的所有对角线上的元素和，考虑到对角线上元素为B[i,i], 即元素下标相等，则可以用for 循环语句实现：for i=1:1:nsum=sum+B(i,i);end这样，整个M文件的内容就确定下来，只需在软件中先编写好该文件，再保存，即可在命令窗口直接调用该文件。

下面为M文件的截图如下图所示，设置n=2，则系统显示错误信息如下图所示，设置n=7.，程序运行结果如下2、一个50Hz的简单电力系统如下图所示，试在Simulink中建立仿真模型研究该系统性能。

k1GTLD2LLD1系统建模要求如下：（1）发电机G 采用“Synchronous Machine pu Fundamental ”模型，变压器T 采用“Three-Phase Transformer (Two Windings)”模型，输电线路L 采用“Three-Phase Series RLC Branch ”模型，负荷LD1、LD2采用“Three-Phase Parelell RLC Load ”模型。

matlab整理

作业一：P55：1.在一个MATLAB命令中，6+7i和6+7*i有何区别？i和I有何区别？解：在MATLAB中6+7i是一个复数常量，6+7*i则是一个表达式。

i是虚数单位，而I是单位向量。

4.要产生均值为3，方差为1的500个正态分布的随机序列，写出相应的表达式。

解：y=3+sqrt(1)*randn(500)5.求下列矩阵的对角线元素、上三角矩阵、逆矩阵、行列式的值、秩、迹。

(1) A =1 -12 35 1 -4 23 0 5 211 15 0 9(2) B =0.43 43 2-8.9 4 21解：（1）A=[1,-1,2,3;5,1,-4,2;3,0,5,2;11,15,0,9]主对角元素：D=diag(A)上三角矩阵：B=triu(A)下三角矩阵：C=tril(A)逆矩阵：X=inv(A)行列式的值：E=det(A)秩：F=rank(A)逆：G=trace(A)运行结果：A =1 -12 35 1 -4 23 0 5 211 15 0 9主对角元素：D =1159上三角矩阵：B =1 -12 30 1 -4 20 0 5 20 0 0 9下三角矩阵：C =1 0 0 05 1 0 03 0 5 011 15 0 9逆矩阵：X =-0.1758 0.1641 0.2016 -0.0227 -0.1055 -0.1016 -0.0391 0.0664 -0.0508 -0.0859 0.1516 0.0023 0.3906 -0.0313 -0.1813 0.0281 行列式的值：E =1280秩：F =4逆：G =16（2）B = [0.43,43,2;-8.9,4,21]主对角元素：D=diag(B)上三角矩阵：Y=triu(B)下三角矩阵：C=tril(B)逆矩阵：X=pinv(B)行列式的值：E=det(B)秩：F=rank(B)迹：G= trace(B)运行结果:B =0.4300 43.0000 2.0000-8.9000 4.0000 21.0000主对角元素：D =0.43004.0000上三角矩阵：Y =0.4300 43.0000 2.00000 4.0000 21.0000下三角矩阵：C =0.4300 0 0-8.9000 4.0000 0逆矩阵：X =0.0022 -0.01750.0234 -0.0017-0.0035 0.0405行列式的值：E = 1.2526e+003秩：F=3迹：G =5.43006. 当A=[34,NaN,Inf,-Inf,-pi,eps,0]时，求函数all（A）、any(A)、isnan(A)、isinf(A)、isfinite(A)的值。