matlab练习

%把数字矩阵转化为单元矩阵
1 1 1 + +"+ ，当 n=100 时，求 y 的值. 3 5 2n − 1
程序如下： y=0;n=100; for i=1:100 y=y+1/(2*i-1); end y 29. while 语句 鸡兔同笼问题：计算笼子中动物的数量及脚的数量，求鸡兔各多少？（36,100） 程序如下： t=input('输入个数: '); j=input('输入脚数: '); i=1; while i if rem(j-i*2,4)==0& (i+(j-i*2)/4)==t break; % 终止循环 end i=i+1; end x=i, y=(j-i*2)/4 30. break 语句和 continuous 语句 求【100,200】之间第一个背 21 整除的整数。 for n=100:200 if rem(n,21)~=0 continue %跳过循环体中下面尚未执行的语句，接着执行下一次循环. end break end n 31. 函数文件： 另一种形式的 M 文件. function 输出形参表=函数名(输入形参表) 函数文件的建立：chap3_li6.m
2(1 + 0.34245 × 10−6 )
b+c
w=sqrt(2)*(1+0.34245e-6) 2) x =
2π a +
π + abc tan(b + c ) + a
Fra Baidu bibliotek− e2
, 其中a = 3.5, b = 5, c = −9.8
a=3.5;b=5;c=-9.8; x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a) 3) 2π a 2 ⎜ 1 − ⎢
⎡⎛ ⎣⎝
π 4
π ⎞ ⎛ ⎟ b − ⎜ 0.8333 − 4 ⎠ ⎝
⎞ ⎤ ⎟ a ⎥ , 其中a = 3.32, b = −7, 9 ⎠ ⎦
a=3.32; b=-7.9; y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a)
Matlab 程序设计
22. M 文件：命令文件和函数文件 命令文件可以直接运行，在 Matlab 命令窗口输入命令文件的名字，就会执行命令文件中的 命令。 函数文件不能直接运行，而要用函数调用的方式运行. 1) 建立命令文件 exp_1.m 3 中方式建立 M 文件 菜单操作，命令操作(edit), 工具栏操 作. 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') 建立命令文件 f2c.m clear; %清楚工作空间中的变量 f=input('Input Fahrenheit temperature:'); %稍后解释‘input’ c=5*(f-32)/9 保存为 f2c.m， 命令窗口中输入 f2c……….输入 73 命令 whos 查看工作空间中的变量，c，f 保留在工作空间中 建立函数文件 f2c1.m function c=f2c(f) c=5*(f-32)/9 在命令窗口中调用该函数： clear; y=input('Input Fahrenheit temperture:'); x=f2c1(y) 输入 70 whos c,f 未被保留在工作空间中，而 x，y 保留在工作空间中.
function [x,y]=chap3_li6(t,j) i=1; while i if rem(j-i*2,4)==0& (i+(j-i*2)/4)==t break; end i=i+1; end x=i;y=(j-i*2)/4;
运行结果为： >> [x,y]=chap3_li6(36,100) 练习 32. (P81) 从键盘输入一个 3 为整数，将它方向输出. 程序如下： t=input('请输入 3 位整数: '); a=fix(t/100); b=fix((t-a*100)/10); c=t-a*100-b*10; f=c*100+b*10+a; disp(f) 33. (P80)输入一个百分成绩，要求输出成绩等级 A,B,C,D,E. 其中 90~100 分为 A，80~89 分 为 B, 70~79 分为 C，60~69 分为 D，60 分以下为 E. 要求： (1): 分别用 if 语句和 switch 语句实现. (2): 输入百分制成绩后判断该成绩的合理性，对不合理的成绩输出出错信息. if 语句 t=input('请输入一个百分制成绩：'); if t<=100&t>=90 s='A'; elseif t<90&t>=80 %容易出错 s='B'; elseif t<80&t>=70 s='C'; elseif t<70&t>=60 s='D'; else t<60 s='E'; end disp(s) swith 语句 t=input('请输入一个百分制成绩：'); switch t case num2cell(90:100) s='A'; case num2cell(80:89) s='B'; case num2cell(70:79) s='C';
初识 Matlab:
1. x=[0:0.5:360]*pi/180; plot(x,sin(x),x,cos(x)); 2. p=[3,7,9,0,-23]; x=roots(p) 3. B=[1,2,3; 4,5,6;7,8,9] 4. a=2.7; b=13/25; c=[1,2*a+i*b, b*sqrt(a);sin(pi/4), a+5*b, 3.5+1] 5. x=linspace(1,5,10) 6. x=0:pi/10:22*pi; y=sin(x) %在工作空间中选择 y，绘制 y 的图形. 7. x=1+2i；y=3-sqrt(17); z=cos(abs(x+y)) 8. save me x y >> load me >> x+y 9. 建立M文件：mymatrix.m mymat=[101,102,103,104;105,106,107,108] x=linspace(1,5,10) mymat([1,2], end) mymat([1],2:end) mymat(:, 1:3:5) mymat(:, 1:3:4) mymat(:, [1,3,4]) 10. zeros(3) 11. A=[1,2,3;4,5,6]; 12. zeros(size(A)) 13. x=20+(50-20)*rand(5) 14. y=0.6+sqrt(0.1)*randn(5) 15. M=100+magic(5) %魔方矩阵 16. A=vander([2;4;6]) %范德蒙矩阵 17. x=0.1:0.3:1; y=sin(x).*cos(x) 18. A=fix(10+(90-10)*rand(5)) 19. A=[1,2,3;1,4,9;1,8,27]; b=[5,-2,6]'; % b 转置 x=inv(A)*b det(A) rank(A) trace(A) [V,D]=eig(A) % D 特征值构成的对角阵；V 相应特征向量构成的矩阵. 20. xm='University of Science and Technology Beijing' 21. 输出表达式的值：1) w =
⎧ x+ π , ⎪ ⎪ e2 ⎨ ⎪ 1 ln( x + 1 + x 2 ), ⎪ ⎩2
x≤0 x>0
程序如下： x=input('请输入 x 的值: ') if x<=0 y=(x+sqrt(pi))/exp(2); else y=log(x+sqrt(1+x*x))/2; y 26 输入一个字符，若为大写字母，则输出对应的小写字母；若为小写字母，则输出对应的 大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出. 程序如下 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); end 27. switch 语句 某商场对顾客所购买的商品实行打折销售，标准如下 price<200 没有折扣 200<=price<500 3%折扣 500<=price<1000 5%折扣 1000<=price<2500 8%折扣 2500<=price<5000 10%折扣 5000<price 14%折扣 输入所售商品的价格，求其实际销售价格. 程序如下： price=input('请输入商品的价格:'); switch fix(price/100) case {0,1} rate=0;
case num2cell(60:69) s='D'; otherwise s='E'; end disp(s) 判断成绩的合理性：不会 34.(P80)请输入 20 个数，求其中最大数和最小数。要求分别用循环结构和调用 Matlab 的 Max 函数、min 函数来实现.（不会） 程序如下： a=input('请输入 20 个数的向量:'); t=a(1,1);s=a(1,1); for i=1:20 if t>=a(1,i); t=a(1,i); end end for i=1:20 if s<=a(1,i); s=a(1,i); end end t s 35. (P81)编写程序，定义向量 N=[1,2,3,4,5], 然后利用向量 N 产生下列向量： (1) 2,4,6,8,10 (2) 1/2,1,3/2,2,5/2 (3) 1,1/2,1/3,1/4,1/5 (4) 1,1/4,1/9,1/16,1/25 程序如下 (2) (3) (4) (1) N=[1,2,3,4,5]; syms b N=[1,2,3,4,5]; syms b N=[1,2,3,4,5]; syms b N=[1,2,3,4,5]; for i=1:5 for i=1:5 for i=1:5 for i=1:5 b(1,i)=1/N(1,i); b(1,i)=1/(N(1,i)^2); b(1,i)=2*N(1,i); b(1,i)=N(1,i)/2; end end end end disp(b) disp(b) disp(b) disp(b) 36 (P81)编写程序，产生 20 个两位随机整数，输出其中小于平均值得偶数. 程序如下; A=fix(89*rand(1,20)+10) t=0; for i=1:20 t=t+A(1,i); end
程序设计
23. (数据的输入，输出). 输入 x，y 的值，并将它们的值互换后输出 程序如下： x=input('Input x please.'); y=input('Input y please.'); z=x;x=y;y=z;disp(x);disp(y) 输入 x 值：[12,3,3;43,6,3]，y 值：1:8 24. 求一元二次方程 ax 2 + bx + c = 0 的根.
a=input('a=?'); b=input('b=?'); c=input('c=?');d=b*b-4*a*c; x=[(-b-sqrt(d))/(2*a), (-b+sqrt(d))/(2*a)]; disp(['x1=',num2str(x(1)), ',x2=',num2str(x(2))]) 25.if 语言
case{2,3,4} rate=3/100; case num2cell(5:9) rate=5/100; case num2cell(10:24) rate=8/100; case num2cell(25:49) rate=10/100; otherwise rate=14/100; end price=price*(1-rate) 28. 已知 y = 1 +