matlab习题解答
更新时间:2024-01-20 09:53:01 阅读量: 教育文库 文档下载
- mATLAB运动仿真推荐度:
- 相关推荐
上机练习题一
班级: 姓名: 学号:
1.建立起始值=3,增量值=5.5,终止值=44的一维数组x 答案: x=(3:5.5:44)
o
2.写出计算 Sin(30)的程序语句. 答案: sin(pi*30/180) 或 sin(pi/6)
3.矩阵
?323??A??426????781??,矩阵
?111??B??222????333??;分别求出A?B及A与
B中对应元素之间的乘积.
答案:A = [3,2,3; 4,2,6; 7,8,1]
B = [1,1,1; 2,2,2; 3,3,3] A*B;A.*B
4计算行列式的值A?432326。答案:det(A)
7815对矩阵
?323??A??426????781??进行下述操作。
(1)求秩。答案:rank(A)
(2)求转置。答案:A' (3) 对矩阵求逆,求伪逆。答案:inv(A) ,pinv(A) (4) 左右反转,上下反转。答案:fliplr(A),flipud(A)
(5) 求矩阵的特征值. 答案:[u,v]=eig(A) (6) 取出上三角和下三角. 答案:triu(a) tril(a)
(7)以A为分块作一个3行2列的分块矩阵。答案:repmat(a)
6 计算矩阵
?535??374?????798??与
?242??679?????836??之和。
>> a=[5 3 5;3 7 4;7 9 8];
>> b=[2 4 2;6 7 9;8 3 6]; >> a+b 7 计算a????6293?75??与b????2441?68??的数组乘积。
>> a=[6 9 3;2 7 5];
>> b=[2 4 1;4 6 8]; >> a.*b ans =
12 36 3
8 42 40 8 已知:
?123??a???456???789??,分别计算a的数组平方和矩阵平
方,并观察其结果。
>> a=[1 2 3;4 5 6;7 8 9]; >> a.^2 ans =
1 4 9 16 25 36 49 64 81 >> a^2 ans =
30 36 42 66 81 96 102 126 150
上机练习题二
班级: 姓名: 学号:
1 对于AX?B,如果
?492??A???764???357??,
?37??B???26???28??,求解X。
>> A=[4 9 2;7 6 4;3 5 7];
>> B=[37 26 28]’; >> X=A\\B X =
-0.5118 4.0427 1.3318
2 角度x??304560?,求x的正弦、余弦、正切和余切。 >> x=[30 45 60]; >> x1=x/180*pi; >> sin(x1) ans =
0.5000 0.7071 0.8660 >> cos(x1) ans =
0.8660 0.7071 0.5000 >> tan(x1) ans =
0.5774 1.0000 1.7321 >> cot(x1) ans =
1.7321 1.0000 0.5774
43 将矩阵a???5?2?7??1??59?c?、b????7和?62?组合成两个新矩阵: 83????(1)组合成一个4?3的矩阵,第一列为按列顺序排列的a矩阵元素,第二列为按列顺序排列的b矩阵元素,第三列为按列顺序排列的c矩阵元素,即 ?75?
?4?586???219??
?732?? ( 2 )按照 a、b、c的列顺序组合成一个行矢量,即?452778135692? 答案:a=[4 2;5 7]; b=[7 1;8 3]; c=[5 9;6 2]; % (1)
>> d=[a(:) b(:) c(:)] d =
4 7 5 5 8 6 2 1 9 7 3 2 % (2)
>> e=[a(:);b(:);c(:)]' e =
4 5 2 7 7 8 1 3 6 9 2
或利用(1)中产生的d >> e=reshape(d,1,12)
5
ans =
4 5 2 7 7 8 1 3 5 6 9 2
4求解在x=8时多项式(x-1)(x-2) (x-3)(x-4)的值。 >> p=poly([1 2 3 4]); >> polyvalm(p,8) ans = 840
5求方程3x?7x?9x?23?0的全部根。
p=[3,7,9,0,-23]; %建立多项式系数向量 x=roots(p) %求根
432
上机练习题三
班级: 姓名: 学号:
1、 设x是数组,求均值和方差 解:函数文件如下:
function [xx,s]=func1(x) n=length(x); xx=sum(x)/n;
s=sqrt((sum(x.^2)-n*xx^2)/(n-1)); 命令窗口:
>> x=[1 2 3 4 5];[xx,s]=func1(x) 2、求满足?mln(1?n)?100的最小m值
n?0s=0; n=0;
while(s<=100) s=s+log(1+n); n=n+1; end n,s
3、用循环语句形成
Fibonacci
F1?F2?1,F1?5k?Fk?1?Fk?2,k?3,4,....。并验证极限FkF? (提示:
k?12 计算至两边误差小于精度1e-8为止) 解: 求Fibonacci数列的函数文件: function f=fun(n)
if n<=2 f=1; else
f=fun(n-1)+fun(n-2); end
数列
验证极限的函数文件:
function [k,a]=funTest(e) a=abs(1-(1+sqrt(5))/2); k=2;
while(a>e) k=k+1;
a=abs(fun(k)/fun(k-1)-(1+sqrt(5))/2); end 命令行:
>> [k,a]=funTest(10^-8) k = 21 a =
9.7719e-009 或者M文件如下:
clear; F(1)=1;F(2)=1;k=2;x=0; e=1e-8; a=(1+sqrt(5))/2; while abs(x-a)>e k=k+1; F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1); end a,x,k 4、分别用for和while循环结构编写程序,求出K??23,
106i?1i并考虑一种避免循环语句的程序设计,比较各种算法的运行时间。
解:for循环结构:M文件loop.m
k=0;
for i=1:10^6
k=k+sqrt(3)*2^-i; end k
while循环结构:M文件loop1.m
k=0;i=1;
while i<=10^6
k=k+sqrt(3)*2^(-i); i=i+1; end k
非循环结构:M文件nonLoop.m i=1:10^6;
x=sqrt(3)*(2.^-i); k=sum(x) 速度比较:>>tic;loop;toc %循环结构的执行时间
k =
1.7321
Elapsed time is 1.813000 seconds.
>> tic;nonLoop;toc %非循环结构的执行时间
k =
1.7321
Elapsed time is 1.094000 seconds.
上机练习题四
班级: 姓名: 学号:
1、作图描述气温变化 >> x=0:24;
>>
y=[15,14,14,14,14,15,16,18,20,22,23,25,28,31,32,31,29,27,25,24,22,20,18,17,16]; >> plot(x,y)
2、作出下列函数图形
(1)y?xsin(x?x?2) ?2?x?2 (分别使用plot和fplot完成) 解:>> fplot('x^2*sin(x^2-x-2)',[-2 2]) %fplot方法 >>
x=-2:0.1:2;y=x.^2.*sin(x.^2-x-2);plot(x,y) %plot方法 如图(4.1)
(2)x4?y9?1 (椭圆 提示:用参数方程)
2222解:>> r=-pi:0.1:pi;x=2*cos(r);y=3*sin(r);plot(x,y) % 如图(4.2) 解法二
x=-2:1/100:2;
y1=3*sqrt(1-x.^2/4); y2=-3*sqrt(1-x.^2/4);
plot(x,y1,'r-',x,y2,'r-'); axis equal tight;
图(4.1) 图(4.2)
(3)z?x?y (抛物面) x?3,y?3 解:(错误)>> x=[-3:0.1:3];y=[-3:0.1:3];z=x.^2+y.^2; plot3(x,y,z) % 如图(4.31) (正确)>> xa=-3:0.1:3;ya=-3:0.1:3;[x,y]=meshgrid(xa,ya); % 如图(4.32)
>> z=x.^2+y.^2;mesh(x,y,z); >> surf(x,y,z)
22
图(4.31)error 图(4.32)
(4)曲面z?x?3x?y?2x?2y?2xy?6,x?3,?3?y?13
解: >> xa=linspace(-3,3,100);ya=linspace(-3,13,100); >> [x,y]=meshgrid(xa,ya);
>> z=x.^4+3*x.^2+y.^2-2*x-2*y-2*x.^2.*y+6; >> mesh(x,y,z)
4222>> surf(x,y,z)
(5)空间曲线x?sint,y?cost,z?cos(2t),0?t?2
解:>> t=linspace(0,2,50);x=sin(t);y=cos(t);z=cos(2*t); >> plot3(x,y,z)
(6)半球面x?2sin?cos?,y?2sin?sin?,z?2cos?,0???360,0???90 解: >>a=linspace(0,2*pi,50);b=linspace(0,pi/2,50); >> [a,b]=meshgrid(a,b);
>> x=2*sin(a).*cos(b);y=2*sin(a).*sin(b);z=2*cos(a); >> surf(x,y,z)
(7)三条曲线合成图y?sinx,y?sintsin(10x),y??sinx,0?x?? 解: >>x=linspace(0,pi,50);y1=sin(x); >> plot(x,y1);hold on; >> y2=sin(x).*sin(10*x); >> plot(x,y2); >> y3=-sin(x); >> plot(x,y3); >> hold off;
??123
3、作下列分段函数图
?1.1x?1.1?y??xx?1.1??1.1x??1.1?
x=-5:0.1:5;
for i=1:length(x) if x(i)>1.1 y(i)=1.1; elseif x(i)<-1.1 y(i)=-1.1; else
y(i)=x(i); end end
plot(x,y); grid on;
4、用MATLAB函数表示下列函数,并作图。
?0.5457exp(?0.75y2?3.75x2?1.5x),x?y?1?p(x,y)??0.7575exp(?y2?6x2),?1?x?y?1?0.5457exp(?0.75y2?3.75x2?1.5x),x?y??1?
解:建立M文件pxy如下:
xa=-2:0.05:2;ya=xa;
nx=length(xa);ny=length(ya); [x,y]=meshgrid(xa,ya); z=zeros(nx,ny);
[a1,b1]=find(x+y>1); %第a1列b1行对应的
x+y>1 %值是第a1 列对应的(x对应列;x值是y对应行xa(a1)) ;第
b1行对应的yz((a1-1)*ny+b1)=0.5457*exp(-0.75*ya(b1).^2-3.75*ya(b1)
xa(a1).^2-1.5*xa(a1)); [a2,b2]=find(x+y<=1&x+y>-1);
z((a2-1)*ny+b2)=0.7575*exp(-ya(b2).^2-6*xa(a2).^2);
[a3,b3]=find(x+y<=-1);
z((a3-1)*ny+b3)=0.5457*exp(-0.75*ya(b3).^2-3.75*xa(a3).^2+1.5*xa(a3)); 命令窗口:surf(x,y,z); 运行结果如右图:>> pxy
或者 M 文件如下: (xa,ya);
clear;close;
xa=-2:0.1:2;ya=-2:0.1:2;[x,y]=meshgrid (k1).^2-1.5*x(k1));
z=zeros(size(x)); k1=find(x+y>1);
z(k1)=0.5457*exp(-0.75*y(k1).^2-3.75*x ;
k2=find(x+y<=1&x+y>-1);
z(k2)=0.7575*exp(-y(k2).^2-6*x(k2).^2) k3=find(x+y<-1);
(k3).^2+1.5*x(k3)); z(k3)=0.5457*exp(-0.75*y(k3).^2-3.75*x mesh(x,y,z);
上机练习题五
号:班级: 姓名:
学1解:、运行>>demo
demo 2、查询trapz的功能、用法、目录、程序结构、相同目录下其它文件
解: >> help trapz ――功能用法
>> type trapz――程序结构,源码 >> which trapz――所在目录
>> help C:\\MATLAB6p5\\toolbox\\matlab\\datafun――该目录下其它文件
3在[0,4pi]画sin(x),cos(x)(在同一个图象中); 其中cos(x)图象用红色小圆圈画.并在函数图上标注 “y=sin(x)”, “y=cos(x)” ,x轴,y轴,标题为“正弦余弦函数图象”. x=linspace(0,4*pi,100); y=sin(x); plot(x,y);
gtext('y = sin(x)'); % 图形注解,注意要用鼠标定位 hold on; y=cos(x);
plot(x,y,'ro');
gtext ('y = cos(x)'); % 图形注解 xlabel('x轴'); % x轴注解 ylabel('y轴'); % y轴注解
title('正弦余弦函数图象'); % 图形标题
4从键盘输入若干个数,当输入0时结束输入,求这些数的平均值和它们之和。 sum=0;
cnt=0;
val=input('Enter a number (end in 0):'); while (val~=0) sum=sum+val; cnt=cnt+1;
val=input('Enter a number (end in 0):'); end
if (cnt > 0) sum
mean=sum/cnt end
5若一个数等于它的各个真因子之和,则称该数为完数,如6=1+2+3,所以6是完数。求[1,500]之间的全部完数。 for m=1:500 s=0;
for k=1:m/2 if rem(m,k)==0 s=s+k; end end
if m==s
disp(m);
end end
上机练习题六
班级: 姓名: 学号:
1 假定数据点来源为:f(x)?(x?3x?5)esinx,试根据生成的数据进行插值处理,得出较平滑的曲线。
x=0:.12:1;y=(x.^2-3*x+5).*exp(-5*x).*sin(x);plot(x,y,x,y,'o')
x1=0:.02:1;y0=(x1.^2-3*x1+5).*exp(-5*x1).*sin(x1); y1=interp1(x,y,x1);
>> y2=interp1(x,y,x1,'cubic'); >> y3=interp1(x,y,x1,'spline'); >> y4=interp1(x,y,x1,'nearest'); >> plot(x1,[y1;y2;y3;y4],':',x,y,'o',x1,y0)或者plot(x1,[y1',y2',y3',y4'],':',x,y,'o',x1,y0)
2 用不同插值的方法计算sin(x)在pi/2的值 clear clc
x=0:0.2:pi;%产生包含被插值点的采样点(做成一个向量)
y=sin(x);%求出各采样点对应的样本值
y1=interp1(x,y,pi/2);%用默认的'linear'方法计算sin(pi/2) y2=interp1(x,y,pi/2,'nearest');%用默认的'linear'方法计
2?5x算sin(pi/2)
y3=interp1(x,y,pi/2,'cubic');%用三次多项式插值方法计算sin(pi/2)
y4=interp1(x,y,pi/2,'spline');%用三次样条插值方法计算sin(pi/2)
y5=spline(x,y,pi/2);%直接用'spline'方法计算sin(pi/2),功能与y4相同
disp('各种方法的插值结果:') out=['y1=', num2str(y1) 'y2=',num2str(y2); 'y3=',num2str(y3)] out2=['y4=',num2str(y4) 'y5=',num2str(y5)] 3求
?dx?dt?4x?2y?dy??2x?y?dt 的通解
[x,y]=dsolve('Dx=4*x-2*y','Dy=2*x-y', 't') %P149(3)的求解
[x,y]=dsolve('Dx=a*y','Dy=-b*x ', 'x(0)=x0', 'y(0)=y0', 't')
y=dsolve('Dy=(b*x)/(a*y)', 'y(0)=y0','x'); 4
?'y2?t?2?y?, 0?t?1求?的数值解,并与精确解y(t)?t?1?1相比较。4(t?1)?y(0)?2?
方法一:
%微分方程的M函数funt.m文件 function y=funt(t,y) y=(y^2-t-2)/4/(t+1);
%以下是求解的脚本m函数,可自由取名,然后在command窗口调用该函数求解. ts=[0,10];%自变量的求解区间 y0=2;%初值条件
[t,y]=ode23('funt',ts,y0);%用2,3阶龙格库塔方法求‘funt’文件里的微分方程
y1=sqrt(t+1)+1;%求精确解对应点上的函数值
[t,y,y1]%以三列的形式显示结果,其中第一列为采样点t对应的值,第二列为t对应的数值解,第二列为t对应的精确解。
方法二:
%或者直接用内建函数编写待求的微分方程,在一个M脚本文件里执行求解。
funt=inline('(y^2-t-2)/4/(t+1)') %内建函数编写待求的微分方程
ts=[0,10];%自变量的求解区间 y0=2;%初值条件
[t,y]=ode23(funt,ts,y0);%用2,3阶龙格库塔方法求‘funt’文件里的微分方程
y1=sqrt(t+1)+1;%求精确解对应点上的函数值 [t,y,y1] 5
lorenz模型的状态方程表示为?1(t)???x1(t)?x2(t)x3(t)x?2(t)???x2(t)??x3(t)x?3(t)??x2(t)x1(t)??x2(t)?x3(t)x8取??10,??28,??,求解微分方程组,并绘制系统相平面图象。3
第一种方法建立lorenz函数模型的状态方程
function xp=lorenz(t,x)%建立lorenz函数模型的状态方程
xp=[-8/3,0,x(2);0,-10,10;-x(2),28,-1]*x%表明x是一个三维向量,前面是一个线性方程组的系数矩阵,xp%是一个三维的输出表示x(1),x(2),x(3)的一阶导数向量。
第二种方法建立lorenz函数模型的状态方程
function xp=lorenz1(t,x)%建立lorenz函数模型的状态方程
xp=[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(2)*x(1)+2
8*x(2)-x(3)];
%求解微分方程 clear;clc;
x0=[0,0,eps]';%三个初值条件构成的向量 [t,x]=ode23('lorenz',[0,100],x0); [t,x]
plot(t,x),grid,pause%绘制解x(1),x(2),x(3)各自相对于变量t的图象,按任意键后继续下面的程序
figure;plot3(x(:,1),x(:,2),x(:,3));%绘制解x(1),x(2),x(3)的关于系统的相平面图象 axis([10,40,-20,20,-20,20]);
上机练习题七
班级: 姓名: 学号:
2
1有一组测量数据如下表所示,数据具有y=x的变化趋势,用最小二乘法求解y。
1.2.3.4.x 1 2 3 4 5 5 5 5 5 -12.5.8.12161826y 3 .4 7 9 4 .2 .6 .8 .2
>> x=[1 1.5 2 2.5 3 3.5 4 4.5 5]'
>> y=[-1.4 2.7 3 5.9 8.4 12.2 16.6 18.8 26.2]' >> e=[ones(size(x)) x.^2]
>> c=e\\y
>> x1=[1:0.1:5]';
>> y1=[ones(size(x1)),x1.^2]*c; >> plot(x,y,'ro',x1,y1,'k')
2求下列线性方程组的解 (1)
??41?1??x1??9??32?6????x?2?????2?
??53???1???x3????1??解1];x=a\\b :>> a=[4 1 -1;3 2 -6;1 -5 3];b=[9; x =
%唯一解 2.0213 2.3830 1.4894 (2)
??4?33??x1??????1?32?6x????2
????2?????1?53???x3????1??解3];b=[-1;-2;1];x=a\\b :>> a=[4 2
-6;1
x =
-3 %唯一解3;3 -0.4706 -0.2941 0 (3)??41??1?32???x1????x???1? ??1?5??2?????1??解
:
>>
a=[4
1;3
-2; -5
2;1
-5];b=[1;1;1];x=a\\b %最小二乘近似解 x =
0.3311 -0.1219 (4)
?x1?21?11?????1??121?1??x2???2????x???3??1121???x???3???4?,求通解
解:>> a=[2 1 -1 1;1 2 1 -1; 1 1 2 1];b=[1;2;3];
>> rank(a),rank([a,b]) ans = 3 ans =
3 %说明有无穷多解 >> rref([a,b]) %行最简化 ans =
1.0000 0 0 1.5000 1.0000
0 1.0000 0 -1.5000 0
0 0 1.0000 0.5000 1.0000
%通解为:x??1.5x?1,x?1.5x,x??0.5x?1
3、求下列矩阵的行列式、逆、特征值和特征向量
142434(1)
?41?1???32?6???1?53???
-5
解:>> a=[4 1 -1;3 2 -6;1 3];det(a),inv(a),[v,d]=eig(a) ans =
-94 行列式 ans =
0.2553 -0.0213 0.0426
0.1596 -0.1383 -0.2234 矩阵的逆
0.1809 -0.2234 -0.0532 v =
0.0185 -0.9009 -0.3066
-0.7693 -0.1240 -0.7248 特征向量
-0.6386 -0.4158 0.6170 d =
-3.0527 0 0
0 3.6760 0 特征值
0 0 8.3766 (2)
1??43???33?5??1?53???
解:>> a=[4 3 1;3 3 -5;1 -5 3];det(a),inv(a),[v,d]=eig(a) ans =
-124 行列式 ans =
0.1290 0.1129 0.1452
0.1129 -0.0887 -0.1855 矩阵的逆
0.1452 -0.1855 -0.0242 v =
0.3757 -0.8583 0.3496
-0.6881 -0.0056 0.7255 特征向量
-0.6208 -0.5131 -0.5927 d =
-3.1480 0 0
0 4.6176 0 特征值
0 0 8.5304 (3)
?5765????71087??68109????57910???
解:>> a=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10];det(a),inv(a),[v,d]=eig(a) ans =
1 行列式 ans =
68.0000 -41.0000 -17.0000 10.0000
-41.0000 25.0000 10.0000 -6.0000 矩阵的逆
-17.0000 10.0000 5.0000 -3.0000 10.0000 -6.0000 -3.0000 2.0000
v =
0.8304 0.0933 0.3963 0.3803
-0.5016 -0.3017 0.6149 0.5286 特征向量
-0.2086 0.7603 -0.2716 0.5520 0.1237 -0.5676 -0.6254 0.5209 d =
0.0102 0 0 0 0 0.8431 0 0
0 0 3.8581 特征值
0 0 0 30.2887 n=5;
tria = 6*eye(n-1); %构造上三角 rown = zeros(1,n-1); coln = zeros(n,1); tria = [tria;rown]; tria = [coln,tria]
eyeb = 5*eye(n); %构造对角阵
tric = 1*eye(n-1); %构造下三角 tric = [rown;tric];
0 (4) tric = [tric,coln]
A=tria+eyeb+tric
4求多项式的根,并分析误差大小 (1)x?x?1
解:>> A=[1,1,1];roots(A) ans =
-0.5000 + 0.8660i -0.5000 - 0.8660i
>> y=polyval(A,ans) %验证 y =
1.0e-015 * 0.3331 0.3331 (2)3x?4x?2x?1
解:>> A=[3 0 -4 0 2 -1];x=roots(A) x =
-0.9479 + 0.3845i -0.9479 - 0.3845i 1.0000 0.4479 + 0.3435i 0.4479 - 0.3435i
>> y=polyval(A,x) %验证 y =
1.0e-013 *
-0.0144 - 0.1138i -0.0144 + 0.1138i -0.0888 0.0044 - 0.0083i 0.0044 + 0.0083i (3)5x?6x?8x?5x
解: >> A=[5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6 8 0 0 0 -5 0 0];x=roots(A)
>> y=polyval(A,x) (4) (2x?3)?4 解:p1=[2 33];p2=conv(p1,p1); p=conv(p2,p1); %(2x+3)
n=length(p); for (i=1:n-1) q(i)=p(i); end
q(n)=p(n)-4; 3
%(2x+3)-4
x=roots(q)
253237623结果:x =
-1.8969 + 0.6874i -1.8969 - 0.6874i -0.7063
验证:>> y=polyval(q,x) y =
1.0e-014 *
-0.7105 - 0.6217i -0.7105 + 0.6217i 0
上机练习题八
班级: 姓名: 学号:
1有一正弦衰减数据y=sin(x).*exp(-x/10),其中x=0:pi/5:4*pi,用三次样条法进行插值。 >> x0=0:pi/5:4*pi;
>> y0=sin(x0).*exp(-x0/10); >> x=0:pi/20:4*pi; >> y=spline(x0,y0,x);
>> plot(x0,y0,'or',x,y,'b')
2 编制一个解数论问题的函数文件:取任意整数,若是偶数,则用2除,否则乘3加1,重复此过程,直到整数变为1。
function c=collatz(n) % collatz
% Classic “3n+1” Ploblem from number theory c=n;
while n>1
if rem(n,2)==0 n=n/2; else
n=3*n+1; end
c=[c n]; end
3有一组学生的考试成绩(见表),根据规定,成绩在100分时为满分,成绩在90~99之间时为优秀,成绩在80~89分之间时为良好,成绩在60~79分之间为及格,成绩在60分以下时为不及格,编制一个根据成绩划分等级的程序。
学生王 张 刘 李 陈 杨 于 黄 郭 赵 姓名 1成 7859896560绩 2 3 6 4 0 8 6 8 4 5
Name=['王','张','刘','李','陈','杨','于','黄','郭','赵']; Marks=[72,83,56,94,100,88,96,68,54,65];
% 划分区域:满分(100),优秀(90-99),良好(80-89),及格(60-79),不及格(<60)。
n=length(Marks); for i=1:n
a{i}=89+i; b{i}=79+i; c{i}=69+i; d{i}=59+i; end; c=[d,c];
% 根据学生的分数,求出相应的等级。 for i=1:n
switch Marks(i)
case 100 %得分为100时 Rank(i,:)=' 满分';
case a %得分在90~99之间
Rank(i,:)=' 优秀';
case b %得分在80~89之间 Rank(i,:)=' 良好';
case c %得分在60~79之间 Rank(i,:)=' 及格';
otherwise %得分低于60。 Rank(i,:)='不及格'; end end
% 将学生姓名,得分,级等信息打印出来。 disp(' ')
disp(['学生姓名 ',' 得分 ',' 等级']); disp('--------------------------') for i=1:10;
disp([' ',Name(i),' ',num2str(Marks(i)),' ',Rank(i,:)]); end
学生姓名 得分 等级 ----------------------------------------- 王 72 及格 张 83 良好 刘 56 不及格 李 94 优秀 陈 100 满分 杨 88 良好 于 96 优秀 黄 68 及格 郭 54 不及格 赵 65 及格
4 编写一段程序,能够把输入的摄氏温度转化成华氏温度,也能把华氏温度转换成摄氏温度。
k=input('选择转换方式(1--摄氏转换为华氏,2--华氏转换为摄氏):'); if k~=1 & k~=2
disp('请指定转换方式') break end
tin=input('输入待转变的温度(允许输入数组):');
if k==1
tout=tin*9/5+32; % 摄氏转换为华氏 k1=2; elseif k==2
tout=(tin-32)*5/9; % 华氏转换为摄氏 k1=1; end
str=[' ?C';' ?F'];
disp(['转换前的温度', ' ', '转换后的温度']) disp([' ',num2str(tin),str(k,:), ' ', num2str(tout),str(k1,:)])
正在阅读:
matlab习题解答01-20
情景交融高考诗歌鉴赏题05-26
机电一体化专业实践现状分析05-20
河南省郑州外国语学校2009-2010高二下期期中数学试题(文)(含答案04-24
维修电工中级理论知识试卷四(含答案)01-13
旅游路上作文400字07-07
塔吊四桩基础的计算书(3AQTZ50)01-25
中国共产党纪律处分条例测试题库05-19
超外差焊接收音机课设论文06-22
Westernblot实验操作步骤05-02
- exercise2
- 铅锌矿详查地质设计 - 图文
- 厨余垃圾、餐厨垃圾堆肥系统设计方案
- 陈明珠开题报告
- 化工原理精选例题
- 政府形象宣传册营销案例
- 小学一至三年级语文阅读专项练习题
- 2014.民诉 期末考试 复习题
- 巅峰智业 - 做好顶层设计对建设城市的重要意义
- (三起)冀教版三年级英语上册Unit4 Lesson24练习题及答案
- 2017年实心轮胎现状及发展趋势分析(目录)
- 基于GIS的农用地定级技术研究定稿
- 2017-2022年中国医疗保健市场调查与市场前景预测报告(目录) - 图文
- 作业
- OFDM技术仿真(MATLAB代码) - 图文
- Android工程师笔试题及答案
- 生命密码联合密码
- 空间地上权若干法律问题探究
- 江苏学业水平测试《机械基础》模拟试题
- 选课走班实施方案
- 习题
- 解答
- matlab
- 交换技术作业及演讲题目
- 浙江绍兴一中2014—2015学年高二下学期期中考试 英语(Word版含答案)
- 南京市险房翻建规划管理办法
- 机车联控用语
- 毕业论文(机械手)
- 2015监理工程师网络继续教育专业课试题(附答案)
- 总体开工报告
- 基于单片机的汽车超载控制系统的设计 - 图文
- 基于UG的齿轮泵三维设计与仿真运动 - 图文
- 中高端珠宝大型落地展会营销活动策划方案 - 图文
- 我县城乡居民社会养老保险试点简报四期7.16 - 图文
- 2016年下半年中小学教师资格考试 物理学科知识与教学能力试题(初级中学)试卷真题及答案
- 2013-2014北京初一数学十个区期末考试题(含标准答案)
- 有效积累小学数学活动经验的策略
- 中国人民政协会议的情况简介
- 180012007内审员试卷
- 初中生物花的结构和类型知识点及复习题
- 最新苏教版数学一(下)教案
- 小学语文教学中的阅读情感熏陶策略-精品教育文档
- 数据结构实验报告4(电大)