MATLAB)课后实验答案
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
解：(1) 结果：(2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5. 下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦ch =123d4e56g9(1) 求方程的解。
(2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解： M 文件如下：实验三 选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用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文件下实验四循环结构程序设计1. 根据2222211116123nπ=++++L，求π的近似值。
当n分别取100、1000、10000时，结果是多少？要求：分别用循环结构和向量运算（使用sum函数）来实现。
解：M文件如下：运行结果如下：2. 根据11113521yn=++++-L，求：(1) y<3时的最大n值。
(2) 与(1)的n值对应的y值。
解：M—文件如下：3. 考虑以下迭代公式：1n nax b x +=+ 其中a 、b 为正的学数。
(1) 编写程序求迭代的结果，迭代的终止条件为|x n+1-x n |≤10-5，迭代初值x 0=1.0，迭代次数不超过500次。
(2) 如果迭代过程收敛于r ，那么r 的准确值是242b b a-+，当(a,b)的值取(1,1)、(8,3)、(10,0.1)时，分别对迭代结果和准确值进行比较。
解：M 文件如下：运算结果如下；5. 若两个连续自然数的乘积减1是素数，则称这两个边疆自然数是亲密数对，该素数是亲密素数。
例如，2×3-1=5，由于5是素数，所以2和3是亲密数，5是亲密素数。
求[2,50]区间内：(1) 亲密数对的对数。
(2) 与上述亲密数对对应的所有亲密素数之和。
解：M文件：实验五 函数文件4. 设2411()(2)0.1(3)0.01f x x x =+-+-+，编写一个MATLAB 函数文件fx.m ，使得调用f(x)时，x 可用矩阵代入，得出的f(x)为同阶矩阵。
解： 函数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;5. 已知(40)(30)(20)fyf f=+(1) 当f(n)=n+10ln(n2+5)时，求y的值。
(2) 当f(n)=1×2+2×3+3×4+...+n×(n+1)时，求y的值。
(2).实验八数据处理与多项式计算2. 将100个学生5门功课的成绩存入矩阵P中，进行如下处理：(1) 分别求每门课的最高分、最低分及相应学生序号。
(2) 分别求每门课的平均分和标准方差。
(3) 5门课总分的最高分、最低分及相应学生序号。
(4) 将5门课总分按从大到小顺序存入zcj中，相应学生序号存入xsxh。
提示：上机调试时，为避免输入学生成绩的麻烦，可用取值范围在[45,95]之间的随机矩阵来表示学生成绩。
运行结果：3. 某气象观测得某日6:00~18:00之间每隔2h的室内外温度（0C）如实验表1所示。
实验表1 室内外温度观测结果（0C）时间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各点的近似温度（0C）。
解：M文件：运行结果：4. 已知lgx 在[1,101]区间10个整数采样点的函数值如实验表2所示。
实验表2 lgx 在10个采样点的函数值x 1 11 21 31 41 51 61 71 81 91 101lgx 0 1.0414 1.3222 1.4914 1.6128 1.7076 1.7853 1.8513 1.9085 1.9510 2.0043试求lgx 的5次拟合多项式p(x)，并绘制出lgx 和p(x)在[1,101]区间的函数曲线。
解：5. 有3个多项式P1(x)=x4+2x3+4x2+5，P2(x)=x+2，P3(x)=x2+2x+3，试进行下列操作：(1) 求P(x)=P 1(x)+P 2(x)P 3(x)。
(2) 求P(x)的根。
(3) 当x 取矩阵A 的每一元素时，求P(x)的值。
其中 ：1 1.2 1.40.752 3.505 2.5A --⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦(4) 当以矩阵A 为自变量时，求P(x)的值。
其中A 的值与第(3)题相同。
解：M 文件：实验九 数值微积分与方程数值求解1. 求函数在指定点的数值导数。
实验六 高层绘图操作3. 已知201ln(02x x e y x x ⎧≤⎪⎪=⎨⎪>⎪⎩ 在-5≤x ≤5区间绘制函数曲线。
解：M 文件：2. 用数值方法求定积分。
(1) 210I π=⎰的近似值。
(2) 2220ln(1)1x I dt x π+=+⎰ 解：M 文件：运行结果：3. 分别用3种不同的数值方法解线性方程组。
6525494133422139211x y z u x y z u x y z u x y u +-+=-⎧⎪-+-=⎪⎨++-=⎪⎪-+=⎩ 解：M 文件：运行结果：4. 求非齐次线性方程组的通解。
1234123412342736352249472x x x x x x x x x x x x +++=⎧⎪+++=⎨⎪+++=⎩解：M 文件：。
5. 求代数方程的数值解。
(1) 3x +sin x -e x =0在x 0=1.5附近的根。
(2) 在给定的初值x 0=1，y 0=1，z 0=1下，求方程组的数值解。
23sin ln 70321050y x y z x z x y z ⎧++-=⎪+-+=⎨⎪++-=⎩解：M 文件：运行结果：6. 求函数在指定区间的极值。
(1) 3cos log ()x x x x x f x e++=在(0,1)内的最小值。
(2) 33212112122(,)2410f x x x x x x x x =+-+在[0,0]附近的最小值点和最小值。
8. 求微分方程组的数值解，并绘制解的曲线。
123213312123'''0.51(0)0,(0)1,(0)1y y y y y y y y y y y y =⎧⎪=-⎪⎨=-⎪⎪===⎩解： 令y1=x,y2=y,y3=z; 这样方程变为:'''0.51(0)0,(0)1,(0)1x yz y xz z xy x y z =⎧⎪=-⎪⎨=-⎪⎪===⎩,自变量是t实验十 符号计算基础与符号微积分 一、1. 已知x=6,y=5，利用符号表达式求z =提示：定义符号常数x=sym(‘6’)，y=sym(‘5’)。
解：M 文件：运行结果：2. 分解因式。
(1) x 4-y 4(2) 5135解：M 文件：5. 用符号方法求下列极限或导数。
22sin tan 3013222220,1(1)2(1)(1)lim (2)lim sin 1cos(2)(3),',''(4),,,cos ln (5)(,)(2),,x x x x x x y xy x y x e e xa t x dA d A d A y y y A xdx dt dxdt t x x y f f x y x x ex x y +→→----==+--⎡⎤-==⎢⎥⎣⎦∂∂=-∂∂∂求已知分别求已知求解：M 文件：运行结果：6. 用符号方法求下列积分。
482ln 22400(1)(2)11(3)(4)(1)1x x dx x x x dx e e dx x +∞+++++⎰⎰⎰解：M文件：。
MATLAB)课后实验答案[1]
![MATLAB)课后实验答案[1]](https://360docs-1393072327.cos.ap-beijing.myqcloud.com/adaedec65727a5e9846a61d9.webp)
实验一 MATLAB 运算基础1。
先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z e =+（2) 21ln(2z x =，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩,其中t =0：0.5:2.5 解：4. 完成下列操作：（1）求［100，999]之间能被21整除的数的个数. （2) 建立一个字符串向量,删除其中的大写字母。
解：(1) 结果：(2）。
建立一个字符串向量例如：ch=’ABC123d4e56Fg9'；则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦,其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A OS +⎡⎤=⎢⎥⎣⎦。
解： M 文件如下；5。
下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦ch =123d4e56g9（1） 求方程的解。
（2) 将方程右边向量元素b 3改为0。
53再求解，并比较b 3的变化和解的相对变化。
（3) 计算系数矩阵A 的条件数并分析结论。
解： M 文件如下：实验三 选择结构程序设计1. 求分段函数的值.2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现,分别输出x=-5.0，—3.0,1.0,2。
数学实验(MATLAB)课后习题答案
数学实验练习2.1画出下列常见曲线的图形。
（其中a=1,b=2,c=3）1、立方抛物线3xy=解：x=-5:0.1:0;y=(-x).^(1/3);y=-y;x=0:0.1:5;y=[y,x.^(1/3)];x=[-5:0.1:0,0:0.1:5];plot(x,y)2、高斯曲线2x e=y-解：fplot('exp(-x.^2)',[-5,5])3、笛卡儿曲线)3(13,1333222axy y x t at y t at x =++=+=解：ezplot('x.^3+y.^3-3*x*y',[-5,5])xyx.3+y.3-3 x y = 0或t=-5:0.1:5; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)4、蔓叶线)(1,1322322xa x y t at y t at x -=+=+=解：ezplot('y.^2-x.^3/(1-x)',[-5,5])xyy.2-x.3/(1-x) = 0或t=-5:0.1:5; x=t.^2./(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)5、摆线)cos 1(),sin (t b y t t a x -=-= 解：t=0:0.1:2*pi;x=t-sin(t); y=2*(1-cos(t)); plot(x,y)6、星形线)(sin ,cos 32323233a y x t a y t a x =+== 解：t=0:0.1:2*pi; x=cos(t).^3; y=sin(t).^3;plot(x,y)或ezplot('x.^(2/3)+y.^(2/3)-1',[-1,1])xyx.2/3+y.2/3-1 = 07、螺旋线ct z t b y t a x ===,sin ,cos 解：t=0:0.1:2*pi; x=cos(t); y=2*sin(t); z=3*t; plot3(x,y,z) grid on8、阿基米德螺线θa r = 解：x =0:0.1:2*pi; r=x; polar(x,r)902701809、对数螺线θa e r = 解：x =0:0.1:2*pi; r=exp(x); polar(x,r)90270180010、双纽线))()((2cos 22222222y x a y x a r -=+=θ 解：x=0:0.1:2*pi; r=sqrt(cos(2*x)); polar(x,r)90270或ezplot('(x.^2+y.^2).^2-(x.^2-y.^2)',[-1,1]) grid onxy(x.2+y.2).2-(x.2-y.2) = 011、双纽线)2)((2sin 222222xy a y x a r =+=θ 解：x=0:0.1:2*pi; r=sqrt(sin(2*x)); polar(x,r)90270或ezplot('(x.^2+y.^2).^2-2*x*y',[-1,1]) grid onxy(x.2+y.2).2-2 x y = 012、心形线)cos 1(θ+=a r 解：x =0:0.1:2*pi; r=1+cos(x); polar(x,r)90270练习2.21、求出下列极限值。
2019年MATLAB)课后实验答案
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:: 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
解：(1) 结果：(2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5. 下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解。
(2) 将方程右边向量元素b 3改为再求解，并比较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解： M 文件如下：123d4e56g9实验三 选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=,,,,,,时的y 值。
MATLAB)课后实验答案
实验一 MATLAB 运算基础1、 先求下列表达式得值，然后显示MATLAB 工作空间得使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =+，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0、5:2、5 解：4、 完成下列操作：(1) 求[100,999]之间能被21整除得数得个数。
(2) 建立一个字符串向量，删除其中得大写字母。
解：(1) 结果：(2)、 建立一个字符串向量 例如：ch='ABC123d4e56Fg9';则要求结果就是：实验二 MATLAB 矩阵分析与处理1、 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵与对角阵，试通过数值计算验证22E R RS A OS +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5、 下面就是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程得解。
(2) 将方程右边向量元素b 3改为0、53再求解，并比较b 3得变化与解得相对变化。
(3) 计算系数矩阵A 得条件数并分析结论。
解： M 文件如下：实验三 选择结构程序设计1、 求分段函数得值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5、0,-3、0,1、0,2、0,2、5,3、0,5、0时得y 值。
MATLAB程序设计与应用课后实验答案
实验一 MATLAB 运算基础欧阳家百(2021.03.07)0<r<lz 4 =<r 2 -11 <r < 2解./Ul •M 文件：z 1 =2*sin(85*pi/l 80)/( 1 +exp(2))x=[2 1+2衍;-.45 习;z2= l/2*log(x+sqrt(l +x A 2))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.A 2)+(t>=1 &t<2)/(t.A 2-1 )+(t>=2&t<3) ?(t.A 2-2*t+1)运算结果:z1=2*sin(85*pi/180)/(1 +exp(2)) x=[2 1+2*i ；-.45 5]; z2=1/2*log(x+sqrt(1 +x A 2))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.A 2)+(t>=1 &t<2).*(t.A 2-1 )+(t>=2&t<3) .*(t.A 2-2*t+1) Z1 =1・先求下列表达式的值, 用情况并保存全部变量。
然后显示MATLAB 工作空间的使 2 sin 85° Z\ = ；-(1) 2 Zi= — ln(x +)⑵-2\其中严Y 5 = 12-0.45二映+ 0・3"¥1 + 2/ 5 _G = _3Q-29・・・，2・9,3・0t 2-2t + l 2<t<3其中 匸0:0.5250.2375 z2 =0.7114 ・ 0.0253i 0.8968 + 0.3658i0.6474 0.64700.63510.6119Columns 57 through 600.2139 + 0.9343i 1.1541 ・ 0.0044iz3 =Columns 1 through 4 0.7388+ 3.1416i 0.7696 3.1416i0.7871 +3.1416i0.7913 + 3.14161Columns 5 through 8 0.7822 + 3.14161 0.7602 3.1416i0.7254+ 3.1416i0.6784+ 3.1416iColumns 9 through 12 0.6196 + 3.14161 0.5496 3.1416i0.4688+ 3.1416i0.3780 + 3.14161Columns 13 through 16 0.2775+ 3.1416i 0.1680 3.1416i0.0497+ 3.1416i-0.0771 +3.1416iColumns 17 through 20 -0.2124 + 3.14161 -0.3566 3.1416i -0.51043.14161 -0.6752+ 3.1416iColumns 21 through 24 -0.8536 + 3.14161 -1.0497 3.1416i -1.27013.14161 -1.5271 +3.1416iColumns 25 through 28 -1.8436 + 3.14161 -2.2727 3.1416i -2.98373.1416i -37.0245Columns 29 through 32 -3.0017-2.3085Columns 33 through 36 -1.3575-1.1531Columns 37 through 40 -0.6567-0.5151Columns 41 through 44 -0.1374-0.0255Columns 45 through 48 0.26630.3478Columns 49 through 52 0.53790.5815Columns 53 through 56-1.8971-1.5978-0.9723-0.8083-0.3819-0.25610.07920.17660.42060.48410.6145 0.63660.5777 0.5327 0.4774 0.4126 Column 610.3388Z4 =0 0.2500 0 1.2500 1.0000 2.25002.已知：求下列表达式的值：(1)A+6*B和A-B+I (其中I为单位矩阵)(2)A*B 和A.*B(3)A A3 和A.A3(4)A/B 及B\A(5)[A,B]和[A([1,3],:);B7] 解.J UT •M文件：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*BA. *BA"3A. A3A/BB\A[A,B][A([l,3]z：)；B^2]运算结果：A=[12 34 -4；34 7 87；3 65 7]；B=[1 3-1；2 0 3；3 -2 7)；A+6/BA-B+eye(3)A*BA/BA A3A.A3A/BB\A[A,B][A([1,3],：)；BA2]ans =18 52 -1046 7 10521 53 49ans =12 31 ■332 8 840 67 1 ans =68 44 62309 ■72596 154 -5 241 ans =12 102 468 0 2619 - 130 49 ans =37226 247370 78688233824 48604 i 149188 600766454142 118820ans =1728 39304 -6439304 343 65850327 274625 343ans =16.4000 -13.6000 7.600035.8000 -76.2000 50.200067.0000 ■134.000068.0000ans =109.4000 ・131.2000 322.8000-53.0000 85.0000 -171.0000-61.6000 89.8000 -186.2000ans =12 34 ■4 1 3 -134 7 87 2 0 33 65 7 3 - -2 7ans =12 34 ■43 65 74 5 111 0 1920 -5 403•设有矩阵A和B(1)求它们的乘积C。
MATLAB数学实验答案(全)
第一次练习教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。
补充命令vpa(x,n) 显示x 的n 位有效数字，教材102页fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin limx mx mx x →-与3sin limx mx mxx →∞- syms xlimit((902*x-sin(902*x))/x^3) ans =366935404/3limit((902*x-sin(902*x))/x^3,inf)//inf 的意思 ans = 0 1.2 cos1000xmxy e =，求''y syms xdiff(exp(x)*cos(902*x/1000),2)//diff 及其后的2的意思 ans =(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算2211x y edxdy +⎰⎰dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)//双重积分 ans = 2.13941.4 计算4224x dx m x +⎰ syms xint(x^4/(902^2+4*x^2))//不定积分 ans =(91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求//高阶导数syms xdiff(exp(x)*cos(902*x),10) ans =-356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x)1.6 0x =的泰勒展式(最高次幂为4).syms xtaylor(sqrt(902/1000+x),5,x)//泰勒展式 ans =-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3)/91733851-(125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/902 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4,)n n n x x x n --=+=用循环语句编程给出该数列的前20项（要求将结果用向量的形式给出）。
MATLAB课后实验答案
第四章作业4.1、a=input('请输入一个4位数:');while (a<1000|a>9999)a=input('输入错误，请重新输入一个4位数:'); endb=fix(a/1000);c=rem(fix(a/100),10);d=rem(fix(a/10),10);e=rem(a,10);b=b+7;c=c+7;d=d+7;e=e+7;b=rem(b,10);c=rem(c,10);d=rem(c,10);e=rem(e,10);g=b;b=d;d=g;g=c;c=e;e=g;a=1000*d+100*e+10*b+c;disp(['加密后：',num2str(a)])4.2、a=input('请输入a: ');b=input('请输入b: ');c=input('请输入c: ');x=0.5:1:5.5;x1=(x>=0.5&x<1.5);x2=(x>=1.5&x<3.5);x3=(x>=3.5&x<=5.5);y1=a.*(x.^2)+b.*x+c;y2=a*(sin(b)^c)+x;y3=log(abs(b+c./x));y=y1.*x1+y1.*x2+y3.*x3;disp(y)x=fix(rand(1,20)*89)+10;x1=fix(sum(x)/20);disp(['平均数是：',num2str(x1)])m=(rem(x,2)==0&x<x1);n=find(m);disp(['小于平均数的数是：',num2str(x(n))]);4.4、A=input('请输入20个数的一个行向量：');m=A;a=m;b=m;for m=Aif a>=ma=m;elseif b<=mb=m;endenddisp(['最小数是：',num2str(a)])disp(['最大数是：',num2str(b)])4.5、s=0;a=0;for b=1:64c=2^a;a=a+1;s=s+c;enddisp(['2的0次方到63次方的和是：',num2str(s)]) 4.6.1、sum1=0;for n=1:100x=(-1)^(n+1)*(1/n);sum1=sum1+x;disp(['当n取100时：sum=',num2str(sum1)])sum2=0;for n=1:1000x=(-1)^(n+1)*(1/n);sum2=sum2+x;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=0;for n=1:10000x=(-1)^(n+1)*(1/n);sum3=sum3+x;enddisp(['当n取10000时：sum=',num2str(sum3)])4.8、clear alla=input('请输入一个矩阵：');b=input('请再输入一个矩阵：(注意：两矩阵要可以相乘)'); [f1,f2]=juzhenji(a,b);5.1x=-2*pi:pi/1000:2*pi;y=100./(1+x.^2);plot(x,y)x=-2*pi:pi/1000:2*pi;y=(1/2*pi)*exp(-0.5*x.^2);plot(x,y)t=-2*pi:pi/1000:2*pi;x=sin(t);y=cos(t);plot(x,y)t=-2*pi:pi/1000:2*pi;x=t.^2;y=5*t.^3;plot(x,y,'r:')实验31、clearm=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];for x=mif x<0&x~=-3y1=x^2+x-6;disp(['y= ',num2str(y1)])elseif x>=0&x<5&x~=2&x~=3y2=x^2-5*x+6;disp(['y= ',num2str(y2)])elsey3=x^2-x-1;disp(['y= ',num2str(y3)])endend2、clear allx=input('请输入成绩：');while (x>100|x<0)x=input('输入错误，请重新输入成绩：'); endswitch fix(x/10)case {9,10}disp('A')case {8}disp('B')case {7}disp('C')case {6}disp('D')case {0,1,2,3,4,5}disp('E')end3、clear allx=input('请输入工号：');y=input('请输入工作时长：');if y>120a=84*120+(y-120)*84*0.15;disp(['您本月工资是：',num2str(a)])elseif y<60b=84*y-700;disp(['您本月工资是：',num2str(b)])elsec=84*y;disp(['您本月工资是：',num2str(c)])end4、clear alla=fix(rand(1)*89)+10;disp(['a=',num2str(a)])b=fix(rand(1)*89)+10;disp(['b=',num2str(b)])x=input('请输入一个四则运算符号','s');if abs(x)==43c=a+b;elseif abs(x)==45c=a-b;elseif abs(x)==42c=a*b;elseif abs(x)==47c=a/b;enddisp(['a于b的计算结果是：',num2str(c)])5、clear alla=fix(rand(5,6)*89)+10;n=input('请输入一个数（输出结果是一个5*6的矩阵的该行元素）：'); if n>5|n<0b=a(5:5,:);disp(['输入错误，程序输出最后一行的元素：',num2str(b)])elsec=a(n:n,:);disp(['输出结果是：',num2str(c)])end实验41、n=input(“请输入n值：”)；s1=0;for k=1:n;x=1/(k^2);s1=s1+x;pi=sqrt(6*s1);2、clearn=1;y=0;while (y<3)x=1/(2*n-1);n=1+n;y=y+x;enddisp(['y<3时n的最大值是：',num2str(n-2)])disp(['相应的的y值是：',num2str(y-x)])5、clearp=0;l=0;H=[];for m=2:49a=m;b=m+1;c=a*b-1;for k=1:cif rem(c,k)==0l=l+1;endendif l==2disp(['亲密数对是','(',num2str(a),' ,',num2str(b),')'])p=p+1;H=[H,c];endl=0;endsum=sum(H);disp(['亲密数对的个数是：',num2str(p)])disp(['亲密数对的和是：',num2str(sum)])实验61、clear allx=linspace(0,2*pi,101);y=(0.5+3*sin(x)/(1+x.^2)).*cos(x);plot(x,y);title('y=[0.5+3sin(x)/(1+x^{2})]cos(x)的图像如下：') 2/1 clear allx=linspace(-2*pi,2*pi,500);y1=x.^2;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'b-',x,y2,'r--',x,y3,'k:');2//3 clear allx=linspace(-2*pi,2*pi,500);y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(3,3,1);bar(x,y1); %条形图subplot(3,3,2);stairs(x,y1); %阶梯图subplot(3,3,3);stem(x,y1); %杆图subplot(3,3,4);bar(x,y2);subplot(3,3,5);stairs(x,y2);subplot(3,3,6);stem(x,y2);subplot(3,3,7);bar(x,y3);subplot(3,3,8);stairs(x,y3);subplot(3,3,9);stem(x,y3);3、clear allt=linspace(-5,5,200);H=[];for x=tif x<=0y1=(x+sqrt(pi))/exp(2);elseif x>0y1=0.5*log(x+sqrt(1+x.^2));endH=[H,y1];endplot(t,H,'r:');4、clear alla=input('请输入一个数a：');b=input('请输入一个数b：');n=input('请输入一个数n：');theta=linspace(0,2*pi,500);rho=a*sin(b+n.*theta);polar(theta,rho);5、clear allx1=linspace(-5,5,21);y1=linspace(0,10,31);[x,y]=meshgrid(x1,y1);z=cos(x).*cos(y).*exp(-0.25.*sqrt(x.^2+y.^2)) subplot(2,1,1);surf(x,y,z);subplot(2,1,2);contour3(x,y,z,30);6、clear alls=linspace(0,0.5*pi,10); t=linspace(0,1.5*pi,30); [S,T]=meshgrid(s,t);x=cos(S).*cos(T);y=cos(S).*sin(T);z=sin(S);colormap(prism);figure(1)surf(x);figure(2)surf(y);shading flat; figure(3)surf(z);shading interp。
MATLAB程序设计与应用(第二版)课后实验答案
Matlab课后实验题答案实验一 MATLAB运算基础1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。
(1）0 122sin851ze =+(2）21ln(2z x=，其中2120.455ix+⎡⎤=⎢⎥-⎣⎦(3）0.30.330.3sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22a ae e az a a--+=++=--（4)2242011122123t tz t tt t t⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t=0:0.5：2.52. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值:(1) A+6＊B 和A —B+I(其中I 为单位矩阵） （2） A*B 和A.＊B （3） A^3和A.^3 (4） A/B 及B\A（5） ［A,B ］和［A （[1,3］，:);B^2] 解：3. 设有矩阵A 和B123453166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦（1) 求它们的乘积C 。
(2） 将矩阵C 的右下角3×2子矩阵赋给D 。
(3) 查看MATLAB 工作空间的使用情况。
4. 完成下列操作：（1） 求[100，999]之间能被21整除的数的个数。
（2) 建立一个字符串向量，删除其中的大写字母。
解:(1） 结果:（2）.建立一个字符串向量 例如： ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1。
设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A OS +⎡⎤=⎢⎥⎣⎦。
(完整版)MATLAB)课后实验答案[1]
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
解：(1) 结果：(2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5. 下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解。
(2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解： M 文件如下： 123d4e56g9实验三 选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
数学实验(MATLAB)课后习题答案
数学实验练习2.1画出下列常见曲线的图形。
（其中a=1,b=2,c=3）1、立方抛物线3xy=解：x=-5:0.1:0;y=(-x).^(1/3);y=-y;x=0:0.1:5;y=[y,x.^(1/3)];x=[-5:0.1:0,0:0.1:5];plot(x,y)2、高斯曲线2x e=y-解：fplot('exp(-x.^2)',[-5,5])3、笛卡儿曲线)3(13,1333222axy y x t at y t at x =++=+=解：ezplot('x.^3+y.^3-3*x*y',[-5,5])xyx.3+y.3-3 x y = 0或t=-5:0.1:5; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)4、蔓叶线)(1,1322322xa x y t at y t at x -=+=+=解：ezplot('y.^2-x.^3/(1-x)',[-5,5])xyy.2-x.3/(1-x) = 0或t=-5:0.1:5; x=t.^2./(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)5、摆线)cos 1(),sin (t b y t t a x -=-= 解：t=0:0.1:2*pi;x=t-sin(t); y=2*(1-cos(t)); plot(x,y)6、星形线)(sin ,cos 32323233a y x t a y t a x =+== 解：t=0:0.1:2*pi; x=cos(t).^3; y=sin(t).^3;plot(x,y)或ezplot('x.^(2/3)+y.^(2/3)-1',[-1,1])xyx.2/3+y.2/3-1 = 07、螺旋线ct z t b y t a x ===,sin ,cos 解：t=0:0.1:2*pi; x=cos(t); y=2*sin(t); z=3*t; plot3(x,y,z) grid on8、阿基米德螺线θa r = 解：x =0:0.1:2*pi; r=x; polar(x,r)902701809、对数螺线θa e r = 解：x =0:0.1:2*pi; r=exp(x); polar(x,r)90270180010、双纽线))()((2cos 22222222y x a y x a r -=+=θ 解：x=0:0.1:2*pi; r=sqrt(cos(2*x)); polar(x,r)90270或ezplot('(x.^2+y.^2).^2-(x.^2-y.^2)',[-1,1]) grid onxy(x.2+y.2).2-(x.2-y.2) = 011、双纽线)2)((2sin 222222xy a y x a r =+=θ 解：x=0:0.1:2*pi; r=sqrt(sin(2*x)); polar(x,r)90270或ezplot('(x.^2+y.^2).^2-2*x*y',[-1,1]) grid onxy(x.2+y.2).2-2 x y = 012、心形线)cos 1(θ+=a r 解：x =0:0.1:2*pi; r=1+cos(x); polar(x,r)90270练习2.21、求出下列极限值。