MATLAB课后实验答案1.docx
- 文档编号:1627650
- 上传时间:2022-10-23
- 格式:DOCX
- 页数:20
- 大小:228.68KB
MATLAB课后实验答案1.docx
《MATLAB课后实验答案1.docx》由会员分享,可在线阅读,更多相关《MATLAB课后实验答案1.docx(20页珍藏版)》请在冰豆网上搜索。
MATLAB课后实验答案1
实验一MATLAB运算基础
1、先求下列表达式的值,然后显示MATLAB工作空间的使用情
况并保存全部变量
⑵z2
1
—In(x
、1
x2),其中x2
12i
2
0.45
5
(3)z3
e0.3ae
0.3a
sin(a
0.3a
0.3)In,
2
a3.0,2.9,卅,2.9,3.0
2
t2
0t
1
⑷z4
t21
1t
2,其中t=0:
0、
5:
2、5
t22t
1
2t
3
(1)zi
2sin85°
1e2
解:
M文件:
z1=2*sin(85*pi/180”(1+exp
(2))
x=[21+2*i;-、455];
z2=1/2*log(x+sqrt(1+xA2))
a=-3、0:
0、1:
3、0;
z3=(exp(0、3、*a)-exp(-0、3、*a))、/2、*sin(a+0、3)+log((0、
3+a)、/2)
t=0:
0、5:
2、5;
z4=(t>=0&t<1)、*(t、A2)+(t>=1&t<2)、*(t、A2-1)+(t>=2&t<3)、*(t、
A2-2*t+1)
4、完成下列操作:
(1)求[100,999]之间能被21整除的数的个数。
(2)建立一个字符串向量,删除其中的大写字母。
解:
(1)结果:
m=100:
999;
n=find(mod(m,21)==0);length(n)
ans=
43
(2)、建立一个字符串向量例如:
ch二'ABC123d4e56Fg9:
则要求结果就是:
ch二'ABC123d4e56Fg9:
k=find(ch>='A'&ch<='Z');
ch(k)=[]
ch=
123d4e56g9
阵与对角阵,试通过数值计算验证
A2
RRS
S2。
实验二MATLAB矩阵分析与处理
5、下面就是一个线性方程组
1
1
1
2
3
4
%
0.95
1
1
1
0.67
X2
3
4
5
1
1
1
X3
0.52
4
5
6
(1)求方程的解。
⑵将方程右边向量元素b3改为0、53再求解拼比较b3的变化与解的相对变化。
(3)计算系数矩阵A的条件数并分析结论。
解:
M文件如下:
实验三选择结构程序设计
1、求分段函数的值。
2x
x6
x0且x3
2
yx
5x6
0x5且x2及x3
2x
x1
其他
用if语句实现,分别输出x=-5、0,-3、0,1、0,2、0,2、5,3、0,5、0时的y值。
解:
M文件如下:
2、输入一个百分制成绩,要求输出成绩等级A、B、C、D、E。
其中90分~100分为A,80
分~89分为B,79分~79分为C,60分~69分为D,60分以下为E。
要求:
(1)分别用if语句与switch语句实现。
(2)输入百分制成绩后要判断该成绩的合理性,对不合理的成绩应输出出错信息。
解:
M文件如下
3、硅谷公司员工的工资计算方法如下:
(1)工作时数超过120小时者,超过部分加发15%。
(2)工作时数低于60小时者,扣发700元。
(3)其余按每小时84元计发。
试编程按输入的工号与该号员工的工时数,计算应发工资。
解:
M文件下
2
1、根据一
6
1丄±
122232
,求n的近似值。
n分别取100、1000、10000
实验四循环结构程序设计
时,结果就是多少?
要求:
分别用循环结构与向量运算(使用sum函数)来实现。
解:
M文件如下:
运行结果如下
2、根据y1
⑴y<3时的最大
1
3
n值。
汕2^,求:
⑵与⑴的n值对应的y值。
解:
M—文件如下:
Untitled9*
3、考虑以下迭代公式:
a
Xn1
bXn
其中a、b为正的学数。
(1)编写程序求迭代的结果,迭代的终止条件为|Xn+1-Xn|W10-5,迭代初值X0=1、0,迭代次
数不超过500次。
bJb24a
(2)如果迭代过程收敛于r,那么r的准确值就是,当(a,b)的值取(1,1)、
(8,3)、(10,0、1)时,分别对迭代结果与准确值进行比较。
解:
M文件如下:
运算结果如下
5、若两个连续自然数的乘积减1就是素数,则称这两个边疆自然数就是亲密数对,该素
数就是亲密素数。
例如,2X3-1=5,由于5就是素数,所以2与3就是亲密数,5就是亲密素数。
求[2,50]区间内:
(1)亲密数对的对数。
(2)与上述亲密数对对应的所有亲密素数之与。
解:
M文件:
4、设f(x)
1
(x2)20.1
1
(x3)4
编写
0.01
个MATLAB函数文件fx、m,使得
实验五函数文件
调用f(x)时,x可用矩阵代入,得出的f(x)为同阶矩阵。
解:
函数fx、m文件:
functionf=fx(x)
%fxfx求算x矩阵下的f(x)的函数值
A=0、1+(x-2)、A2;
B=0、01+(x-3)、A4;
f=1、/A+1、/B;
k
IIs
liT~
_c
IICMu
IIcou
E
o)
報pg
u
T—
ll_
J£
pu0)
oo
CLu7F
u
e
X
—
llT~
_c
IICMu
IIcou
E
匸
報
oo
CLU7F
u
(L)麗
。
赳星A^-起(L+u)XU+寸xe+eXCXI+CXIXu(u=汕
(2)
。
赳星A^Fm(g+zU)U_OL+U"U=II(L)
(02二(oe二
ol寸二
「-l (XXJL Tux世舉<緯-xnduHx QQ -eH令< y3=g(n3);y=y1/(y2+y3) 实验八数据处理与多项式计算 2、将100个学生5门功课的成绩存入矩阵P中,进行如下处理: (1)分别求每门课的最高分、最低分及相应学生序号。 (2)分别求每门课的平均分与标准方差。 (3)5门课总分的最高分、最低分及相应学生序号。 (4)将5门课总分按从大到小顺序存入zcj中,相应学生序号存入xsxh。 提示: 上机调试时,为避免输入学生成绩的麻烦,可用取值范围在[45,95]之间的随机矩阵来表示学生成绩。 解: M文件: clc; t=45+50*rand(100,5); P=fix(t);%生成100个学生5门功课成绩 [x,l]=max(P) %x为每门课最高分行向量,1为相应学生序号 [y,k]=min(P) %y为每门课最低分行向列,k为相应学生序号mu=mean(P)%每门课的平均值行向量 sig=std(P)%每门课的标准差行向量 s=sum(P,2)%5门课总分的列向量 [X,m]=max(s)%5门课总分的最高分X与相应学生序号m[Y,n]=min(s)%5门课总分的最低分Y与相应学生序号n [zcj,xsxh]=sort(s) %zcj为5门课总分从大到小排序,相应学生序号xsxh 运行结果: 3、某气象观测得某日6: 00~18: 00之间每隔2h的室内外温度(0C)如实验表1所示。 实验表1室内外温度观测结果(°C) 时间h 6 8 10 12 14 16 18 室内温度 t1 18、 0 20、 0 22、 0 25、 0 30、 0 28、 0 24、0 室外温度 t2 15、 0 19、 0 24、 0 28、 0 34、 0 32、 0 30、0 试用三次样条插值分别求出该日室内外6: 30~18: 30之间每隔2h各点的近似温度(°C)。 解: M文件: clc; h=6: 2: 18; t1=[18、020、022、025、030、028、024、0];t2=[15、019、024、028、034、032、030、0];T1=interp1(h,t1,'spline')%室内的3次样条插值温度 T2=interp1(h,t2,'spline')%室外的3次样条插值温度 运行结果: 4、已知lgx在[1,101]区间10个整数采样点的函数值如实验表2所示。 实验表2lgx在10个采样点的函数值 x 101 1112131 41 51 61 718191 lgx 9085 01、04141、3222 1、95102、0043 1、49141、 6128 1、7076 1、78531、85131、 试求lgx的5次拟合多项式 解: M文件: p(x),并绘制出 lgx与p(x)在[1,101]区间的函数曲线。 x=1: 10: 101;y=lg10(x); P=polyfit(x,y,5)y仁polyval(P,x);plot(x,y,': o',x,y1,'-*') 5、有3个多项式P1(x)=x4+2x3+4x2+5,P2(x)=x+2,P3(x)=x2+2x+3,试进行下列操作 (1)求P(x)=P1(x)+P2(x)P3(x)。 (2)求P(x)的根。 ⑶当x取矩阵A的每一元素时,求P(x)的值。 其中: 1 1.2 1.4 A 0.75 2 3.5 0 5 2.5 (4)当以矩阵A为自变量时,求P(x)的值。 其中A的值与第(3)题相同。 解: M文件: clc;clear; p1=[1,2,4,0,5]; p2=[1,2]; p3=[1,2,3];p2=[0,0,0,p2]; p3=[0,0,p3]; p4=conv(p2,p3);%p4就是p2与p3的乘积后的多项式 np4=length(p4); np1=length(p1); p=[zeros(1,np4-np1)p1]+p4 %求p(x)=p1(x)+p2(x) x=roots(p) %求p(x)的根 A=[-11、2-1、4;0、7523、 5;052、5]; y=polyval(p,A) %x取矩阵A的每一兀素时的p(x)值 实验九数值微积分与方程数值求解 1、求函数在指定点的数值导数。 实验六高层绘图操作 3、已知 x、、 2 e 丄1n(x.1x2) 2 在-5wxw5区间绘制函数曲线。 解: M文件: clc; 、*(x<=0)+0、5*log(x+sqrt(1+x 、A2))、*(x>0); x=-5: 0、01: 5; y=(x+sqrt(p
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 课后 实验 答案
![提示](https://static.bdocx.com/images/bang_tan.gif)