matlab上机-第二次.docVIP

上机目的： 1、数值计算的练习； 2、掌握循环语句使用； 3、简单m文件的编写、调试； 4、matlab优化工具箱的应用； 上机内容： 1、数值计算 加减乘除+，- ，*，/，\　左除:　2\3 = 1.5000；　^　幂: x=2; x^3; x^(-3) sqrt(x) 开平方 　 abs(x) 绝对值　abs(3-4i)　 exp(x) ex; log(x) 以e为底，x的对数 例：Z=4^2+4*2+2/2+4-2 Z=sqrt(4) 矩阵的运算： A=[4 -2 2;-3 0 5;1 5 3];B=[1 3 4;-2 0 -3;2 -1 1]; 加减运算：计算C=A+B，D=A-B，E=A+3? ?乘法运算：计算F=A*B?,G=A*3? ?除法运算：计算A/B,A\B? 2、循环语句练习 程序分析：可以看出while ... end循环的循环次数由表达式来决定，例2当n=101就停止循环，例1为自定义函数。 3、m文件练习 （1）参见PPT“ matlab-m文件与m函数”，并练习该PPT内容 （2）附加练习，编写成m文件练习，不要在命令窗口练习： 设 ，把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线。 x=linspace(0,2*pi,101); y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); plot(x,y,r) %plot画图命令,r代表红色，b蓝色 求下列联立方程的解 a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; b=[4;4;9;4]; c=a\b c = 5.2226 4.4570 1.4718 1.5994 （3）编写下面M文件，并调试 微分：f(x)=ax2+bx+c?分别以x,a作为自变量求f(x)的一次和二次微分?f syms a b c x;%定义符号 f=a*x^2+b*x+c; %表达式 diff(f,x,1) %对x一次求导 diff(f,x,2) %对x二次求导 diff(f,a,1) diff(f,a,2) （4）求N的阶乘?自定义如下函数，并保存为five2.m function?y=five2(n)? y=1;? for?m=1:n???? ?y=y*m;? end? y=y? ?在命名窗口输入five2(10) %自动调用函数 结果?y?=?3628800? 4、matlab优化工具箱的练习（练习PPt“现代设计方法课件第6节”内容） （1）线性规划问题 （2）二次规划问题 （3）无约束非线性规划问题 对比练习 6.5 约束非线性规划问题