第一章电子资料.ppt

 本章小结 简要介绍了控制系统仿真的基本概念及仿真工具MATLAB的语言特点；(1.1) 详细讲述了在系统仿真过程中可能用到的MATLAB基本功能，从数值计算功能入手，介绍了MATLAB中的数组、矩阵、绘图、函数、M文件及编程控制等基础知识。 希望通过本章的学习能够掌握MATLAB的基本用法，为后续利用MATLAB进行仿真打下基础。 * 习题：1. 查看MATLAB的目录结构，并检查计算机中安装了 哪些MATLAB工具箱。2. 使用help命令查找函数plot和plot3的帮助信息。3. 已知矩阵A为四阶魔方矩阵，矩阵B =[1 2 3 4;3 4 5 6;5 6 7 8;7 8 9 0]，求矩阵A和矩阵B的矩阵乘积和数 组乘积。 %magic(4) 4. 用符号计算验证三角等式。 % expand() * syms x y;expand(sin(x-y)) 5. 求f = sin(x)^2的导数，求 f = sin(x)^2 + cos(y)^2对y 的导数。 % diff() 6. 求积分 。%int(f,t,a,b) 7. 解方程x^3-6*x^2+11*x-5=1 。%solve()8. 绘图：f = (x^2 + y^2)^4 - (x^2 - y^2)^2， 其中–1 x 1。9. 编程求1^2+2^2+3^2+...+100^2，并计算运算时间。 %(tic toc)10. Fibonacci数组的元素满足规则：(k=1,2,…)， ；且 。试编程求该数组中 第一个大于10000的元素及序号。 * * 6. syms tao omiga t A;int(A*exp(-i*omiga*t),t,-tao/2,tao/2) 8. x=-1:0.01:1; y=x; [X,Y]=meshgrid(x,y); F=(X.^2+Y.^2).^4-(X.^2-Y.^2).^2; mesh(X,Y,F); xlabel(x);ylabel(y);title(f的图像) 9. tic a(1)=1^2; sum(1)=a(1); for i=2:100 a(i)=i^2; sum(i)=sum(i-1)+a(i); end sum_100=sum(100) toc 10. a(1)=1;a(2)=1; for i=3:50 a(i)=a(i-1)+a(i-2); if a(i)10000 element=a(i) order=i break end end * * 例如：解方程组 输入命令： [x,y,z]=solve(x+y+z=1,x-y+z=2,2*x-y-z=1)返回结果为：x = 2/3，y =-1/2，z =5/6 [x,y,z]=solve(x+y+z-1,x-y+z-2,2*x-y-z-1) ② 符号微分方程求解符号微分方程求解指令：dsolve格式：dsolve( eq1 , eq2 ,..., cond1 , cond2 ,..., v ) 说明：eq1，eq2，…为微分方程（组），可多至12个微分方程的求解；cond1，cond2，...为初始条件；v为指定自变量，默认时为t； 微分方程的各阶导数项以大写字母D表示， 如：y的一阶导数 可表示为：Dy， y的二阶导数 可表示为：D2y， y的n阶导数 可表示为：Dny。 * 例如： 求微分方程 的解  在命令窗口输入： y=dsolve( D2y+2*Dy+2*y=0 , y(0)=1 , Dy(0)=0 )运行返回结果为：y =cos(t)/exp(t) + sin(t)/exp(t)若再输入命令“ezplot(y)”，可进行符号函数绘图。 * 1.4 程序、文件和函数 Matlab不仅具有强大的数值处理和符号运算功能，而且可以像计算机高级语言一样进行程序设计。用matlab编程语言编写的程序以.m为扩展名，简称M文件。 M文件的打开方法： 单击 File-New-M-file/Script *  1.4.1 m文件编程M文件根据调用方式的不同分为命令文件和函数文件。 命令式文件