05级通信工程、电子信息工程《MATLAB语言程序设计》考试卷
《MATLAB语言及应用》期末考试试题
c.三次多项式内插法(method='cubic')。
d.多项式拟合法：直接利用6次多项式去通过7个资料点。
请将这些内插的结果及原先的资料点画在同一个图。
共2页第1页
6. (10分)假设有一组实测数据
x
0.1
0.2
0.3
0.4
0.5
0.6
0.7
b.请用各种你可以想到的方法，在上述图形中加入各种解释和说明的文字。
4.(10分)试用roots指令算出x3+2x2+x+1=0的根。
5.(10分)假设一曲线资料点为
x = 0:2:4*pi;
y = sin(x).*exp(-x/5);
试将x的间距调成0.1，并用下列方法进行内插：
a.线性内插法(method='linear')。
0.8
0.9
1
y
2.3201
2.6470
2.9070
3.2885
3.6008
3.9090
4.2147
4.5191
4.8232
5.1275
(1)绘制出样条插值算法下的拟合结果；
(2)假设已知该数据可能满足的原型函数为 ，试求出满足下面数据的最小二乘解的值。
7.(10分)求下列联立方程的解
8.(10分)设
命题方式：
********学院2003—2004学年第二学期
《MATLAB语言及应用》课程期末考试试题
专业、班级：姓名：学号：
题号
一
二
三
四
五
六
七
八
《Matlab编程技术》考试题答案
要求：给出以下问题的matlab 实现方法（程序段）和运行结果。
1．求出满足∑=>mi i 1100000的最小m 值。
clc clear all sum=0; i=0;while sum<100000 i=i+1; sum=sum+i; end sum i Result ： sum =100128 i = 4472. 分别用循环结构和向量化方式求解级数求和问题i i i 21110012+∑=。
循环结构 clc clear allfor i=1:100 j=1/i^2+1/2^i; s=s+j; end 向量化： clc clear all a=zeros(1,100); for i=1:100a(i)=1/i^2+1/2^i; end s=sum(a) 结果s =2.63503. 试用函数plot 和函数fplot 绘制)tan(sin )sin(tan x x y -=在[]ππ,-∈x 区间内的曲线。
Plot ： clc clear all pi=3.14; for i=1:2*pi*100 u(i)=-pi+0.01*(i-1);y(i)=sin(tan(u(i)))-tan(sin(u(i))); end plot(u,y); Fplot ： clc clear all pi=3.14;y='sin(tan(x))-tan(sin(x))';-4-3-2-101234-3-2-11234. 若函数)cos(5)sin(),(2222y x e y y x x y x f y x ++++=+，试求其偏导数),(y x f y '，并绘制它们的曲线。
clear all clc syms x y pi=3.14;f=x^2*sin(x+y^2)+y^2*exp((x+y)+5*cos(x^2+y)); df=diff(f,y) 结果 df =2*x^2*cos(x+y^2)*y+2*y*exp(x+y+5*cos(x^2+y))+y^2*(1-5*sin(x^2+y))*exp(x+y+5*cos(x^2+y))x=0:pi/50:pi;y=x;f=x.^2.*sin(x+y .^2)+y .^2.*exp((x+y)+5.*cos(x.^2+y));df=2*x.^2.*cos(x+y .^2).*y+2*y .*exp(x+y+5*cos(x.^2+y))+y .^2.*(1-5.*sin(x.^2+y)).*exp(x+subplot(1,2,1) plot3(f,x,y) grid on title('f') subplot(1,2,2) plot3(df,x,y) grid on title('df')-5510x 10502401234-112x 106024012345.假设某概率密度函数有以下分段函数表示：⎪⎩⎪⎨⎧-≤+≤+<->+=+-------1,5457.011,7575.01,5457.0),(215.175.575.0216215.175.575.02112122212212122x x e x x e x x e x x p x x x x x x x x ，试以三维曲面的形式来表示这一函数，并在同一窗口显示该三维曲面的三视图。
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的最小值。
电子信息工程系2008级《MATLAB语言》试卷A参考答案
电子信息工程系2008级《MATLAB语言》试卷A参考答案1.Fill in the blanks ( 2’ * 5 = 10’ )(1)letter / 字母(2)0(3)imag( )(4)toc(5)eye (100, 102) 或[ eye (100), zeroes (100, 2) ]2.True or False ( 2’ * 5 = 10’)(1)T(2)T(3) F(4) F(5)T3.Short Ans wers ( 5’ * 3 = 15’ )(1)a) Add semicolons to the end of the statements.在句末加分号；b) Preallocate all the arrays used in loops.对所有循环体中用到的矩阵进行预分配空间；c) V ectorization向量化编程；d) Call a function by its handle.通过函数句柄调用函数。
注：答对一条给2分，答对两条给4分，答对三条给5分。
(2)a) Clearly state the problem that you are trying to solve.清晰的描述待解决的问题；b) Define the inputs required and the outputs to be produced by theprogram.定义程序需要的输入变量和返回的输出变量；c) Design the algorithm that you intend to implement in the program.设计程序中需要用到的算法；d) Turn the algorithm into MATLAB statements.将算法转化为MATLAB语句；e) Test the resulting MATLAB program.测试生成的MATLAB程序。
注：回答时只需回答出黑体部分的关键句即可，每条1分。
MATLAB期末考试试题(全12套)
MATLAB期末考试试题一、填空（30分）1. 表达式 (3>2)*(5~=5)的类型是（double）。
2. 表达式 (5<2)*120的值是（ 0 ）。
3. 表达式 (5>2)*(6~=5)的值是（ 1 ）。
4. 表达式 char(65)=='A' 的值是（ 1 ）。
5. 表达式 char(65)+1 的值是（66 ）。
6. 表达式 'A'+1的值是（ 66 ）。
7. 表达式 'A'+'B' 的值是（131 ）。
8. 存储double型数据占用内存（8 ）字节。
9. 存储single型数据占用内存（ 4 ）字节。
10. 清除命令窗口内容的命令是（ clc ）。
11. 删除工作空间中保存的变量x的命令是（clearx ）。
12. 将双精度实数的显示格式设置成15位定点小数格式的命令是（ format long ）。
13. 将横坐标轴标签设置成“时间（秒）”的语句是（xlabel('时间（秒）') ）。
14. 设置图例的Matlab库函数名是（ legend ）。
15. 绘制三维线图的Matlab库函数名是（ plot3 ）。
二、选择题（30分）1. 执行语句x=55后，Matlab将创建变量x，为其分配的存储空间的大小为（C）A)2字节 B)4字节 C)8字节 D)16字节2. 执行语句y=66后，Matlab将创建变量y，其类型为（D ）A)int8 B)int16 C)single D)double3. 下列整数类型中，不能参与任何运算的类型为（ D ）A)int8 B)int16 C)int32 D)int644. 设已执行语句x=3>2; y=x>0后，下面表达式中错误的是（ D ）A)x+y B)x-y C)x*y D)x/y5. 下列的数组写法中错误的是（C）A)[1:9] B)1:9 C)[1:2:9;2:2:8] D)[1:3;4:6;7:9]6. 设有数组定义：x=[1,2,3,4,5,6], y=x' ，下列表达式中正确的是（ D）A)y+x B)y-x C)y./x B)y*x7. 执行语句for x=1:2:10, disp(x), end，循环体将执行几次（ B）A)10次 B)5次 C)1次 D)0次8. 函数首部格式为function [out1,out2]=myfunc(in1,in2)，不正确的调用格式是（C ）A)[x,y]=myfunc() B)myfunc(a,b) C)[x,y]=myfunc(a)D)x=myfunc(a,b)9. 语句 x=-1:0.1:1;plot([x+i*exp(-x.^2);x+i*exp(-2*x.^2);x+i*exp(-4*x.^2)]' )，绘制（B ）A)1条曲线 B)3条曲线 C)21条曲线 D)0条曲线10. 下列哪条指令是求矩阵的行列式的值 （ C ）A) inv B) diag C) detD) eig三、解答题（40分）1．已知多项式323)(2345+++-=x x x x x f ，1331)(23--+=x x x x g ，写出计算下列问题的MATLAB 命令序列（1）)(x f 的根解：>> p1=[3,-1,2,1,3];>> x=roots(p1)x =0.6833 + 0.9251i0.6833 - 0.9251i-0.5166 + 0.6994i-0.5166 - 0.6994i（2）)(x g 在闭区间[-1,2]上的最小值解：>> [y,min]=fminbnd(@(x)((1/3)*x.^3+x.^2-3*x-1),-1,2)y =1.0000min =-2.66672.已知 ax -ax e -ex +ay =sin(x +a)+a ln 22, 写出完成下列任务的MATLAB 命令序列。
Matlab考试题库+答案
Matlab考试题库+答案填空题1、标点符号; %—用来表示该行为注释行。
可以使命令行不显示运算结果，2、x为0～4pi，步长为0.1pi的向量，使用命令 x=0:0.1*pi:4*pi 创建。
3、输入矩阵A=，使用全下标方式用A(2,2)取出元素“-5”，使用单下标方式用A(5)取出元素“-5”。
4、符号表达式sin(2*a+t)+m中独立的符号变量为t。
5、M脚本文件和M函数文件的主要区别是M脚本文件没有函数定义和M函数文件有函数定义_______。
6. 设x是一维数组，x的倒数第3个元素表示为x(_end-2_)设y为二维数组，要删除y的第34行和48列，可使用命令y(34,:)=[] ；y(:,48)=[]；7. 将变量x以Ascii文本格式存储到文件fname.txt，应使用命令save _x ；8. 在while 表达式, 语句体, End 循环语句中，表达式的值非零时表示循环条件为真，语句体将被执行，否则跳出该循环语句；9.要从键盘读入一个字符串并赋值给变量x，且给出提示“Who is she?”，应使用命令x=input(…Who is she??,?s?)；10．设A=和B=和C=均为m*n矩阵，且存在于WorkSpace中，要产生矩阵D=，可用命令D=(A-C)/B.^C，计算可用命令det(inv(A?*B)11. 在MATLAB命令窗口中的“>>”标志为MATLAB的命令行提示符，“│”标志为输入提示符。
12.已知A=[1 2 3；4 5 0；7 8 9]；B=[1 0 3；1 5 0；0 1 2]；写出下列各指令运行的结果。
A+B；A.*B；A==B ；ans= [2,2,6;5,10,0;7,9,11]ans= [1,0,9;4,25,0;0,8,18]ans= [1,0,1;0,1,1;0,0,0]13.已知A是矩阵，求A的对角矩阵函数是 diag ；求A的下三角矩阵函数是 tril。
Matlab期末考试题库(共12套卷)(word文档良心出品)
第一套华东交通大学2015-2016学年第一学期Matlab期末考试一、填空题(4*6=24分)1、在MATLAB^令窗口中的“ >>”标志为MATLAB勺_提示符，标志为提示符。
2、符号表达式sin(2*a+t)+m 中独立的符号变量为__ 。
3、在通常情况下，左除x=a\b是—的解，右除x=b/a是—的解，一般情况下，_______________ 。
4、为了使两个plot的图形在同一个坐标显示，可以使用 _______ 令进行图形保持；可以使用—命令为图形添加网格。
5、倘若要是对x进行赋值，从5到25,间隔是0.1,为_________ ;倘若要是对x进行赋值，从2到20，中间间隔100个点，为 ______6、A=[1,2;3,4];B=[1,0;0,1]; A*B= _ , A.*B=_ _ ;二、选择题(4*6=24分)7、如果x=1: 2 : 8,则x(1)和x⑷分别是()A. 1，8B. 1, 7C. 2, 8D. 2, 78、运行如下程序后，输入9回车，命令窗口( comma nd win dows)显示的结果为()c=input('请输入一个字符','s');if c>='A' & c<='Z'disp(setstr(abs(c)+abs('a')-abs('A')));elseif c>='a'& c<='z'disp(setstr(abs(c)- abs('a')+abs('A')));elseif c>='0' & c<='9'disp(abs(c)-abs('0'));else disp(c); endA. 9B. 8C. 7D. 69、MATLA表达式2*2A3A2的结果是()A. 128B. 4096C. 262144D. 25610、在循环结构中跳出循环，执行循环后面代码的命令为()(A) return (B) break (C) con ti nue (D) keyboard11、在图形指定位置加标注命令是()A. title(x,y, ^y=sin(x) ?B. xlabel(x,y, ?y=s in(x)C. text(x,y, ?y=sin(x) ?);D. lege nd(x,y, ?y=s in(x)12、下列哪个变量的定义是不合法的()。
MATLAB期末考试试卷及答案
页 01 共 页 7 第
__]6;5;4;3;2[__ = f ']6:2[=f�1 .stnemetats BALTAM gniwollof eht fo stluser eht etirW
。数函名匿和数函有私、 。 。 。容内的中口窗令命除清令命 、 、 和 、数函套嵌有型类数函的中 BALTAM .4 有件文 M 的中 BALTAM 3 、符算运学数括包符算运的中 BALTAM .2 用可中口窗令命 BALTAM .1
。 �
。果结的行运令指各列下出写;]0 5 1;3 4 1[=B;]0 3 1;6 9 0[=A 知已 .1 �分 03 共总�分 2 空每�题空填、一
。上卷题答在写笔用部全)等案答�行令指(答解 。料资和书的生考他其用借许 允不也�论讨许允不但�料资他其和书的用备己自用利可�卷开�同不卷试生考座领 .2 �网上许允不也�盘动移用许允不后考开 ,盘硬到拷序程备准将可前考开 .1 �明说试考机上
))01(pxe(gol �5 _]1 1;1 1;1 1[_ =sna
)2,3(seno �4 _]1 0 0;3 2 1;0 0 1[_ = a ;b=):,2(a;]3 2 1[=b;)3,3(eye=a�3 __]11 11;01 01;9 9;8 8;7 7[__ = n ]’m ’m[=n ;11:7=m �2
为值的 yarra 后行执序程的面下 .8 ))x(dar2ged(nis )D( )x(nis )C( )x(NIS )B( �)x(dar2ged�NIS )A(
为算运的数函弦正其算计� � 06 54 03� � x 度角 .7
页 01 共 页 4 第
。身本数该于等和之方立字数位各其�数位三个一指是� ”数花仙水“谓所。数花仙水的有所出印打 �分 51�程编行进言语 BALTAM 用使、四 号星为点记标�线虚色蓝为 2y�圈
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
安徽农业大学2006―2007学年第二学期 《MATLAB 语言程序设计》试卷（A 卷）
考试形式: 闭卷笔试，2小时
适用专业： 05级通信工程、电子信息工程
一、填空题(本题满分30分，每空3分)
1． 设有程序
A=[1,2,3,4;4,3,2,1;1, -2 ,1, -2];B=[3,1,-1;0,4,2];A1=A(:,[1 3]); M=size(A1)+ length(B(:,2))
将下列命令的运行结果填在横线上
M= ;A1+B' = .
2. A=[1,2,3,4;2,3,4,5;3,4,5,6];
A1=sum(A<4,2);A([2,3],:)=[],A2=A
A1=
; A2= ;
3.P=[1,2,3,4;3,0,1,2];max(P)= ;mean(P,2)=
院： 专业班级： 姓名： 学号：
装 订 线
二、(本题满分12分)试编写计算程序
V;
与特征向量
U
的特征值
B
2
A
求
(4).
E;
5
B
A
2B
BX
解矩阵
.
(3)
);
A
B)
(B(A
R
R
的秩
A
(2).
|;
B)
2
B)(A
(A
|
D
(1).
,
6
2
1
5
7
2
,
2
1
5
1
2
3
3
2
1
1
1
+
+
=
+
=
+
-
=
⎪
⎪
⎪
⎭
⎫
⎝
⎛
-
=
⎪
⎪
⎪
⎭
⎫
⎝
⎛
=
-
-
T
T
T
B
A
方程
设
三、(本题满分8分)
的程序组写出求解超定线性方程⎪⎪⎩⎪⎪⎨
⎧=++-=++=--=++1
21212:321
321321
321x x x x x x x x x x x x
四、(本题满分16分)运用符号运算功能写出下列各题Matlab 程序
解求方程03.12=-x e x
Adx dx A d x wx x x xe A x ⎰⎥⎦
⎤⎢⎣⎡=-1
022,,cos )sin(sin .2并计算
生成符号矩阵
3. 解求微分方程t x x x sin 23=+'-''
4. ).()(,)(2s F Laplace t f te t f t 变换函数的求设-=
五、（本题满分12分）
下列这组数据是美国1900—2000年人口的近似值（单位：百万）。
（1） 若.bt ae y t y =的经验公式为与试编写程序计算出上式中的a 、b ; （2） 画出数表中的散点图（红色圆圈o ）,以及拟合曲线图.bt
ae y = (蓝色实心线).
（3） 图形标注要求:无网格线,横标注“时间t ”，纵标注“人口数(百万)”，图
形标题“美国1900—2000年的人口数据”。
六、（本题满分12分）
编写程序解决下列问题：
（1）试产生80 4阶矩阵使其元素在0—100之间的整数并服从正态分布N (10，102 ),可将此矩阵理解为80名同学的四门课程(课程1, 课程2, 课程3, 课程4)
的成绩.
（2）计算出各门课程的平均成绩以及各名同学的平均成绩;
（3）统计出总的补考人次(60分以下为不及格需要补考);
（4）统计出四门课程成绩至少有三门课程小于60分的人数.
七、（本题满10分）
).
20,10(]4,0[)2())
,(:(;)1(,
12)12sin()(1
==--=∑
=n n x fourier y function k x
k x f n
k 分别取内图形用子图画出该函数在程序的第一句为返回该函数并保存试编写函数式文件程序设π
《MATLAB语言程序设计》参考程序
一、(30分)
1．M = 5 4
A1+B' =
4 3
5 6
0 3
2．
A1 =
3
2
1
A2 =
1 2 3 4
3．
max(P)=[ 3 2 3 4]
mean(P,2)=[2.5；1.5]
4．
y=20
5．
a =
0 1 2
0 2 4
0 0 6
6．
b =
90 80 70 65 80
7．228
二、(12分)
A=[1 2 3;3 2 1;5 1 2]
B=[2 0 0;7 5 0; 1 2 -6]
D=det((A'-B)*inv(A+2*B)) %(1)
R=rank(B*(A+B)* A') %(2)
X=pinv(B)*(2*B'*inv(A)*B+5*eye(3)) %(3)
[V,U]=eigs(A+2*B) %(4)
三、(8分)
A=[1,1 ,1;1,-1,-1;1,2,1;1,1,2];
b=[2,1,-1,1]'
x=pinv(A)*B
四、(15分)
%1
solve('exp(x)-3*x^2=0')
%2
syms x w
A=[x*exp(-x),x^2*sin(x);sin(w*x),cos(x)]
diff(A,x,2)
int(A,x,0,1)
%3
syms t x
dsolve('D2x-3*Dx+2*x=sin(t) ')
%4
syms t
f= t*exp(-2*t)
laplace (f, 's')
五、(10分)
t=[1900:10:2000];
y=[76 92106123132151179203227250281]; y1=log(y);
A=[ones(11,1),t']
B=A\y1';
a=exp(B(1))
b=B(2)
ye=a*exp(b*t);
plot(t,y,'ro', t, ye, 'b-')
grid off
xlabel('时间t')
ylabel('人口数(百万)')
title('美国1900—2000年的人口数据')
六、(12分)
%参考程序
cj=round(normrnd(70,10,80,4));
for i=1:80
for j=1:4
if cj(i,j)>100 cj(i,j)=100;
elseif cj(i,j)<0 cj(i,j)=0;
end
end
end
kmean=round(mean(cj))
disp('课程1 课程2 课程3 课程4')
disp(sprintf(' %4d%8d%8d%8d' ,kmean(1),kmean(2),kmean(3),kmean(4))) studentmean=round(mean(cj,2));
studentnumber=[1:80]';
[studentnumber, studentmean]
c=sum(sum(cj<60))
d=sum(sum(cj<60,2)>=3)
七、(10分)
(1).function y=fourier(x,n)
y=0;
for i=1:n
y=y+sin((2*i-1)*x)/(2*i-1);
end
y
(2).x=0:pi/20:4*pi;
y1=fourier(x,10);
y2= fourier(x,20);
subplot(211),plot(x,y1)
subplot(212),plot(x,y2)。