MATLAB上机作业
数值分析上机作业(MATLAB)
将系数矩阵 A 分解为：A=L+U+D
Ax=b
⇔ (D + L +U)x = b ⇔ Dx = −(L + U )x + b ⇔ x = −D −1(L + U )x + D −1b x(k +1) = −D −1 (L + U ) x(k ) + D −1b
输入 A，b 和初始向量 x
迭代矩阵 BJ , BG
否
ρ(B) < 1?
按雅各比方法进行迭代
否
|| x (k+1) − x(k) ||< ε ?
按高斯-塞德尔法进行迭代
否
|| x(k+1) − x (k ) ||< ε ?
输出迭代结果
图 1 雅各布和高斯-赛德尔算法程序流程图
1.2 问题求解
按图 1 所示的程序流程，用 MATLAB 编写程序代码，具体见附录 1。解上述三个问题 如下
16
-0.72723528355328
0.80813484897616
0.25249261987171
17
-0.72729617968010
0.80805513082418
0.25253982509100
18
-0.72726173942623
0.80809395746552
0.25251408253388
0.80756312717373
8
-0.72715363032573
0.80789064377799
9
-0.72718652854079
MATLAB期末上机试题带答案
MATLAB 期末上机考试试题带答案版姓名： 学号： 成绩：1.请实现下图：xyy=sin(x)x=linspace(0,8*pi,250); y=sin(x); plot(x,y) area(y,-1) xlabel('x') ylabel('y') title('y=sin(x)') 2.请实现下图：x=linspace(0,2*pi,100); y1=sin(x);subplot(2,2,1)plot(x,y1,'k--')grid onxlabel('x')ylabel('y')title('sin(x)')legend('y=sin(x)')y2=cos(x);subplot(2,2,2)plot(x,y2,'r--')grid onxlabel('x')ylabel('y')title('cos(x)')legend('y=cos(x)')y3=tan(x);subplot(2,2,3)plot(x,y3,'k-')grid onxlabel('x')ylabel('y')title('tan(x)')legend('y=tan(x)')y4=cot(x);subplot(2,2,4)plot(x,y4)grid onxlabel('x')ylabel('y')title('cot(x)')legend('y=cot(x)')3.解方程组：a=[3 2 1;1 -1 3;2 4 -4]；b=[7;6;-2] ；x=a\b4.请实现下图：yxx=linspace(0,4*pi,1000);y1=sin(x);y2=sin(2*x);plot(x,y1,'--',x,y2,'b*')grid onxlabel('x');ylabel('y');title('耿蒙蒙')legend('sin(x)','sin(2*x)')5.请在x，y在（-2，2）内的z=xexp (-x2-y2) 绘制网格图[x,y]=meshgrid(-2:0.1:2);z=x.*exp (-x.^2-y.^2);mesh(x,y,z)6.请实现peaks函数：-55xPeaksy[x,y]=meshgrid(-3:1/8:3); z=peaks(x,y); mesh(x,y,z) surf(x,y,z) shading flataxis([-3 3 -3 3 -8 8])xlabel('x');ylabel('y');title('Peaks')7.请在x=[0,2],y=[-0.5*pi,7.5*pi],绘制光栅的振幅为0.4的三维正弦光栅。
2010专第一次上机作业： MatLab练习题-基本运算
一、书上命令练习：二、MatLab 练习1注意: 1,请按要求作题,写出程序并把结果附上,并保存文件.2,仔细填写 “学生信息”.3,把此word 文档的文件名更名为你的 “学号 姓名(*)”-*表示第几次作业.在下课前5分钟提交.1, (5*2+1.3-0.8)*10^2/25答案：2，y =23log )10sin(4)43.0(2+⋅+-eπ 答案：3，z = 10sin(π/3) ...cos(π/3)答案：4，x = sin(223/3)， y = x^2，z = y*10 ；求x+2y-5z答案：5，i e z i z i z ⋅==+=π23215,2,72,计算3221z z z z z += 答案：6, 建立起始值=3,增量值=5.5,终止值=44的一维数组x(填写程序语句)答案：7, 建立等差一维数组x ：首项为0,末项为π,项数为15 (填写程序语句) 答案：8, 计算行列式的值187624323=A答案：9, 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=187624323A ,矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333222111B ;分别求出A x B 及A 与B 中对应元素之间的乘积的程序语句。
答案10,写出计算 Sin(36o )的程序语句.答案：三、MatLab 练习2-矩阵及数组运算1.输入矩阵。
2.输入一个与A同阶的随机矩阵B。
要求元素为整数。
3.输入一个数组b=（1 3 5 7 9 11 ）4.输入与A同阶的单位矩阵E5.输入与B同阶的零矩阵Z6.构建一个大矩阵7.将 46-70这25个整数填入一个五行五列的矩阵数表G中，使其各行、各列以及主对角线和次对角线的和等相等。
8.分别求d1=|A|、d2=|B|、d3=|D|、d4=|G|9.对以上四个矩阵A、B、D、G中的可逆矩阵求逆。
10.(1) 计算d5=A+B(2) 计算d6=A*B11.取A的1，3，5行做一个子阵a1.12.取A的2，4，6行，4，5，6列上的交叉元素做一个子阵a2.13.由a2产生列向量a3.14.取a3绝对值大于4的元素构成列向量a4,15.改变A的元素A（4，5）为0 值。
MATLAB上机实验实验报告
MATLAB上机实验一一、实验目的初步熟悉MATLAB 工作环境，熟悉命令窗口，学会使用帮助窗口查找帮助信息。
命令窗口二、实验内容(1) 熟悉MATLAB 平台的工作环境。
(2) 熟悉MATLAB 的5 个工作窗口。
(3) MATLAB 的优先搜索顺序。
三、实验步骤1. 熟悉MATLAB 的5 个基本窗口①Command Window (命令窗口)②Workspace (工作空间窗口)—③Command History (命令历史记录窗口)④Current Directory (当前目录窗口)⑤Help Window (帮助窗口)(1) 命令窗口(Command Window)。
在命令窗口中依次输入以下命令：>>x=1>> y=[1 2 34 5 67 8 9];>> z1=[1:10],z2=[1:2:5];>> w=linspace(1,10,10);>> t1=ones(3),t2=ones(1,3),t3=ones(3,1)>> t4=ones(3),t4=eye(4)x =1z1 =1 2 3 4 5 6 7 8 9 10 t1 =1 1 11 1 11 1t2 =1 1 1t3 =111t4 =1 1 11 1 11 1 1t4 =1 0 0 00 1 0 00 0 1 00 0 0 1思考题：①变量如何声明，变量名须遵守什么规则、是否区分大小写。
答：（1）变量声明1.局部变量每个函数都有自己的局部变量，这些变量只能在定义它的函数内部使用。
当函数运行时，局部变量保存在函数的工作空间中，一旦函数退出，这些局部变量将不复存在。
脚本（没有输入输出参数，由一系列MATLAB命令组成的M文件）没有单独的工作空间，只能共享调用者的工作空间。
当从命令行调用，脚本变量存在基本工作空间中；当从函数调用，脚本变量存在函数空间中。
2.全局变量在函数或基本工作空间内，用global声明的变量为全局变量。
数值分析作业-matlab上机作业
数值分析———Matlab上机作业学院：班级：老师：姓名：学号：第二章解线性方程组的直接解法第14题【解】1、编写一个追赶法的函数输入a，b，c，d输出结果x，均为数组形式function x=Zhuiganfa(a,b,c,d)%首先说明：追赶法是适用于三对角矩阵的线性方程组求解的方法，并不适用于其他类型矩阵。
%定义三对角矩阵A的各组成单元。
方程为Ax=d%b为A的对角线元素(1~n)，a为-1对角线元素(2~n)，c为+1对角线元素(1~n-1)。
% A=[2 -1 0 0% -1 3 -2 0% 0 -2 4 -3% 0 0 -3 5]% a=[-1 -2 -3];c=[-1 -2 -3];b=[2 3 4 5];d=[6 1 -2 1];n=length(b);u(1)=b(1);y(1)=d(1);for i=2:nl(i)=a(i-1)/u(i-1);%先求l(i)u(i)=b(i)-c(i-1)*l(i);%再求u(i)%A=LU,Ax=LUx=d,y=Ux,%Ly=d,由于L是下三角矩阵，对角线均为1，所以可求y(i)y(i)=d(i)-l(i)*y(i-1);endx(n)=y(n)/u(n);for i=(n-1):-1:1%Ux=y,由于U是上三角矩阵，所以可求x(i)x(i)=(y(i)-c(i)*x(i+1))/u(i);end2、输入已知参数>>a=[2 2 2 2 2 2 2];>>b=[2 5 5 5 5 5 5 5];>>c=[2 2 2 2 2 2 2];>>d=[220/27 0 0 0 0 0 0 0];3、按定义格式调用函数>>x=zhuiganfa(a,b,c,d)4、输出结果x=[8.147775166909105 -4.073701092835030 2.036477565178471 -1.017492820111148 0.507254485099400 -0.250643392637350 0.119353996493976 -0.047741598597591]第15题【解】1、编写一个程序生成题目条件生成线性方程组A x=b 的系数矩阵A 和右端项量b ，分别定义矩阵A 、B 、a 、b 分别表示系数矩阵，其中1(10.1;,1,2,...,)j ij i i a x x i i j n -==+=或1(,1,2,...,)1ij a i j n i j ==+-分别构成A 、B 对应右端项量分别a 、b 。
matlab上机练习（附答案）
1.以下两种说法对吗？（1）MATLAB进行数值的表达精度与其指令窗中的数据显示精度相同。
（2）MATLAB指令窗中显示的数据有效位数不超过七位。
2.历史指令窗所记录的内容与diary指令所产生的“日志”内容有什么不同？DIARY filename causes a copy of all subsequent command window inputand most of the resulting command window output to be appended to thenamed file. If no file is specified, the file 'diary' is used.DIARY OFF suspends it.DIARY ON turns it back on.DIARY, by itself, toggles the diary state.Use the functional form of DIARY, such as DIARY('file'),when the file name is stored in a string.3.如何把用户自己的“工作目录”永久地设置在MATLAB的搜索路径上？“位于搜索路径上的目录”与“当前目录”在MATLAB中的功用相同吗？4.如何向MATLAB工作空间输入一个含有100个左右元素的一维或二维数值数组？用直接键入法？用数组编辑器？用M文件编辑器？5.运用数组算术运算符去掉下面程序里的for/end循环：x=11:15for k=1:length(x)z(k)=x(k)^2+2.3*x(k)^0.5;endx=11:15 z1=x.^2+2.3*x.^0.56.不使用数组算术运算符，重写下面的程序代码：x=[2 1 4]z=1./(1+x.^2)x=2;k=1;while i<=4,z2(k)=1/(1+i^2);i=i+1;x=x+1;end7．某公司销售电脑打印机的价格方案如下：（）如果顾客只买一台打印机，则一台的基本价格为$150。
MATLAB上机作业
MA TLAB上机作业学院名称：专业班级：学生姓名：学生学号：年月作业11.用MATLAB可以识别的格式输入下面两个矩阵12342357135732391894A ⎛⎫⎪⎪ ⎪= ⎪⎪ ⎪⎝⎭,144367723355422675342189543i iB i +⎛⎫⎪+⎪= ⎪+ ⎪ ⎪⎝⎭再求出它们的乘积矩阵C ,并将C 矩阵的右下角23⨯子矩阵赋给D 矩阵。
赋值完成后,调用相应的命令查看MATLAB 工作空间的占有情况。
解：A=[1 2 3 4;2 3 5 7;1 3 5 7 ;3 2 3 9 ;1 8 9 4;]B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;] B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;] C=A*B D=C(4:5,4:6); whos; 2.设矩阵16231351110897612414152⎛⎫ ⎪⎪ ⎪ ⎪⎝⎭,求A ,1A -,3A ,12A A -+,1'3A A --,并求矩阵A 的特征值和特征向量。
解：A=[16 2 3 13;5 11 10 8;9 7 6 12;4 14 15 2;] det(A) inv(A) A.^3 2*A+inv(A) 3*A-A' [V,D]=eig(A)abs (A)3.解下列矩阵方程: 01010014310001201001010120X -⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪=-⎪ ⎪⎪⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭ 解：A=[0 1 0;1 0 0;0 0 1;]; B=[1 0 0;0 0 1;0 1 0;]; C=[1 -4 3; 2 0 -1;1 -2 0;]; X=inv(A)*C*inv(B)4.一球从100米高度自由落下,每次落地后反跳回原高度的一半,再落下。
matlab上机习题答案
matlab上机习题答案Matlab上机习题答案在现代科学和工程领域中，计算机编程和数值计算已经成为必不可少的技能。
而Matlab作为一种强大的数值计算软件，被广泛应用于各种领域。
为了帮助学习者更好地掌握Matlab的基本操作和数值计算方法，老师们经常会布置一些上机习题，让学生通过实际操作来加深对Matlab的理解。
下面我们来看一些常见的Matlab上机习题答案：1. 编写一个Matlab程序，计算并输出1到100之间所有奇数的和。
答案：```matlabsum = 0;for i = 1:2:100sum = sum + i;enddisp(sum);```2. 编写一个Matlab程序，计算并输出斐波那契数列的前20个数字。
答案：```matlabfib = zeros(1,20);fib(1) = 1;fib(2) = 1;for i = 3:20fib(i) = fib(i-1) + fib(i-2);enddisp(fib);```3. 编写一个Matlab程序，求解一元二次方程ax^2 + bx + c = 0的根。
答案：```matlaba = 1;b = -3;c = 2;delta = b^2 - 4*a*c;if delta < 0disp('无实根');elseif delta == 0x = -b / (2*a);disp(x);elsex1 = (-b + sqrt(delta)) / (2*a);x2 = (-b - sqrt(delta)) / (2*a);disp(x1);disp(x2);end```通过以上几个例子，我们可以看到，Matlab的语法简洁明了，功能强大。
通过编写程序来解决实际问题，不仅加深了对Matlab的理解，也提高了计算机编程和数值计算的能力。
希望大家在学习Matlab的过程中能够多多练习，不断提高自己的编程水平。
(完整版)matlab上机练习题答案
(完整版)matlab 上机练习题答案1.计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积>〉 a=［6 9 3;2 7 5]; 〉〉 b=[2 4 1;4 6 8］; 〉〉 a 。
*b ans =12 36 3 8 42 402。
对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
〉〉 A=［4 9 2；7 6 4；3 5 7]; 〉> B=[37 26 28］’； >〉 X=A\B X = -0.5118 4.0427 1.33183。
⎥⎦⎤⎢⎣⎡-=463521a ，⎥⎦⎤⎢⎣⎡-=263478b ,观察a 与b 之间的六种关系运算的结果 >〉 a=[1 2 3;4 5 6]; 〉> b=［8 –7 4；3 6 2]； >〉 a 〉b ans =0 1 0 1 0 1 〉> a 〉=b ans =0 1 0 1 0 1 >> a 〈b ans =1 0 1 0 1 0 〉> a<=b ans =1 0 1 0 1 0 >〉 a==b ans =0 0 0 0 0 0 〉〉 a~=b ans =1 1 1 1 1 14计算多项式乘法（x 2+2x +2)（x 2+5x +4)>> c=conv （[1 2 2]，[1 5 4］) c =1 7 16 18 8 5计算多项式除法（3x 3+13x 2+6x +8）/(x +4） 〉〉 d=deconv （［3 13 6 8],［1 4］) d =3 1 26求欠定方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡5865394742x 的最小范数解〉〉 a=［2 4 7 4;9 3 5 6］； 〉> b=[8 5］’; 〉> x=pinv （a ）*b x =—0.2151 0.4459 0。
Matlab上机题及答案
1 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。
输出全部水仙花数。
for m=100:999m1=fix(m/100); %求m的百位数字m2=rem(fix(m/10),10); %求m的十位数字m3=rem(m,10); %求m的个位数字if m==m1*m1*m1+m2*m2*m2+m3*m3*m3disp(m)endend2.从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。
sum=0;n=0;val=input('Enter a number (end in 0):');while (val~=0)sum=sum+val;n=n+1;val=input('Enter a number (end in 0):');endif (n > 0)summean=sum/nend3. 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。
求[1,500]之间的全部完数。
for m=1:500s=0;for k=1:m/2if rem(m,k)==0s=s+k;endendif m==sdisp(m);endend4. 从键盘上输入数字星期，在屏幕上显示对应英文星期的单词。
function weekn=input('input the number:');if isempty(n)errror('please input !!')endif n>7|n<1error('n between 1 and 7')endswitch ncase 1disp('Monday')case 2disp('Tuesday')case 3disp('Wednesday')case 4disp('Thursday')case 5disp('Friday')case 6disp('Saturday')case 7disp('Sunday')end5. 某公司销售电脑打印机的价格方案如下：（）如果顾客只买一台打印机，则一台的基本价格为$150。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
MATLAB 上机作业学院名称： 专业班级： 学生姓名： 学生学号：年 月作业11.用MATLAB 可以识别的格式输入下面两个矩阵12342357135732391894A ⎛⎫⎪⎪ ⎪= ⎪⎪ ⎪⎝⎭,144367723355422675342189543i i B i +⎛⎫⎪+⎪= ⎪+ ⎪⎪⎝⎭再求出它们的乘积矩阵C ,并将C 矩阵的右下角23⨯子矩阵赋给D 矩阵。
赋值完成后,调用相应的命令查看MATLAB 工作空间的占有情况。
解：A=[1 2 3 4;2 3 5 7;1 3 5 7 ;3 2 3 9 ;1 8 9 4;]B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;] B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;] C=A*B D=C(4:5,4:6); whos;2.设矩阵16231351110897612414152⎛⎫⎪⎪⎪⎪⎝⎭,求A,1A-,3A,12A A-+,1'3A A--,并求矩阵A的特征值和特征向量。
解：A=[16 2 3 13;5 11 10 8;9 7 6 12;4 14 15 2;]det(A)inv(A)A.^32*A+inv(A)3*A-A'[V,D]=eig(A)abs (A)3.解下列矩阵方程:解：A=[0 1 0;1 0 0;0 0 1;];B=[1 0 0;0 0 1;0 1 0;];C=[1 -4 3; 2 0 -1;1 -2 0;];X=inv(A)*C*inv(B)4.一球从100米高度自由落下,每次落地后反跳回原高度的一半,再落下。
求它在第10次落地时,共经过多少米?第10次反弹有多高?syms x y sum ;x=100;y=0.5*x;sum=0;for i=2:10, sum=sum+x+y; x=0.5*x; y=0.5*x; end sum y5.用MATLAB 语言实现下面的分段函数 解：function y=f(x) if x>10 y=5; elseif x>=-10&x<=10 y=1/2*x; else y=-5; end end6.分别用for 和while 循环编写程序,求出并考虑一种避免循环的简洁方法来进行求和,并比较各种算法的运行时间。
syms x K ; x=2; K=0;for i=1:63, K=K+x^i; end K7.应用MA TLAB 语言及二分法编写求解一元方程32()1459700f x x x x =-+-=在区间[3，6]的实数解的算法，要求绝对误差不超过0.001。
format longA=[1 -14 59 -70]; a=3; b=6; c=0.001;while c<0.5*(b-a), x1=(a+b)/2; f1=polyval(A,x1);fa=polyval(A,a); fb=polyval(A,b); if f1==0; x=x1; else if f1*fa>0; a=x1;else f1*fb>0; b=x1; end end x=x18.二阶系统的单位阶跃响应为1cos )t y a ζζ-=+，在同一平面绘制ζ分别为0，0.3，0.5，0.707的单位阶跃响应曲线。
要求：(1)四条曲线的颜色分别为蓝、绿、红、黄，线型分别为“——”、“……”、“oooooo”、“++++++”；(2)添加横坐标轴和纵坐标轴名分别为“时间t”和“响应y”，并在平面图上添加标题“二阶系统曲线”和网格；(3)在右上角添加图例（即用对应的字符串区分图形上的线），并分别在对应的响应曲线的第一个峰值处标示“zeta ＝0”、“zeta ＝0.3”、“zeta ＝0.5”、“zeta ＝0.707”。
t=[0:0.1:10]; kos=0;y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos)); plot(t,y,'b-') hold on; kos=0.3;y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos)); plot(t,y,'g.') hold on; kos=0.5;y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos)); plot(t,y,'ro') hold on; kos=0.707;y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos)); plot(t,y,'y+') hold on;xlabel('时间t'); ylabel('响应y'); grid on;title('二阶系统曲线');legend('kos=0','kos=0.3','kos=0.5','kos=0.707') gtext('kos=0')gtext('kos=0.3') gtext('kos=0.5') gtext('kos=0.707')9．绘制如下图所示的图形，把图形窗口分割为2行2列，窗口左上角画一正弦曲线sin(2),[0,2]y t t π=∈；窗口右上角画3条单边指数曲线23,,,[0,2]t t t y e y e y e t ---===∈；窗口左下角画一矩形脉冲信号，脉冲宽度为1，高为2，开始时间为1；窗口右下角画一单位圆。
t=0:0.01:2;subplot(2,2,1);plot(t,sin(2*pi*t));title('plot(x,2*PI*t)');grid on; subplot(2,2,2);plot(t,exp(-t));hold on;plot(t,exp(-2*t),'g');hold on; plot(t,exp(-3*t),'r');hold on;title('exp(-t),exp(-2*t),exp(-3*t)');grid on;x=[0 1 1 2 2 3 4]; y=[0 0 2 2 0 0 0]; subplot(2,2,3); plot(x,y);axis([0 4 -0.5 3]); title('pulse signal');syms x y subplot(2,2,4); ezplot(x^2+y^2-1); axis([-1.5 1.5 -1.2 1.2]); title('circle'); xlabel('')ylabel('') 10.已知函数(,)z f x y ==，试分别应用三维曲线图绘制命令plot3、三维网线图绘制命令mesh 、三维曲面图绘制命令surf 在同一窗口中绘制出3个子图。
x=-2:0.1:2;[X,Y]=meshgrid(x);Z=1./sqrt((1-X).^2+Y.^2)+1./sqrt((1+X).^2+Y.^2); subplot(311); plot3(X,Y,Z); grid on;axis([-1 1 -1 1 -1 9]);subplot(312);mesh(X,Y,Z);axis([-1 1 -1 1 -1 9]); subplot(313); surf(X,Y,Z);axis([-1 1 -1 1 -1 9]); 11.对合适的θ范围选取分别绘制出下列极坐标图形:(1)cos(7/2)ρθ= (2)31cos (7)ρθ=-clctheta=0:0.01:2*pi;rh o1=cos(7*theta./2);pol ar(theta,rho1),title('p=c os(7o/2)') figure(2);rho2=1-cos(7*theta).^3;polar(theta,rho2),title('p=1-(cos(7o))^3')12.求解下面两个方程构成的联立方程组在区间[2,6],[3,5]x y ∈-∈-内的解，并用绘图的方法绘出两曲线在同一坐标上的图，以验证求得的解的正确性。
解22(2)3(1)3x y -+-=, 26(2)y x =- clc[x,y]=solve('(x-2)^2+3*(y-1)^2-3=0','y-6*(x-2)^2=0') ezplot('(x-2)^2+3*(y-1)^2-3');hold on ezplot('y-6*(x-2)^2=0');grid on axis([0 4 -1 4]);作业21.求多项式322()(5713)(53)(3)f s s s s s s s =++++++的乘积并求()f s 的导数。
解A=[1 5 7 13]; B=[1 5 3]; C=[1 3];f=conv(conv(A,B),C); diff(f)2.求当3x =时多项式322()(53513)(453)(313)f x x x x x x x =++++++的值。
解A=[1 5 35 13]; B=[2 45 3]; C=[31 3];f=conv(conv(A,B),C); polyval(f,3)3.试用不同的方法展开多项式23242()(1)(5)(47)P s s s s s =++++,并比较其结果。
解A=[1 0 1]; B=[1 5]; C=[1 0 4 7];f=conv(conv(A.^3,B.^2),C) 4.求下列多项式的根和导数 1)326154g x x x =-+-2)432514103h x x x x =---- 解（1）g=[1 -6 15 -4]; r=roots(g) diff(g)（2）h=[1 -5 -14 -10 -3]; r=roots(h) diff(h)5.对于有理多项式 42210(4567)()(1)(2)(3)s s s s f s s s s ++++=+++ (1)计算该多项式相除的结果;(2)将该多项式展开为部分分式的形式;(3)计算df ds。