2013MATLAB数学实验答案(全)
第一次练习教学要求：熟练掌握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) ans = 0 1.2 cos1000xmxy e =，求''y syms xdiff(exp(x)*cos(902*x/1000),2) ans =(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算221100x 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项（要求将结果用向量的形式给出）。
x=[1,1]; for n=3:20x(n)=x(n-1)+x(n-2); end x x=Columns 1 through 101 123 5 8 13 21 34 55 Columns 11 through 2089 144 233 377 610 987 1597 2584 4181 67651.8 对矩阵211020411000A m ⎛⎫ ⎪- ⎪= ⎪⎪- ⎪⎝⎭，求该矩阵的逆矩阵，特征值，特征向量，行列式，计算6A ，并求矩阵,P D （D 是对角矩阵），使得1A PDP -=。
A=[-2,1,1;0,2,0;-4,1,902/1000];inv(A)ans =0.4107 0.0223 -0.4554 0 0.5000 0 1.8215 -0.4554 -0.9107 eig(A) ans =-0.5490 + 1.3764i -0.5490 - 1.3764i 2.0000 det(A) ans =4.3920[P,D]=eig(A)P = %特征向量0.3245 - 0.3078i 0.3245 + 0.3078i 0.2425 0 0 0.9701 0.8944 0.8944 0.0000 D =-0.5490 + 1.3764i 0 0 0 -0.5490 - 1.3764i 0 0 0 2.0000 P*D^6*inv(P) %A^6的值 ans =15.3661 12.1585 + 0.0000i -5.8531 0 64.0000 0 23.4124 -5.8531 + 0.0000i -1.61961.9 作出如下函数的图形（注：先用M 文件定义函数，再用fplot 进行函数作图）：1202()12(1)12x x f x x x ⎧≤≤⎪⎪=⎨⎪-<≤⎪⎩m 文件:function y=fenduan(x)if x<=1/2 y=2*xelse x<=1 y=2-2*xend end执行函数：fplot('fenduan',[0,1]); grid ontitle('第1.9题图') 得下图：第1.9题图1.10 在同一坐标系下作出下面两条空间曲线（要求两条曲线用不同的颜色表示）（1）cos sin x t y t z t =⎧⎪=⎨⎪=⎩ （2）2cos 2sin x t y t z t =⎧⎪=⎨⎪=⎩t=-10:0.01:10; x1=cos(t); y1=sin(t); z1=t;plot3(x1,y1,z1); hold onx2=cos(2*t); y2=sin(2*t); z2=t;plot3(x2,y2,z2,'m'); grid ontitle('第1.10题图') 得下图：1.11 已知422134305,203153211A B m -⎛⎫⎛⎫ ⎪ ⎪=-=-- ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭，在MATLAB 命令窗口中建立A 、B 矩阵并对其进行以下操作：(1) 计算矩阵A 的行列式的值det()A(2) 分别计算下列各式：1122,*,.*,,,,T A B A B A B AB A B A A --- 解：A=[4,-2,2;-3,0,5;1,5*902,3]; B=[1,3,4;-2,0,3;2,-1,1];det(A) ans =-117288 2*A-B ans =7 -7 0-4 0 7 0 9021 5 A*B ans =12 10 127 -14 -7 -9013 0 13537 A.*B ans =4 -6 86 0 15 2 -4510 3 A*inv(B) ans =1.0e+003 *-0.0000 0 0.0020 0.0000 0.0016 0.0001 1.0311 -0.9016 -1.4167 inv(A)*B ans =0.3463 0.5767 0.5383 0.0005 -0.0006 -0.0005 -0.1922 0.3460 0.9230 A*A ans =24 9012 4-7 22556 9 -13523 13528 22561 A' ans =4 -3 1-2 0 4510 2 5 31.12已知22()2()x f x μσ--=分别在下列条件下画出)(x f 的图形：(1)/600m σ=,μ分别为0,1,1-(在同一坐标系上作图); (2)0μ=,σ分别为1,2,4,/100m (在同一坐标系上作图). (1)x=-5:0.1:5;h=inline('1/sqrt(2*pi)/s*exp(-(x-mu).^2/(2*s^2))'); y1=h(0,902/600,x);y2=h(-1,902/600,x);y3=h(1,902/600,x); plot(x,y1,'b',x,y2,'m',x,y3,'y')grid ontitle('第1.12题')(2) z1=h(0,1,x);z2=h(0,2,x);z3=h(0,4,x); z4=h(0,902/100,x);plot(x,z1,x,z2,'y',x,z3,'m',x,z4, 'g') grid ontitle('第1.12题')z1=h(0,1,x);z2=h(0,2,x);z3=h(0,4,x); z4=h(0,902/100,x);1.13 作出24z mx y =+的函数图形。
x=-10:0.2:10;y=x;[X Y]=meshgrid(x,y);Z=902*X.^2+Y.^4; mesh(X,Y,Z); title('第1.13题')第1.13题1.14对于方程50.10200mx x --=，先画出左边的函数在合适的区间上的图形，借助于软件中的方程求根的命令求出所有的实根,找出函数的单调区间,结合高等数学的知识说明函数为什么在这些区间上是单调的,以及该方程确实只有你求出的这些实根。
最后写出你做此题的体会。
解：作图程序：（注：x 范围的选择是经过试探而得到的）x=-1.7:0.02:1.7;y=x.^5-902/200*x-0.1;plot(x,y);grid on; title('第1.14题')由图形观察，在x=-1.5,x=0,x=1.5附近各有一个实根 solve('x^5-902/200*x-0.1') ans =-1.4516870267499636199995749888894 -0.022172950190557703188753959027919 1.46277510594806546372292321961741.4573364935933870280941533926624*i + 0.0055424354962279297327028641499658 0.0055424354962279297327028641499658 - 1.4573364935933870280941533926624*i 三个实根的近似值分别为：-1.4517,-0.0222,1.4628由图形可以看出，函数在区间(,1)-∞-单调上升，在区间(1,1)-单调下降，在区间(1,)∞单调上升。
syms xdiff('x^5-902/200*x-0.1',x)结果为5*x^4-4.51solve('5*x^4-902/200') ans =-(451^(1/4)*500^(3/4))/500 (451^(1/4)*500^(3/4))/500 -(451^(1/4)*500^(3/4)*i)/500 (451^(1/4)*500^(3/4)*i)/500 vpa(ans) ans =-0.97454440927373918149075795211629 0.97454440927373918149075795211629 -0.97454440927373918149075795211629*i 0.97454440927373918149075795211629*i 得到两个实根：-0.9745与0.9745可以验证导函数在)9745.0,(--∞内为正，函数单调上升 导函数在)9745.0,9745.0(-内为负，函数单调下降 导函数在),9745.0(∞内为正，函数单调上升 根据函数的单调性，最多有3个实根。
数学实验(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、求出下列极限值。
北科大matlab数学实验报告第三次北科大matlab数学实验报告第三次Matlab 数学实验 全部答案 精华版
《数学实验》报告实验名称数学实验第三次报告学院专业班级姓名学号2013年11月学会用MATLAB绘制二维曲线,掌握gtext, legend, title, xlabel, ylabel等指令用法。
二、【实验任务】完成课后第4章第5，6，8，9题。
P79三、【实验程序】1,x=0:pi/50:4*pi;y=exp(x/3).*sin(3*x);plot(x,y ,'b*'),grid onhold onz1=exp(x/3);z2=-exp(x/3);plot(x,z1,'r-.')hold onplot(x,z2,'r-.')2,clf;x1=-pi:pi/50:pi; y1=x1.*cos(x1);x2=pi:pi/50:4*pi;y2=x2.*tan(1./x2).*sin(x2.^3);x3=1:0.1:8; y3=exp(1./x3).*sin(x3);subplot(1,3,1);plot(x1,y1,'r-'),grid ontitle('y1=x1.*cosx1')legend('x1.*cosx1')xlabel('x轴') ,ylabel('y轴')subplot(1,3,2);plot(x2,y2,'g-'),grid ontitle('y2=x2.*tan(1./x2).*sin(x2.^3)')legend('x2.*tan(1/x2).*sinx2')xlabel('x轴') ,ylabel('y轴')subplot(1,3,3);plot(x3,y3,'b-'),grid ontitle('y3=exp(1./x3).*sin(x3)')xlabel('x轴') ,ylabel('y轴')legend('exp(1./x3).*sinx3')axis tightgtext('y1=x1.*cosx1')gtext('y2=x2.*tan(1./x2).*sin(x2.^3)')gtext('y3=exp(1./x3).*sin(x3)')1，2，这次作业明显比第一次轻松了许多，在上课认真听讲后，又研读了书上的讲解和例题，起初是模仿，后来自己就能完全理解并做出题目来。
数学实验2013答案(全部)
A =
1 2 6
2 4 4
1 -5 7
>> B=[1 3 4;3 0 5;2 -1 1]
B =
1 3 4
3 0 5
2 -1 1
>> A+B,2*A-B,A*B,A.*B,A/B,A\B,A.^B
ans =
2 5 10
5 4 9
3 -6 8
ans =
1 1 8
（1） （2）
>> [X,Y]=eig(A)
X =
0 0.3015 0.3015
0 0.3015 0.3015
1.0000 -0.9045 -0.9045
Y =
2 0 0
0 1 0
0 0 1
A的特征值是2，1，1
特征向量是0和0.3015和0.3015
00.30150.3015
1.0000-0.9045-0.9045
>> a=5.3
a =
5.3000
>> b=[1 3; 2 5]
b =
1 3
2 5
>> who
Your variables are:
a b
>> whos
Name Size Bytes Class
a 1x1 8 double array
b 2x2 32 double array
Grand total is 5 elements using 40 bytes
%这是求从100米高空自由落下的皮球到最后静止经过总路程的程序
%其中s是总路程，h是弹起高度，n是弹起次数
matlab参考答案matlab第5次作业至诚2013参考答案
1 将字符串’matlab’和’is string’合并。
1-1）直接水平合并，使用strcat水平合并。
str1='matlab’str2='is string'str3=strcat(str1,str2)str3 =matlabis string1-2）分别使用strvcat,char,str2mat等指令垂直合并。
>> str4=strvcat(str1,str2)>> str5=char(str1,str2)>> str6=str2mat(str1,str2)2 运行下列指令，写出执行结果，比较c和d的区别。
2-1）a={'abcdefg'}b={'hello',1:3}c={a,b}d=[a,b]c ={1x1 cell} {1x2 cell}d ='abcdefg' 'hello' [1x3 double]区别：C为1*2的元胞，包含两个元胞元素：元胞数组a和元胞数组bD为1*3的元胞，由元胞数组a和b的各个元胞单元组成新的一个元胞数组。
2-2) 分别从c，d元胞中提取a中的子串’cde’c{1}{1}(2:4)d{1}(2:4)33.1-1）将每个学生的姓名以字符串形式存储，并将所有的字符串存入一个元胞数组Namename={'Mike';'Nancy';'Emily'}3.1-2）将每个学生的学号以字符串形式存储，并将所有的字符串存入一个字符矩阵Idid=['1001';'1002';'1003']3.1-3）将学生的所有成绩存放在一个二维数组中ScorestuCell={name,id,score}3.1-5）从步骤4）创建的元胞中提取所有学生的数学成绩。
stuCell{3}(:,2)3.2 结构数组32-1) 创建结构数组，将表格信息保存在结构数组student中。
MATLAB数学实验答案（全）
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 lim x 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 计算221100x 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参考答案matlab第7次作业至诚2013参考答案
1 t 在[0，10]范围内，)3010sin(02+=-t e y t ，计算y 的最大值，最小值，平均值和积分值。
t=0:10y=exp(-2*t).*sin(10*t+pi/6)max(y)min(y)mean(y)trapz(t,y)2 数据分析（必做）下表为2.1 >> score=[75 81 73;62 85 79;71 68 60;58 92 75;73 90 81]2.2 统计每位学生的平均分和总分>> mean(score,2) %平均分>> sum(score,2) %总分，注意cumsum 用于求累计求和，此处用sum 函数即可 >> max(score,[],2) %最高分%%方法二：mean(score')sum(score')max(score')2.3 统计每门课程的平均分和总分>> mean(score) %平均分>> sum(score) %总分>> max(score) %最高分2.4 统计每门课程的标准差>> std(score,0)2.4 按照数学成绩的升序序列对矩阵进行重新排序>> [x,ind]=sortrows(score,2)2.5 找出语文成绩最高分的学生序号>> [x,ind]=max(score(:,2))2.6 打印语文成绩最高的学生三门课程的成绩信息： 显示格式如下：×的数学成绩：××分，语文成绩：××分物理成绩：××分[x,ind]=sortrows(score,2)fprintf('%d 的数学成绩是%d,语文成绩是%d,物理成绩是%d\n',ind,x(4,1),x(4,2),x(4,3))%%方法二:[s,ind]=max(score(:,2))fprintf('%d 的数学成绩是%d,语文成绩是%d,物理成绩是%d\n', ind,score(ind,1),score(ind,2),score(ind,3))。
(完整版)MATLAB)课后实验答案[1]
实验一 MATLAB 运算基础1.先求下列表达式的值，然后显示 MATLAB 工作空间的使用情 况并保存全部变量解:4.完成下列操作：(1) 求[100,999] 之间能被21整除的数的个数。
(1) z i 2sin 85° 1 e 2Z 2 卯(x L)，其中x2 0.45 2i 5Z 3 0.3a 0.3a e e 2sin (a 0.3)3.0, 2.9,L ,2.9, 3.0 Z 4t 2t 2 t 2 1 2t 其中 t=0:0.5:2.5(2) 建立一个字符串向量，删除其中的大写字母。
解: (1)结果：m=100:999;n=fin d(mod(m,21)==0);len gth( n)ans =43(2).建立一个字符串向量例如:ch二'ABC123d4e56Fg9：则要求结果是: ch二'ABC123d4e56Fg9：k=fi nd(ch>='A'&ch<='Z'); ch(k)=[] ch =123d4e56g9 实验二MATLAB 矩阵分析与处理1.设有分块矩阵A E 3 3 °2 3 R 3 2S 2 2 ，其中E 、R 、0、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证A 2 E R RS 0 S 22 3 4 1 1 1 3 4 5 1 1 1 4 5 6 x , 0.95 x 2 0.67 x 3 0.52(1) 求方程的解。
(2) 将方程右边向量元素 b 3改为0.53再求解，并比较 b 3的变化和解的相对变化。
(3)计算系数矩阵A 的条件数并分析结论。
解：M 文件如下： 解：M 文件如下; 5.下面是一个线性方程组:实验三选择结构程序设计1. 求分段函数的值。
x2x 6 x 0且x 3y x2 5x 6 0 x 5且x 2 及x 3x2x 1 其他用if语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0 时的y值。
MATLAB数学实验第二版课后练习题含答案
MATLAB数学实验第二版课后练习题含答案课后练习题MATLAB数学实验第二版的课后练习题如下：第一章课后练习题1.编写MATLAB程序，计算并输出下列公式的结果：y = \\frac{1}{\\sqrt{2\\pi\\sigma^2}} e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}其中，x, $\\mu$, $\\sigma$ 分别由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');mu=input('请输入 mu 的值:');sigma=input('请输入 sigma 的值:');y=1/sqrt(2*pi*sigma^2) *exp(-(x-mu)^2/ (2*sigma^2));fprintf('y = %.2f\', y);2.编写MATLAB程序，求解下列方程的解：4x + y = 11\\\\x + 2y = 7答案如下：A= [4,1;1,2];B= [11;7];X=inv(A) *B;fprintf('x = %.2f, y = %.2f\', X(1), X(2));第二章课后练习题1.编写MATLAB程序，计算下列多项式的值：P(x) = x^4 - 2x^3 + 3x^2 - x + 1其中，x 由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');y=x^4-2*x^3+3*x^2-x+1;fprintf('P(%.2f) = %.2f\', x, y);2.编写MATLAB程序，绘制下列函数的图像：f(x) = \\begin{cases} x + 1, & x < 0 \\\\ x^2, & 0 \\leq x < 1 \\\\ 2x - 1, & x \\geq 1 \\end{cases}答案如下：x=-2:0.01:2;y1=x+1;y2=x.^2.* ((x>=0) & (x<1));y3=2*x-1;plot(x,y1,x,y2,x,y3);legend('y1 = x + 1','y2 = x^2','y3 = 2x - 1');总结本文提供了《MATLAB数学实验第二版》的部分课后练习题及其答案。
Matlab实验习题集答案.docx
1：用以上两种形式计算56+sin^ + e 3算术运算结果。
» 5A 6+sin(pi)+exp(3)ans =1.5645e+004» x=5A 6+sin(pi)+exp(3)»A=[1 2;1 2];» B=[l 1;2 2];» C=(A<B)&(A==B)C =0 0 03：对数d = 5 + sin7用五位定点、十五位定点以及有理数形式表示出来。
»a=5+sin(7); format short,a 5.6570» a=5+sin(7);» format long,aa =» A=[l 2 3;4 15 60;7 8 9]A =12 3 415 60 7 8 9 n i n i i i5：输入矩阵U 1 1丿。
％利用MATLAB 命令直接输入矩阵OneMatrix=ones(33J)OneMatrix =8 10 9 101 1 11 1 11 1 1‘0 0 0 0 0、c养」0 0 0 0 0,» OneMatrix=ones(2,5；l);ZeroMatrix=zeros(size(OneMatrix)) ZeroMatrix =0 0 0 0 00 0 0 0 07：生成3阶魔方矩阵。
» magic(3) ans =8 1 63 5 749 28：操作符冒号”：”的应用a）步长为1的等差数列b）步长为2的等差数列c）步长为-2的等差、递减数列» 0:1:10 ans =» 0:2:10 ans =» 10:(-2):0C）取矩阵A 的A （2）, A （3）, A （4）d）取矩阵A的第一行e）取矩阵A的第三列D把矩阵A的第一行第三列元素赋值给变吊: g）把矩阵A的第二行第列元素修改为100» A=[l 23 56;sin(3) 7 9;log(2) 6 1]» A(2,3)fl 3 100) c 二,利用A 与B 生成矩阵V 6 9 0 A ans = ans =23 » A ⑵,A ⑶,A(4)ans =0.141120008059867ans =0.693147180559945 ans =23»A(1,:)ans =1 23 56» A(:,3)ans =5691» x=A(l,3)56 » A(2,1)=100L0000 100.0000 0.6931 23.0000 7.0000 6.000056.00009.00001.00001 3 A =6 9 io ：己知矩阵 L° R _(A 0) D = (A B)曲=〔0 B 丿»A=[13;6 9];B=[1 5;0 8]; » C=[AB(:,1)];» C(l,3)=100c =1 3 1006 9 0» A=[l 3;6 9];B=[1 5;0 8];D=[A B]D =公比为10°" » A=[l 23 56;sin(3) 7 9;log(2) 6 1];»A(1,:)=[]A =0.1411 7.0000 9.00000.6931 6.0000 1.000012： 1 inspace(0, 1, 6) / %给出区间[0, 1]的6个等分点数据logspace(0, 1,6) / %给出区间的6个等比点数据,»linspace(0J,6) ans =» a=l:5 b=3:2:ll c=a>*a a = 1 2b = 3 5c = 1 4» a.*bans = 3 103 4 57 9 119 16 2521 36 55■兀 . sm(k ——),k = ±2, ±1,014：计算 2 的值。
MATLAB全部实验及答案
MATLAB全部实验及答案MATLAB全部实验及答案实验一、MATLAB基本操作实验内容及步骤4、有关向量、矩阵或数组的一些运算（1）设A=15；B=20；求C=A+B与c=a+b？（2）设A=[1 2 3;4 5 6;7 8 9]，B=[9 8 7;6 5 4;3 2 1]；求A*B 与A.*B？A*B就是线代里面的矩阵相乘A.*B是对应位置的元素相乘（3）设a=10，b=20；求i=a/b=0.5与j=a\b=2？（4）设a=[1 -2 3;4 5 -4;5 -6 7]请设计出程序，分别找出小于0的矩阵元素及其位置（单下标、全下标的形式），并将其单下标转换成全下标。
clear,clca=[1 -2 3;4 5 -4;5 -6 7];[x,y]=find(a<0);c=[];for i=1:length(x)c(i,1)=a(x(i),y(i));c(i,2)=x(i);c(i,3)=y(i);c(i,4)=(y(i)-1)*size(a,2)+x(i);endc（5）在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8]；看结果如何？如果改成运行A=[1,2;3,4]+i[5,6;7,8]，结果又如何？前面那个是虚数矩阵，后面那个出错（6）请写出完成下列计算的指令：a=[1 2 3;3 4 2;5 2 3]，求a^2=？,a.^2=？a^2= 22 16 1625 26 2326 24 28a.^2=1 4 99 16 425 4 9（7）有一段指令如下，请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X( : ) 转化为列向量（8）使用三元组方法，创建下列稀疏矩阵2 0 8 00 0 0 10 4 0 06 0 0 0方法一：clear,clcdata=[2 8 1 4 6];ir=[1 1 2 3 4 ];jc=[1 3 4 2 1];s=sparse(ir,jc,data,4,4);full(s)方法二：不用三元组法clear,clca=zeros(4,4);a(1,[1,3])=[2,8];a(2,4)=1;a(3,2)=4;a(4,1)=6;a（9）写出下列指令的运行结果>> A = [ 1 2 3 ]; B = [ 4 5 6 ];>> C = 3.^A>> D = A.^B5、已知+?=-334sin 234πt e y t 若需要计算t ∈[-1,1]，取间隔为0.01，试计算出相对应的y 值。