MATLAB期末上机试题带答案
matlab期末考试试卷及参考答案
一、填空题(本题满分30分，每空3分)1.设有程序A=[2,0,1,0;1,0,1,5;0, 6 ,2, 3];B=[1,2,-1;1,3,1];C=A(1:3,[1,2]);D=size(C)- length(B(2,:));F= C.*B ’将下列命令的运行结果填在横线上0*、0D = [0 -1] ; F= [2 0;2 0;0 6] .2. A=[2,3,2,1;1,2,2,1; 1,2,3,-2];A1=sum(A>A(2,1));A(2,:)=[],A2=A.^2-4A1= ; A2= ;3.P=[1,2,5,4;2,2,1,3];max(P,[ ],1)= [2 2 5 4] ;mean(P,2)= [3 2] .二、(本题满分10分)试编写计算程序画出曲线与曲面图形.22,22,)1((2).);20(cos sin (1).223322≤≤-≤≤-++=≤≤⎪⎩⎪⎨⎧==--y x y x e x z t t y t x y xπ%1 %2t=0:0.1: 2*pi; [x,y]=meshgrid(-2:0.1:2);x= sin(t).^3; z=(x.^2+1).* exp(-x.^2- y.^2+)+x .^2.*y; y= cos (t).^3; mesh(x,y,z)plot(x,y)三、（本题满分12分）编写程序解决下列问题：（1）试产生100⨯4阶矩阵使其元素在0—100之间的随机整数,可将此矩阵理解为100名同学的四门课程(课程1, 课程2, 课程3, 课程4)的成绩.（2）计算出每门课程和每位同学的平均成绩;（3）统计出总的优秀人次(90分及以上为优秀);（4）统计出四门课程成绩至少有两门课程不及格的人数.解：(1) a=fix((100-0+1)*rand(100,4));(2) mean(a) mean(a,2)(3)sum(sum(a>=90,2)==4)(4)sum(sum(a<60,2)>=2)四、（本题满10分）).50,20(]2,0[)2()),(1:(;)1(,sin )1()(11==-=∑=+n n x fun y function kx k x f nk k 分别取内图形用子图画出该函数在程序的第一句为返回该函数并保存试编写函数式文件程序设π(1). function y=fun1(x,n) (2). x=0:0.1:2*pi;y=0; y1= fun1 (x,20);for k=1:n y2= fun1 (x,50);y=y+ ((-1)^(k+1)*sin(k*x))/k; subplot(2,1,1),plot(x,y1)end subplot(2,1,2),plot(x,y2)y五、(本题满分10分)试利用微分方程的数值解法写出下列微分方程的求解步骤和MATLAB 程序.0)0(,2)0(,3='==-'+''x x t x x t x .(t 的变化区间为[0,10]).第一步：先将高阶微分方程转化为一阶微分方程组：选择状态变量x x x x ='=21,，则原方程化为2)0(,0)0(.,32112211==⎩⎨⎧='++-='x x x x t x tx x第二步：建立函数文件function dxdt = fun2(t,x)dxdt = [-3*t*x(1)+ x(2)+t; x(1)];第三步：求解微分方程，命令如下：[t,x]=ode45(@fun2,[0,10],[0;2]);)(cos cos lim .3;,)(,,1ln sin .2;2.1:,sin 022x e x x x A dx x A dx dA x s e x x x A x e 、x x x x --⎥⎥⎦⎤⎢⎢⎣⎡+=+=→-⎰计算极限并计算生成符号矩阵的根求方程下列各题运用符号运算功能计算六 的解函数求解微分方程用1)0()0()0()0(,2.4)4(='''=''='==+x x x x e x x dsolve t六、参考程序：1． solve('exp(x)=x^2+2')2．syms x sA=[x*sin(x),log(x);exp(-x),1/(s+x)]diff(A)int(A)A*A3．syms xlimit((cos(x)-cos(x)^(1/2))/x/(exp(sin(x))-1))4．dsolve('D4x+x=2*exp(t)','x(0)=1','Dx(0)=1','D2x(0)=1','D3x(0)=1')⎪⎩⎪⎨⎧=-++=--++03320224:.2122212212221x x x x x x x x x fsolve 组求解下面的非线性方程函数以及符号运算功能用七解一：function y=fc(x)y(1)=4*x(1).^2+x(2).^2+2*x(1).*x(2)-x(2)-2y(2)=2*x(1).^2+x(2).^2+3*x(1).*x(2)-3y=[y(1) y(2)];x0=[1 1];fsolve(@fc,x0)解二: syms x1 x2[x1 x2]=solve(‘4*x1^2+ x2^2+2* x1* x2- x2-2=0’,’ 2*x1^2+ x2^2+3* x1* x2-3=0’,’ x1’,’ x2’)《MA TLAB程序设计》参考程序一、(30分)1．D =0 -1F =2 02 00 62．A1 =1 3 3 0A2 =0 5 0 -3-3 0 5 03．max(P,[ ],1)=2 2 5 4 mean(P,2)=324．x=155．b =0 1 4-1 0 1-4 -1 06．z =20107．285二、(10分)%1t=0:0.1: 2*pi;x= sin(t).^3;y= cos (t).^3;plot(x,y)%2[x,y]=meshgrid(-2:0.1:2);z=(x.^2+1).* exp(-x.^2- y.^2+)+x.^2.*y; mesh(x,y,z)三、(16分)syms xlimit((x-1)*log(x)/(exp(x-1)-1),x,1)%2solve('exp(x)-3*cos(x)+1=0')%3syms xM=[ exp(-x), x;sin(x),exp(-0.1*x)* cos(3*x)]diff(M,2)int(M,x,0,pi)%4 syms t xdsolve('Dx-6*x=exp(-t) ','x(0)=1')四、(12分)x=0.5:9.5;y=[ 0.16, 0.40, 0.66, 0.75, 0.88, 1.0, 1.1, 1.1, 1.2, 1.3];x1=1./x;y1=1./y;p=polyfit(x1,y1,1);a=p(2); b= p(1);ye=x./(a*x+b);plot(x,y,'ro', x, ye, 'b-')grid onxlabel('自变量x ')ylabel('因变量y ')title('因变量y与自变量x的拟合曲线图')五、(12分)%参考程序cj=round(normrnd(70,10,100,4));for i=1:100for j=1:4if cj(i,j)>100 cj(i,j)=100;elseif cj(i,j)<0 cj(i,j)=0;endendendkmean=round(mean(cj))disp('课程1 课程2 课程3 课程4')disp(sprintf(' %4d%7d%8d%8d ' ,kmean(1),kmean(2),kmean(3),kmean(4))) studentmean=round(mean(cj,2));studentnumber=[1:100]';[studentnumber, studentmean]c=sum(sum(cj>=90))d=sum(sum(cj<60,2)>=2)六、(10分)(1). function y=fun1(x,n)y=0;for k=1:ny=y+ ((-1)^k*sin(k*x))/k;endy(2). x=0:0.1:2*pi;y1= fun1 (x,20);y2= fun1 (x,50);subplot(211),plot(x,y1)subplot(212),plot(x,y2)七、(10分)第一步：先将高阶微分方程转化为一阶微分方程组：选择状态变量x x x x ='=21,，则原方程化为2)0(,0)0(.,32112211==⎩⎨⎧='++-='x x x x t x tx x第二步：建立函数文件function dxdt = fun2(t,x)dxdt = [-3*t*x(1)+ x(2)+t; x(1)];第三步：求解微分方程，命令如下：[t,x]=ode45(@fun2,[0,10],[0;2])补充：1.matlab 中清屏的命令是———clc2.matlab 中用于清除工作空间变量的命令是——clear3.在【0,5】上求函数f=（x-3）^2-1的最小值。
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期末考试题及答案
matlab期末考试题及答案MATLAB期末考试题及答案一、选择题（每题2分，共20分）1. MATLAB中用于创建向量的函数是：A. vectorB. arrayC. linspaceD. ones答案：D2. 下列哪个命令可以计算矩阵的行列式？A. detB. diagC. traceD. rank答案：A3. 在MATLAB中，以下哪个选项是用于绘制三维图形的？A. plotB. plot3C. barD. scatter答案：B4. MATLAB中，用于计算向量范数的函数是：A. normB. meanC. medianD. std答案：A5. 下列哪个命令可以用于创建一个二维数组？A. array2dB. matrixC. create2dD. make2d答案：B6. MATLAB中，用于求解线性方程组的函数是：A. solveB. linsolveC. equationD. linprog答案：A7. 以下哪个函数可以用于生成随机数？A. randB. randomC. randnD. randi答案：A8. MATLAB中，用于实现循环结构的关键字是：A. loopB. forC. whileD. repeat答案：B9. 下列哪个命令可以用于绘制函数图形？A. plotB. graphC. drawD. functionplot答案：A10. MATLAB中，用于计算矩阵特征值的函数是：A. eigB. eigenvalueC. characteristicD. eigen答案：A二、简答题（每题5分，共30分）1. 简述MATLAB中矩阵的基本操作有哪些？答案：矩阵的基本操作包括矩阵的创建、矩阵的加法、减法、乘法、转置、求逆、求行列式等。
2. MATLAB中如何实现条件语句？答案：MATLAB中实现条件语句主要使用if-else结构，也可以使用switch-case结构。
3. 请解释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期末考试试卷汇总
MATLAB期末考试试卷及其参考答案一、填空题1、MATLAB常用操作界面包括、工作空间窗口、、、内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。
2、MATLAB中Inf或inf表示、NaN或nan表示、nargout表示。
3、MATLAB中逗号主要用作；用作输入量与输入量之间的分隔符；。
4、工作空间浏览器主要用于内存变量的、和。
5、MATLAB实现将全下标转换为单下标的指令为、据单下标换算出全下标的指令为。
6、二维数组的标识有、、“逻辑1”标识。
二、简答题1、简述MATLAB历史指令窗的主要作用。
2、简述空数组的功用。
3、简述MATLAB函数的基本结构。
4、简述绘制二维图形的一般步骤。
三、阅读程序并回答问题（每题4分，共28分）1、写出下列指令运行结果。
A=zeros(2,4);A(:)=1:8;s=[2 3 5];A(s)Sa=[10 20 30]'A(s)=Sa2、写出下列指令运行结果。
A=reshape(1:16,2,8)reshape(A,4,4)s=[1 3 6 8 9 11 14 16];A(s)=03、写出下列指令运行结果。
A=[1,2;3,4];B=[-1,-2;2;1];S=3;A.*BA*BS.*AS*B4、下面的函数主要完成什么功能？function f=factor(n)if n<=1f=1;elsef=factor(n-1)*n;end5、写出下列指令运行结果。
ch=‘ABc123d4e56Fg9’;subch=ch(1:5)revch=ch(end:-1:1)k=find(ch>=‘a’&ch<=‘z’);ch(k)=ch(k)-(‘a’-‘A’);char(ch)6、写出下列指令运行结果。
A(1,1)={'this is cell'};A{1,2}={[1 2 3;4 5 6]};A{2,1}=[1+2*i];A{2,2}=A{1,2}{1}+(A{1,2}{1}(1,1)+A{1,2}{1}(2,2));celldisp(A)7、下面的程序完成功能是什么？t=0:pi/50:4*pi;y0=exp(-t/3);y=exp(-t/3).*sin(3*t);plot(t,y,'-r',t,y0,':b',t,-y0,':b')xlabel(‘\bf\it t’);ylabel(‘\bf\it y’);grid on;四、编程题（32分）1、1、在同一图上分别用红色实线和绿色虚线绘制y1=sin(x)和y2=cos(x)在区间[0，4*pi]的曲线，并用星号*标出两条曲线的交点以及建立图例。
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。
MATLAB期末考试试卷及其参考答案.doc
一、填空题（每空1分，20分）1、MATLAB常用操作界面包括、工作空间窗口、、、内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。
2、MATLAB中Inf或inf表示、NaN或nan表示、nargout表示。
3、MATLAB中逗号主要用作；用作输入量与输入量之间的分隔符；。
4、工作空间浏览器主要用于内存变量的、和。
5、MATLAB实现将全下标转换为单下标的指令为、据单下标换算出全下标的指令为。
6、二维数组的标识有、、“逻辑1”标识。
、在一个元胞数组A中寻访第2行第3列元胞元素用；寻访数组第2行第3列元胞中所存的内容用。
8、4、MATLAB中clf用于、clc用于、clear用于。
二、简答题（每题5分，共20分）1、简述MATLAB历史指令窗的主要作用。
2、简述空数组的功用。
3、简述MATLAB函数的基本结构。
4、简述绘制二维图形的一般步骤。
三、阅读程序并回答问题（每题4分，共28分）、写出下列指令运行结果。
=zeros(2,4);(:)=1:8;=[2 3 5];(s)a=[10 20 30]'(s)=Sa、写出下列指令运行结果。
=reshape(1:16,2,8)eshape(A,4,4)=[1 3 6 8 9 11 14 16]; (s)=0、写出下列指令运行结果。
=[1,2;3,4];=[-1,-2;2;1];=3;.*B*B.*AS*B4、下面的函数主要完成什么功能？function f=factor(n)f n<=1f=1;lsef=factor(n-1)*n;nd、写出下列指令运行结果。
(1,1)={'this is cell'};{1,2}={[1 2 3;4 5 6]};{2,1}=[1+2*i];{2,2}=A{1,2}{1}+(A{1,2}{1}(1,1)+A{1,2}{1}(2,2)); elldisp(A)、下面的程序完成功能是什么？=0:pi/50:4*pi;0=exp(-t/3);=exp(-t/3).*sin(3*t);lot(t,y,'-r',t,y0,':b',t,-y0,':b')label(‘\bf\it t’);label(‘\bf\it y’);rid on;四、编程题（32分）price<200 没有折扣200≤price<500 3%折扣500≤price<1000 5%折扣1000≤price<2500 8%折扣2500≤price<5000 10%折扣000≤price 14%折扣入所售商品的价格，求其实际销售价格。
大学专业试题-Matlab期末考试及答案
大学专业试题-Matlab期末考试及答案一、填空题（一题3分，共15分）1、标点符号可以使命令行不显示运算结果。
答案：%2、x为0〜4，步长为0.1pi的向量，使用命令创建。
答案：x=0:0.1*pi:4。
3、输入矩阵A=[1 2 3;7 -5 8;9-1-7],使用全下标方式用A（2,2）取出元素“-5 ”，使用单下标方式用取出元素“-5 ”。
答案：A （5）。
4、在while表达式,语句体,End循环语句中，表达式的值__时表示循环条件为真，体语句将被执行，否则跳出该循环语句。
答案：非零。
5、要从键盘读入一个字符串并赋值给变量x,且给出提示“Hello"，应使用命令—答案：x=input（‘Hello,；s'）b二、选择题（一题3分，共15分）1、下列表示可以作为MATLAB的变量的是(D )A. abcd-2B. xyz_2#C. @hD. X_1_a2、如果x=1: 2 : 8,则U x(1)和x(4)分别是(B )A. 1, 8B. 1,7C. 2,8D. 2,73、指出下列错误的指令（B）。
A.symsa b;B. symsa, b;C.syms('a',’b');D. syms('a','b','positive');4、下列哪条指令是求矩阵的行列式的值(C)。
A .inv B. diag C. det D.eig5、清空Matlab工作空间内所有变量的指令是(C)。
A. clcB.clsC.clearD. clf三、程序分析题（一题10分，共20分）1、分析下面程序运行后s1,n1,s2,n2,m的值。
s1=0;s2=0;n1=0;n2=0；x=[1,-4,-8,3,10,-9,7,-3,10, 8,-5,-2,2,0 ]; m=length(x); fori=1:mif x(i)<0s1=s1+x(i); n1=n1+1; elses2=s2+x(i); n2=n2+1;end ends1,n1,s2,n2,m请回答s1,n1,s2,n2,m的结果。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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 on xlabel('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.请实现下图：xyx=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 (-x 2-y 2) 绘制网格图 [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的三维正弦光栅。
x=0:0.1:2;y=-0.5*pi:0.01*pi:7.5*pi;[x,y]=meshgrid(x,y)z=sin(y);mesh(x,y,z)surf(x,y,z)shading flataxis([0 2 -0.5*pi 7.5*pi -6 6])8.请用ezplot函数绘制sin（x2），x的区间为0到8*pi。
ezplot('x','sin(x^2)',[0,8*pi])9.样本点; x=[0 0.25*pi 0.5*pi 0.75*pi pi 1.25*pi 1.5*pi 1.75*pi 2*pi];y=[0 0.5*2.^0.5 1 0.5*2.^0.5 0 -0.5*2.^0.5 -1 -0.5*2.^0.5 0]; 对样本点进行spline插值。
并将样本点和插值后的数据进行绘图。
x0=[0 0.25*pi 0.5*pi 0.75*pi pi 1.25*pi 1.5*pi 1.75*pi 2*pi];y0=[0 0.5*2.^0.5 1 0.5*2.^0.5 0 -0.5*2.^0.5 -1 -0.5*2.^0.5 0]; plot(x0,y0,'o')hold onx=[0:0.01:2*pi];y=interp1(x0,y0,x,'spline')plot(x,y)10.请实现下图：利用多项式求根方法求解x^3-x^2-3=0。
y=[1 -1 0 -3];r=roots(y)11. A = 2 i + 5j + 7 kB = 8 i + 4j + 6 k求C=A×BA=[2 5 7];B=[8 4 6];C=conv(A,B)12. A = 2 i + 5j + 7 kB = 8 i + 4j + 6 k求C=A●BA=[2 5 7];B=[8 4 6];C=A*B'13. 用不同标度在同一坐标内绘制曲线y1=e-0.3x cos(2x)及曲线y2=10e-1.5x。
x=0:pi/180:2*pi;y1=exp(-0.3*x).*cos(2*x);y2=10*exp(-1.5*x);plotyy(x,y1,x,y2)14.请实现下图：x=linspace(0,8*pi,1400);y=sin(x);plot(x,y)area(y,0)15.n=dblquad('exp(-(x.^2)/2).*sin(x.^2+y)',-1,1,-2,2)16.请实现: t=0:0.1:10*pi,x=tcos(t),y=tsin(t),z=t,三维曲线。
ezplot3('t.*cos(t)','t.*sin(t)','t',[0,10*pi])17.a=[1 6 11 6]; r=roots(a) poly(r)18.syms xf=5*x^3+6*x^2+3*x+9; diff(f,x,1)19. 已知样本点x=[-2.8 -1 0.2 2.1 5.2 6.8]; y=[3.1 4.6 2.3 1.2 2.3 -1.1];求其三次拟合，并绘出样本点和拟合图像。
x=[-2.8 -1 0.2 2.1 5.2 6.8]; y=[3.1 4.6 2.3 1.2 2.3 -1.1]; plot(x,y,'o') hold onp=polyfit(x,y,3) x0=-3:0.01:7; y0=polyval(p,x0) plot(x0,y0,'r-')20.构建内联函数y=sin(x)exp(x^2);并求出x=[1 4 2 5 8]的y 值。
y=inline('sin(x).*exp(x.^2)');xi=[ 1 4 2 5 8]; yi=y(xi)21.请实现从距离地面20米高处，以水平速度5m/s 跳下的实际运动轨迹。
解：212h gt =-，0x v t =，得2202g h x v =-，即20.2h x =-（0<x<10m ）x=0:0.1:10;h=-0.2*x.^2;plot(x,h,'-') grid onxlabel('x/m');ylabel('h/m');title('ƽÅ×Ô˶¯') 22.请绘出斜抛运动的实际轨迹。
初速度为10m/s ，与地面的夹角为300。
解：2012y h gt v t =--,0x x v t =,得2115h x x =- x=0:0.1:13;h=-15^(-1)*x.^2-3^(-1/2)*x;plot(x,h) grid onxlabel('x/m');ylabel('h/m');title('бÅ×Ô˶¯') 23.请求出df(x)/dx=ax 3+x 2-bx-c syms x a b cf=a*x^3+x^2-b*x-c; diff(f,x,1)24.[x,y]=meshgrid(-3:0.1:3);z=1./((x+1).^2+(y+1).^2+1)-1./((x-1).^2+(y-1).^2+1);mesh(x,y,z)25.x=-10:0.01:10subplot(1,2,1)plot(x,sin(2*x).*cos(3*x))xlabel('x');ylabel('y');title('sin(2x)cos(3x)')subplot(1,2,2)plot(x,0.4*x)xlabel('x');ylabel('y');title('0.4x')26.x=0:0.01:25;y1=2.6*exp(-0.5*x).*cos(0.6*x)+0.8;y2=1.6*cos(3*x)+sin(x);plot(x,y1,'b--',x,y2,'r-')legend('y1=2.6exp(-0.5x)cos(0.6x)+0.8','y2=1.6cos(3x)+sin(x)') grid on27.y=int(sin(x)+2,x,0,pi/6)28.solve('sin(x)+tan(x)+1=0','x')29.syms xy=dsolve('Dy=(x+y)*(x-y)','x')30.解：插值法拟合法31. 请用三种方法求解sin(x)在0~pi之间的积分。
1：a=quad('sin(x)',0,pi)2：x=linspace(0,pi,1000);y=sin(x);a=trapz(x,y)3：y=int(sin(x),x,0,pi)32.[x,y]=meshgrid(-2:0.1:2);z=x.^2.*exp(-x.^2-y.^2); mesh(x,y,z)surf(x,y,z)33.22yxxez--=，当x和y的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗口中绘制出网线图、表面图和去网格效果的表面图。