Matlab作业
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)的值。
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实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案
“MATLAB”练习题要求:抄题、写出操作命令、运行结果，并根据要求，贴上运行图。
1、求230x e x -=的所有根.（先画图后求解）（要求贴图）>> solve('exp(x)—3*x^2',0)ans =—2*lambertw （—1/6＊3^（1/2））-2＊lambertw(—1,—1/6＊3^(1/2）)—2*lambertw （1/6*3^(1/2))3、求解下列各题：1）30sin lim x x x x ->->> sym x ；〉> limit((x-sin （x))/x^3)ans =1/62） (10)cos ,x y e x y =求>> sym x;>> diff （exp(x ）＊cos(x)，10）ans =（-32）＊exp(x)*sin （x)3）21/20(17x e dx ⎰精确到位有效数字）〉〉 sym x;〉〉 vpa((int(exp(x^2），x,0,1/2)），17)ans =0.544987104183622224）42254x dx x+⎰〉> sym x ；>〉 int （x^4/(25+x^2）,x)ans =125*atan （x/5) - 25＊x + x^3/35)求由参数方程arctan x y t⎧⎪=⎨=⎪⎩dy dx 与二阶导数22d y dx 。
〉> sym t;>> x=log(sqrt （1+t^2)）;y=atan(t)；〉> diff （y ，t ）/diff （x ，t)ans =1/t6)设函数y =f (x ）由方程xy +e y = e 所确定,求y ′(x ).>> syms x y ；f=x ＊y+exp(y ）—exp （1);〉> -diff(f,x ）/diff （f,y)ans =-y/（x + exp （y))7)0sin 2x e xdx +∞-⎰>〉 syms x ；>〉 y=exp(-x)＊sin(2*x ）；〉> int(y ，0,inf ）ans =2/58） 08x =展开（最高次幂为）〉> syms xf=sqrt(1+x);taylor(f,0，9)ans =— (429＊x^8）/32768 + (33*x^7）/2048 — (21＊x^6)/1024 +（7*x^5）/256 - (5＊x^4)/128 + x^3/16 - x^2/8 + x/2 + 19) 1sin (3)(2)x y e y =求〉> syms x y ；>〉 y=exp(sin （1/x)）;>〉 dy=subs （diff(y,3）,x ，2）dy =—0.582610)求变上限函数2x x ⎰对变量x 的导数.>> syms a t ；>〉 diff （int(sqrt(a+t),t,x ，x^2）)Warning: Explicit integral could not be found 。
MATLB实验作业
实验一MATLAB运算基础1、用逻辑表达式求下列分段函数的值。
t=0:0.5:2.5y=t.^2.*((t>=0)&(t<1))+(t.^2-1).*(( t>=1)&(t<2))+(t.^2-2*t+1).*((t>=2 )&(t<3))2、求[100，999]之间能被21整除的数的个数。
p=rem([100:999],21)==0;sum(p)3、建立一个字符串向量，删除其中的大写字母。
ch='KdDdfKaWdsfCI',k=find(ch>=' A'&ch<='Z'),ch(k)=[]4、输入矩阵，并找出A中大于或等于5的元素。
A=[1 2 3;4 5 6;7 8 9],[m,n]=find(A>=5),for j=1:length(m)x(j)=A(m(j),n(j))xend5、求矩阵的行列式值、逆和特征根。
a11=input('a11='),a12=input('a12 ='),a21=input('a21='),a22=input('a22 ='),A=[a11,a12;a21,a22],DA=det(A),IA=inv(A),EA=eig(A) 6、不采用循环的形式求出和式的数值解。
sum(2.^[0:63])实验二1、1行100列的Fibonacc数组a，a(1)=a(2)=1，a(i)=a(i-1)+a(i-2)，用for循环指令来寻求该数组中第一个大于10000的元素，并指出其位置i。
n=100;a=ones(1,n);for i=3:na(i)=a(i-1)+a(i-2);if a(i)>10000a(i),break;end;end,i2、编写M脚本文件，定义下列分段函数，并分别求出当（）、（）和（）时的函数值。
MATLAB与控制系统仿真大作业
>> G1=tf(num,den);
>> G1=zpk(G1)
G1 =
5 (s+0.6) ----------------(s+3) (s+2) (s+1)
>> num=[2 1]; den=[1 2.9 1]; G2=tf(num,den); >> G2=zpk(G2)
G2 =
2 (s+0.5)
>> xlabel('x') ylabel('exp(-2.*x)')
三、simulink 建模（每题 10 分，共 20 分）
1、已知单位负反馈的开环传递函数为 G(s) =
2 s2 + 4s
，试利用
simulink
建立系在单位阶跃输入作用下的模型。
要求答案包括：（1）simulink 建模结构图；
（2）在同一个坐标中的阶跃信号和响应曲线图。
2、已知系统的开环传递函数为
G1(s)
=
2s2 + 5s + 6 s2 + 2s + 3
、G2 (s)
=
s2
s+6 + 7s +1
，H
(s)
=
5(s + 2) s +10
求：建立 Simulink 仿真模型，并求出其系统在单位阶跃响应；
要求答案包括：（1）simulink 建模结构图；
x=
-2.9709
0.5491
3.6000
0.0509
2. 已知下列矩阵
2 3 1
−1 3 5
A
=
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 作业答案1.绘出函数f(x)=sin x x ，在[0,4]上的图形解：在M 文件输入：x=0:pi/100:4;y=x.*sin(x);plot(y)运行2. 求3x +2x +5 = 0的根解：在命令窗口输入：>> solve('x^3+2*x+5=0')ans =((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3) - 2/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3))1/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) - ((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)/2 -(3^(1/2)*i*(2/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) + ((108^(1/2)*707^(1/2))/108 -5/2)^(1/3)))/21/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) - ((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)/2 +(3^(1/2)*i*(2/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) + ((108^(1/2)*707^(1/2))/108 -5/2)^(1/3)))/23.321436min x x x z ++=120..321=++x x x t s301≥x5002≤≤x203≥x解：运用单纯形法计算此题，首先把约束条件化成标准形式：,,,,,205030120654321635241321≥=-=+=-=++x x x x x x x x x x x x x x x（1）在M 文件输入SimpleMthd 函数：function [x,minf] = SimpleMthd(A,c,b,baseVector)sz = size(A);nVia = sz(2);n = sz(1);xx = 1:nVia;nobase = zeros(1,1);m = 1;for i=1:nViaif (isempty(find(baseVector == xx(i),1)))nobase(m) = i;m = m + 1;else;endendbCon = 1;M = 0;while bConnB = A(:,nobase);ncb = c(nobase);B = A(:,baseVector);cb = c(baseVector);xb = inv(B)*b;f = cb*xb;w = cb*inv(B);for i=1:length(nobase)sigma(i) = w*nB(:,i)-ncb(i);end[maxs,ind] = max(sigma);if maxs <= 0minf = cb*xb;vr = find(c~=0 ,1,'last');for l=1:vrele = find(baseVector == l,1);if (isempty(ele))x(l) = 0;elsex(l)=xb(ele);endendbCon = 0;elsey = inv(B)*A(:,nobase(ind));if y <= 0disp('不存在最优解!');x = NaN;minf = NaN;return;elseminb = inf;chagB = 0;for j=1:length(y)if y(j)>0bz = xb(j)/y(j);if bz<minbminb = bz;chagB = j;endendendtmp = baseVector(chagB);baseVector(chagB) = nobase(ind);nobase(ind) = tmp;endendM = M + 1;if (M == 1000000)disp('找不到最优解!');x = NaN;minf = NaN;return;endend（2）在命令窗口输入：clear allA=[1 1 1 0 0 0;1 0 0 -1 0 0;0 1 0 0 1 0;0 0 1 0 0 -1];c=[6 3 4 0 0 0];b=[120;30;50;20];[xm,mf]=SimpleMthd(A,c,b,[3 4 5 6])xm =0 50 70mf =4304.计算下面函数在区间(0，1)内的最小值。
matlab作业
matlab作业第一次作业[实验目的]1。
了解并熟悉Matlab，并在此基础上学习基本操作2.掌握Matlab 的基本操作和常用命令3。
了解Matlab常用的函数、运算符和表达式4.掌握Matlab的工作方法和M文件的相关知识5.在Matlab中学习矩阵和数组的运算1，[实验任务]1。
P11:例1_4写一个函数，求方程AX 2+Bx+C = 0的解2.P16: 1。
计算表达式e 12+23 3log25/tan 21.3。
写m命令文件找到？k +？1/k的值是k？1k？15010 4.编写函数文件并计算？k！并找出k=20时表达式的值k？1n3。
例2_2使用函数生成矩阵A =100020003，B = 00200000000104。
例2_3 (1)生成三阶魔方矩阵A；(2)生成一个四阶单位矩阵B 5。
示例2_10数组操作示例6。
例2_11计算sin(k？/2)(k=？2，？1，0)值7。
P27:2。
矩阵A=456，B=556，计算A*B，A.*B，并比较差值789123468323。
已知矩阵A=5，212，B=，计算A>B，A==B，A2)12，【实验步骤和结果】例1 _ 4函数y =洁(A，B，c)如果(绝对值(A)>洁(1，2，2) 9有复数根: >洁(A，B，c)(exp(1)12)+(23 3)* log2(5)/tan(21) ans =21.4426 e+0053。
函数s = J(m)s = 0；k = 1；而k v=[1 2 3] A=diag(v，0) B=diag(v，1)v =1 2 33A =1 0 02 0 03 B =100 000 200 000 300 000A =魔法(3)A =8 1 6 3 5 7 4 9 2> > B =眼睛(4)B =1 0 0 0 0 0 0 0 0 0 0 1例2-11。
a = 1:5 a =1 2 3 4 5> > b = 3:2:11 b =43 5 7 9 11> > a.ans =1 4 9 16 25> > a . * b94 5 6；7 8 9]A =1 2 3 4 5 6 7 8 9> > B =[4 6 8；5 5 6；3 2 2]B = 4 6 8 5 5 6 3 2> > A * Bans =23 22 26 59 61 74 95 100 122> > A . * Bans =4 12 245。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《Matlab/Simulink电力系统建模与仿真》
上机实验报告
班级：15电气工程及其自动化二班
学号：154139240096
姓名：汤嘉旺
实验一：Powergui在简单电力系统潮流计算中的应用实例
1.1实验内容与要求
完成2机5节点电力系统的潮流计算，以2机5节点电力系统为模型进行Matlab/Simulink电力系统建模与仿真。
并完成电力系统元件的模型选择、模型参数的计算及设置、计算结果及比较。
2机5节点电力系统图
1.2 Simulink建模原理图和主要模块参数设置
1、发电机模型
在该系统中的两台发电机均选用p.u.标准同步电机模块
“Synchronous Machine pu Standard”,该模块使用标幺值参数，以转子dq 轴建立的坐标系为参数，定子绕组为星形连接。
2、变压器模块
系统中的两台变压器均选用三相双绕组变压器模块
“Three-phase Transformer(Two Windings)”,采用Y-Y连接方式。
3、线路模块
系统中带有地导纳的线路选用三相“II”形等值模块
“Three Phase PI Section Line”,没有对地导纳的线路选用三相串联RLC支路模块“Three Phase Series RLC Branch”。
4、负荷模块
在SimPowerSystems库中，利用R、L、C的串联或并联组合，提供了两个静态三相负荷模块。
这两种模块是用恒阻抗支路模拟负荷，在仿真时，在给定的频率下负荷阻抗为常数。
5、母线模块
选择带有测量元件的母线模型，及三相电压电流测量元
件“Three-Phase V-I Measurement”来模拟系统中的母线，同时方便测量流过线路的潮流，在线路元件的两端也设置了该元件。
2机5节点电力系统潮流计算仿真模型图
电机运行参数图
线路运行参数图
变压器运行参数图
负载参数图
1.3 Simulink仿真和结果分析
潮流计算结果图
在完成所有的设置工作后，在“潮流计算和电机初始化”窗口中单击“更行潮流*（Update Load Flow）”,就能得到潮流计算的结果。
将利用Powergui得到的各节点电压向量与利用MATPOWER程序得到的节点电压向量进行比较，结果见下表，若以MATPOWER程序计算结果为基准，由Powergui得到的电压幅值的最大差值仅为1.067%，角度的最大差值为-3.607%，可见其结果是正确的。
稳态电压电流计算结果图
潮流计算结果比较表
本次实验利用了Powergui计算简单电力系统潮流的方法，正确地计算和设置各模型的参数是得到准确仿真结果的前提和保障。
掌握该方法不但能够更深入地理解电力系统潮流计算的过程和Powergui的功能，而且是利用Powergui进行电力系统暂态过程及稳定性的坚实基础。
实验二：电力系统暂态稳定性Simulink仿真
1.1实验内容与要求
1.网络参数及运行参数计算
2.系统转移电抗和功率特性计算
3.系统极限切除角计算
4.发电机摇摆曲线δ-t计算
5.Simulink模型及仿真结果
单机无穷大系统图
序网及短路时的等值电路图
a)负序网络b)零序网络c)短路时的等值电路
1.2 Simulink建模原理图和主要模块参数设置
电力系统暂态稳定性Simulink仿真模型图
电机参数
断
路器参数
故障器参数
负载参数
故障0.1S后切除线路，发电机转速变化曲线
仿真结果分析：
改变故障模块中的短路类型，就可以仿真系统在发生各种短路时的暂态稳定性，同样改变系统中的元件参数就可以研究各种参数对系统的暂态稳定性的影响。
总结：本次实验利用了Simulink仿真电力系统暂态稳定性，了解当发电机突然短路时，定子各绕组电流将包含基频分量、倍频分量，具体如上图所示。