数学实验及模型实验01MATLAB基础(一).pptVIP

4 矩阵元素的访问、修改 例如利用rand随机生成一个矩阵，程序为：A=rand(4,4) 结果如下： A= 0.9501 0.8913 0.8214 0.9218 0.2311 0.7621 0.4447 0.7382 0.6068 0.4565 0.6154 0.1763 0.4860 0.0185 0.7919 0.4057 A(i)：表示矩阵A的第i个元素. 在MATLAB中，所有矩阵都是以列进行存储的. 如 在上面的例子中，A(1)表示0.9501，A(2)表示0.2311，A(3)表示0.6068，A(7)表示0.4565，… A(i,j)：表示矩阵A的第i行第j列元素，括号内的逗号前的i表示行，逗号后的j表示列. 如：A(i,:)表示A的第i行的所有元素，A(:,j)表示A的第j列的所有元素. A([2,1],2:4)表示取矩阵A的第2、1行，第2至4列的元素，结果如下： ans = 0.7621 0.4447 0.7382 0.8913 0.8214 0.9218 我们可以随意更改、增加或删除矩阵以及向量中的元素： y=[1,2,4,6,9]; %表示对y赋值 y(3)=2 %表示将y的第3个元素更改为2 结果为： y=1 2 2 6 9 MATLAB也可取出向量的一个元素或一部份来做运算，如程序： 2*y(2:4)-3 %取出y的第2至第4个元素来做运算 结果为：ans=1 1 15 1.3 例题 例1 计算 . 解 程序为：sin((8+5*log2(4))/(abs(3-7))^2) 结果为：ans =0.9023 例2 计算 ，以符号的形式输出. 解 程序为：a=2*sym(sqrt(8))或者2*sym(power(8,1/2)) 结果均为：ans = 4*2^(1/2) 例3 用MATLAB计算 能得到–2吗？ 解 用三种不同方法计算： （1）程序为：a=-8;r=a^(1/3) 结果为：r =1.0000 + 1.7321i （2）程序为：sym(power(-8,1/3)) 结果为：ans = (1)+(sqrt(3))*i （3）构造 ，然后用roots函数求根 程序为：p=[1,0,0, 8];R=roots(p) 结果为： R= -2.0000 1.0000+1.7321i 由此可得用不同的MATLAB计算程序可能得到不同的结果. 例4 已知 ，计算 . 解 程序为：a=3+4*i;b=2-i;c=2*exp(i*pi/6); x=a*b/c 结果为：x =5.5801 - 0.3349i 例5 输出字符串“MATLAB6.5”. 解 程序为：s=MATLAB6.5 结果为：s=MATLAB6.5 例6 输出多项式 . 解 程序为：syms x; y=x^2+5*x+7 结果为：y=x^2+5*x+7 例7 已知多项式 ， 计算 ， . 解 程序为： f1=sym(x^2+4),f2=sym(2*x^3-3*x^2+2*x+3) f=f1+f2,g=f1*f2,expand(g) 结果为：f1=x^2+4 f2=2*x^3-3*x^2+2*x+3 f=-2*x^2+7+2*x^3+2*x g=(x^2+4)*(2*x^3-3*x^2+2*x+3) ans=2*x^5-3*x^4+10*x^3-9*x^2+8*x+12 例8 求解方程组 . 解 程序为：[x,y]=solve(x^2+2*x*y+y^2-4,x-y-3) 结果为：x=[5/2,1/2],y=[-1/2,-5/2] 所以方程组有两组解，分别为： 1.4 实验练习 实验01 MATLAB基础（一） 数学实验及模型 实验要求 1、按时上下课，不迟到、不早退； 2、爱护设备，不作非法操作；下课时摆放好凳子，关好计算机； 3、每次上课时带好笔、纸及相关资料，上课之前关闭所有与本课程无关的程序、网页等； 4、上课期间遵守实验室纪律，不得关闭学生端程序、不得修改学生端程序设置，不做与上课无关的事情（包含但不仅限于扫雷、纸牌、开心农场、看