matlab实验答案刘卫国,MATLAB程序设计与应用(刘卫国编)课后实验答案.pptx

[摘要]实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示 MATLAB 工作空间的使用情 况并保存全部变量。

;z1 = 0.2375;-0.3819;A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*B A-B+eye(3) A*B A.*B A^3 A.^3 A/B B\A [A,B] [A([1,3],:);B^2]; ans =;;实验二 MATLAB 矩阵分析与处理;ans =;Hp = 1 Th = 4.7661e+005 Tp = 8.5175e+003;4. 已知;将方程右边向量元素b3 改为 0.53 再求解，并比较 b3 的变化和解的相对变化。

计算系数矩阵 A 的条件数并分析结论。

解： M 文件如下：;;实验三 选择结构程序设计;14;输入一个百分制成绩，要求输出成绩等级 A、B、C、D、E。

其中 90 分~100 分为 A， 80 分~89 分为 B，79 分~79 分为C，60 分~69 分为D，60 分以下为 E。

要求： 分别用 if 语句和 switch 语句实现。

输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。

解：M 文件如下;score=123 错误：输入的成绩不是百分制成绩;M 文件如下；;132;实验四 循环结构程序设计;运行结果如下： K>> %循环结构计算pi 值 y=0; n=input('n='); for i=1:n y=y+1/i/i; end pi=sqrt(6*y) n=100 pi = 3.1321 n=1000 pi = 3.1406 n=10000 pi = 3.1415 %向量方法计算 Pi 值 n=input('n='); i=1./(1:n).^2; s=sum(i); pi=sqrt(6*s) n=100 pi = 3.1321 n=1000 pi = 3.1406 n=10000 pi =3.1415;;n = 57 n = 56;23;解：M—文件;解： M 文件：;实验五 函数文件;27;输入 m1=1 输入 m2=1 输入 theta=30 X = 7.8400 3.3948 6.7896 15.6800;函数 fx.m 文件： function f= fx(x) %fx fx 求算 x 矩阵下的 f(x)的函数值 A=0.1+(x-2).^2; B=0.01+(x-3).^4; f=1./A+1./B; 命令文件： clc; x=input('输入矩阵 x='); f=fx(x);30;实验六 高层绘图操作;y2=cos(2*x); y3=y1.*y2; plot(x,y1,'b-',x,y2,'r:',x,y3,'k--');;;;;;;subplot(2,2,4); surf(x,y,z);shading interp; title('shading interp');;实验七 低层绘图操作;;3. 利用曲面对象绘制曲面 v(x,t)=10e-0.01xsin(2000πt-0.2x+π)。

解：M 文件：;;;mu=mean(x) %求这 30000 个均匀分布随机数的平均值 sig=std(x) % 求 其 标 准 差 σ1 y=length(find(x>0.5)); %找出大于 0.5 数的个数 p=y/30000 %大于 0.5 的所占百分比;;Column 7 52.000000000002444 T2 = Columns 1 through 3 34.000000000000284 42.000000000000902 52.000000000002444 Columns 4 through 6 60.000000000004512 72.000000000009408 68.000000000007503 Column 7 64.000000000005883;;p = 0 0 0 0 1 3 8 7 11 x = -1.3840 + 1.8317i -1.3840 - 1.8317i -0.1160 + 1.4400i -0.1160 - 1.4400i y = 1.0e+003 *;f3=dx(length(g)-1) %x=3 的数值倒数;51;;1289/682 (2). M 文件： function F=fun(X) x=X(1); y=X(2); z=X(3); F(1)=sin(x)+y^2+log(z)-7; F(2)=3*x+2-z^3+1; F(3)=x+y+z-5; X=fsolve('myfun',[1,1,1]',optimset('Display','off')