matlab教程课后作业.docx

【例1.3-5】图示复数的和。z1=4+3*i;z2=1+2i;z12=z1+z2clf,hold onplot(0,z1,z12,-b,LineWidth,3)plot(0,z12,-r,LineWidth,3)plot(z1,z12,ob,MarkerSize,8)hold off,grid onaxis equalaxis (0,6,0,6)text(3.5,2.3,z1)text(5,4.5,z2)text(2.5,3.5,z12)xlabel(real)ylabel(image)shg z12 = 5.0000 + 5.0000i 【例1.3-6】 用MATLAB计算能得到 2 吗？a=-8;r_a=a(1/3)p=1,0,0,-a;R=roots(p)MR=abs(R(1);t=0:pi/20:2*pi;x=MR*sin(t);y=MR*cos(t);plot(x,y,b:),grid onhold onplot(R(2),.,MarkerSize,30,Color,r)plot(R(1,3),o,MarkerSize,15,Color,b)axis(-3,3,-3,3),axis squarehold off r_a = 1.0000 + 1.7321iR = -2.0000 1.0000 + 1.7321i 1.0000 - 1.7321i 【例1.3-10】画出衰减振荡曲线，的取值范围是。t=0:pi/50:4*pi;y=exp(-t/3).*sin(3*t);plot(t,y,r,LineWidth,2)axis(0,4*pi,-1,1)xlabel(t),ylabel(y) 习题31.在MATLAB中，先运行指令A=magic(3), B=1,2,1;3,4,3;5,6,7, C=reshape(1:6,3,2)生成阵列，然后根据运行结果回答以下问题： （1）计算A*B, B*A，这两个乘积相同吗？（2）计算AB, B/A，左除、右除结果相同吗？（3）计算B( : ,1,2).*C和C.*B( : , 1,2)，这两个乘积相同吗？（4）计算AA和A.A，这两个计算结果相同吗？（5）计算Aeye(3)和inv(A)，这两个计算结果相同吗？（提示：根据对计算结果的目测回答问题）A=magic(3), B=1,2,1;3,4,3;5,6,7, C=reshape(1:6,3,2) A = 8 1 6 3 5 7 4 9 2B = 1 2 1 3 4 3 5 6 7C = 1 4 2 5 3 6（1）A*B, B*A ans = 41 56 53 53 68 67 41 56 45ans = 18 20 22 48 50 52 86 98 86 结果不同，所以A*B不等于B*A（2）AB, B/A ans = 0.0333 0.1000 0.1611 0.5333 0.6000 0.7444 0.0333 0.1000 -0.1722ans = 0.0056 0.0889 0.1722 0.1389 0.2222 0.30560.2333 0.7333 0.2333 结果不同，所以AB,不等于B/A（3）B( : ,1,2).*C,C.*B( : , 1,2) ans = 1 8 6 20 15 36ans = 1 8 6 20 15 36 结果相等，所以B( : ,1,2).*C等于C.*B( : , 1,2)（4）AA,A.Aans = 1 0 0 0 1 0 0 0 1ans = 1 1 1 1 1 1 1 1 1 结果不同，所以AA等于A.A（5）Aeye(3),inv(A) ans = 0.1472 -0.1444 0.0639 -0.0611 0.0222 0.1056 -0.0194 0.1889 -0.1028ans = 0.1472 -0.1444 0.0639 -0.0611 0.0222 0.1056 -0.0194 0.1889 -0.1028 结果相同，所以Aeye(3)等于inv(A)5. 5，要求在闭区间上产生具有10个等距采样点的一维数组。试用两种不同的指令实现（提示：冒号生成法，定点生成法）(1) 冒号生成法：A=0:2*pi/9:2*pi A = 0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 4.8869 5.5851 6.2832 (2) 定点生成法： A=linspace(0,2*pi,10) A = 0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 4.8869 5.5851 6.2832 6. 由指令rng(default),A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5的元素的位置，分别求出它们的“全下标”和“单序号”。（提示：find和sub2ind）rng(default),A=rand(3,5) A = 0.8147 0.9134 0.2785 0.9649 0.9572 0.9058 0.6324 0.5469 0.1576 0.4854 0.1270 0.0975 0.9575 0.9706 0.8003 B=A=0.5 B = 1 1 0 1 1 1 1 1 0 0 0 0 1 1 1 B中所有的1为A中所有大于0.5的元素C=B(1:15) C = 1 1 0 1 1 0 0 1 1 1 0 1 1 0 1 C为A中所有大于0.5的元素的单序号7. 先运行clear, format long, rng(default), A=rand(3,3)，然后根据A写出两个矩阵：一个对角阵B。其相应元素由A的对角元素构成；另一个矩阵C，其对角元素全为0，而其余元素与对应的A阵元素相同。（提示：diag）clear, format long, rng(default), A=rand(3,3) A = 0.814723686393179 0.913375856139019 0.278498218867048 0.905791937075619 0.632359246225410 0.546881519204984 0.126986816293506 0.097540404999410 0.957506835434298 B=diag(A) B = 0.814723686393179 0.632359246225410 0.957506835434298 B为A的对角元素构成C=A-diag(B) C = 0 0.913375856139019 0.278498218867048 0.905791937075619 0 0.546881519204984 0.126986816293506 0.097540404999410 0 C为A对角元素全为0，而其余元素与对应的A阵元素相同的矩阵10，下面有一段程序，企图用来解决如下计算任务：有矩阵，当依次取10, 9, 8, 7, 6, 5, 4, 3, 2, 1时，计算矩阵“各列元素的和”，并把此求和结果存放为矩阵Sa的第k行。例如时，A阵为，此时它各列元素 的和是一个行数组，并把它保存为Sa的第3行。问题：该段程序的计算结果对吗？假如计算结果不正确，请指出错误发生的根源，并改正之。（提示：本题专为揭示sum对行数组的作用而设计。仔细观察以下程序运行后所得到的Sa正确吗？for k=10:-1:1; A=reshape(1:10*k,k,10); Sa(k,:)=sum(A); end; Sa ）for k=10:-1:1; A=reshape(1:10*k,k,10); Sa(k,:)=sum(A); end; Sa Sa = 55 55 55 55 55 55 55 55 55 55 3 7 11 15 19 23 27 31 35 39 6 15 24 33 42 51 60 69 78 87 10 26 42 58 74 90 106 122 138 154 15 40 65 90 115 140 165 190 215 240 21 57 93 129 165 201 237 273 309 345 28 77 126 175 224 273 322 371 420 469 36 100 164 228 292 356 420 484 548 612 45 126 207 288 369 450 531 612 693 774 55 155 255 355 455 555 655 755 855 955 由得出的结果可以看出，第一行的值全部为55，是因为k=1的时候，数组为一维数组，求和之后的结果为55，而按照题意，第一行的结果就是k=1的时候的数组。改正后的程序：for k=10:-1:1; A=reshape(1:10*k,k,10); Sa(k,:)=sum(A); end; SaSa(1, : )=A习题6请分别写出用for 和while 循环语句计算的程序。此外，还请写出避免循环的数值、符号计算程序。（提示：sum和“指数采用数组”配合； tic, toc可用以记录计算所花的时间。）for语句：for k=10:-1:1; A=reshape(1:10*k,k,10); Sa(k,:)=sum(A);if k=1Sa(k,:)=A;end;end;Sa Sa = 1 2 3 4 5 6 7 8 9 10 3 7 11 15 19 23 27 31 35 39 6 15 24 33 42 51 60 69 78 87 10 26 42 58 74 90 106 122 138 154 15 40 65 90 115 140 165 190 215 240 21 57 93 129 165 201 237 273 309 345 28 77 126 175 224 273 322 371 420 469 36 100 164 228 292 356 420 484 548 612 45 126 207 288 369 450 531 612 693 774 55 155 255 355 455 555 655 755 855 955 习题6（p264）：1，2请分别写出用for 和while 循环语句计算的程序。此外，还请写出避免循环的数值、符号计算程序。（提示：sum和“指数采用数组”