matlab编程题整理

作业一1输出x,y两个中值较大的一个值x=input(x);y=input(y);if xy xelse yend2输入x，计算y的值。计算函数的值y=x+1,x0,y=2x-1,x0x=input错误！未指定书签。(x);if x100 | x=90 y=优秀elseif x=80 y=良好elseif x=70 y=中等elseif x=60 y=及格else y=不合格emd4某超市节日期间举办购物打折的促销活动，优惠办法是：每位顾客当天一次性购物在100元以上者，按九五折优惠；在200元以上者，按九折优惠；在300元以上者，按八五折优惠；在500元以上者，按八折优惠。x=input(x);if x=500y=x*0.8elseif x=300y=x*0.85elseif x=200y=x*0.9elseif x=100y=x*0.95else y=xend5编程计算：s=1+2+3+100sum=0;for i=1:100 sum=sum+i;endsum引申 1!-2!+3!-4!+5!-+99!sum=0;for i=1:99 pdr=1; for k=1:i pdr=pdr*k; endsum=sum+pdr*(-1)(i-1);endsum引申 1*2*3*4*100 sum=1；for i=1:100 sum=sum*iendsum6计算1100的奇数和sum=0;for i=1:2:100 sum=sum+i;endsum7百元买百鸡问题。假定小鸡每只5角，公鸡每只2元，母鸡每只3元。编程列出所有可能的够鸡方案for x=1:33 for y=1:50for z=1:200if x+y+z=100 & 3*x+2*y+0.5*z=100 x,y,z end endendend8我国有13亿人口，按人口年增长0.8%计算，多少年后我国人口超过26亿n=1;m=13;while m=3），他的每个位上的数字的n次幂之和等于它本身。求出1000以内for i=100:999 a(1)=fix(i/100); a(2)=fix(i-a(1)*100)/10); a(3)=i-a(1)*100-a(2)*10;ifa(1)3+a(2)3+a(3)3=i i endend10鸡和兔子关在一个笼子里，已知共有头36个，脚100个，求笼内关了多少只兔子和多少只鸡？for i=1:50 for j=1:36ifi+j=36 & i*2+j*4=100 i,j endendend作业二1.用matlab编程求算式xyz+yzz=532中x， y， z的值（其中xyz和yzz分别表示一个三位数）。clearfor x=1:9 for y=1:9 for z=1:9i=100*x+10*y+z+100*y+10*z+z;if i=532 x,y,z end endendend2.用matlab编程完成用一元人民币换成一分、两分、五分的所有兑换方案（即输出所有的组合方式）。clearfor i=0:20 for j=0:50 k=100-i*5-j*2; if(k=0) I,j,k end endend3.有1020个西瓜，第一天卖一半多两个，以后每天卖剩下的一半多两个，问几天后可以卖完，请用matlab编程计算。day=0; x1=1020while (x1)x2=x1/2-2; x1=x2; day=day+1;endday4.有一堆零件（100到200之间），如果分成4个零件一组的，则多2个零件；若分成7个零件一组的，则多3个零件；若分成9个零件一组的，则多5个零件。用matlab编程求这堆零件总数。for x=100:200 if mod(x,4)=2 if mod(x,7)=3 if mod(x,9)=5 x end end endend5.编写程序，求1000至9999之间的回文数回文数是指正读与反读都一样的数，如1221.for i=1000:1000:9999 for j=0:9 n=i+j*100+j*10=i/1000 endend6.作函数y=sinx,x0,2的图象，用蓝色五角星表示点。 x=linspace(0,2*pi);y=sin(x);plot(x,y,bp);作业三1: 编写M函数 , 计算function y=fun1(x) y=(x3-2*x2+x-6.3)/(x2+0.05*x-3.14);fun1(1)*fun1(2)+fun1(3)*fun1(3) 2：编写M函数，计算，function y=fff(x)if x=2 & x8 & x=20y=4*x-5;elsey=cos(x)+sin(x);end3：根据pi*pi/6=1/12+1/22+1/32+1/n2,求pi的近似值。当n分别取100，1000，10000时，结果是多少？ y=0;n=100;for i=1:ny=y+1/i/i;endypi=sqrt(6*y)y=0;n=1000;for i=1:ny=y+1/i/i;endypi=sqrt(6*y)y=0;n=10000;for i=1:ny=y+1/i/i;endypi=sqrt(6*y)4:编程:求100，200之间第一个能被21整除的整数。for n=100:200if mod(i,21)=0continuebreakendend5：编写函数文件求半径为r的圆的面积和周长。function s,p=fcircle(r)s=pi*r*r;p=2*pi*r;6：根据y=1+,求：（1）y3时的最大n值。（2）与（1）的n值对应的y值。y=1;n=1;while(y0) x=(x+1)*2; day=day-1;endx作业五上机题1.给定一个矩阵a = 1 2 3 4 5 6编写一个M函数，要求输入是a，输出各元素之和。function y=f(a)r,q=size(a);y=0;for i=1:r for j=1:q y=y+a(i,j); endend2. A 是一个mn的矩阵. 写一M函数, 算出A中有多少个零元素。function y=f(a)m,n=size(a);y=0;for i=1:m for j=1:n if a(i,j)=0 y=y+1; end endend3. 矩阵乘法运算要求两矩阵的维数相容，否则会出错。先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。A=input(A);B= input(B);tryC=A*B;catchC=A.*B;endC4. 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。for i=100:999 a(1)=fix(i/100);a(2)=fix(i-a(1)*100)/10);a(3)=i-a(1)*100-a(2)*10;if a(1)3+a(2)3+a(3)3=i iendend5.fibonccci数组的元素满足fibonacci规则：ak+2=ak+ak+1,(k=1,2,);且a1=a2=1。现要求求出该数组中第一个大于10000的元素。a(1)=1;a(2)=1;i=1;while a(i)=10000a(i+2)=a(i)+a(i+1);i=i+1;enda(i)6.计算1！+2！+3！+4！+5！6！+n！值的程序。（n可以自由输入）Sum=0N=input(n)For i=1 to nSum1=1For j=1 to iSum1=sum1*iEndSum=sum+sum1end作业六1.求解一阶微分方程（1） y=dsolve(Dy=a*y+b,x)（2）f,g=dsolve(Df=f+g,Dg=g-f,Df(0)=1,Dg(0)=1,x)(3) 求解微分方程y=dsolve(Dy=-y+t+1,y(0)=1)(4) 求方程的通解.y=dsolve(Dy=2*x*y) (5) 求方程的通解.y=dsolve(Dy=-(exp(x+y)-exp(x)/exp(x+y)+exp(y),x)（6）求方程 的通解y=dsolve(Dy=-(x-y)/(x+y) 2.求解二阶常微分方程（1）y=dsolve(D2y=sin(2*x)-y,y(0)=0,Dy(0)=1,x) （2）计算微分方程 通解。y=dsolve(D2y=-2*Dy-exp(x),x) y=dsolve(D2y=-a*Dy+b*y,x)（4）y=dsolve(D2y=2*Dy-y,x) 3求解微分方程组的解3.解微分方程组4.解微分方程组1. x,y=dsolve(D2y=-Dy+x,Dx=-y,x(0)=1,y(0)=0,t)2. x,y=dsolve(Dy=-2*Dx+4*x+y+exp(t),Dx=-y-3*x,x(0)=3/2,y(0)=0,t)3. x,y=dsolve(Dy=3*y-2*z,Dz=2*y-z,x) 4. x,y=dsolve(D2x=-Dy+x+exp(t),D2y=-Dx-y,t)作业七一不定积分1 int(3x*exp(x),x) 2 int(cos(2*x)/(cos(x)2)*(sin(x)2),x) 3 int(1-1/