Matlab数学实验报告

./WORD格式整理数学实验报告姓名：班级：学号：第一次实验任务过程：a=1+3i;b=2-i;结果：a+b=3.0000+2.0000ia-b=-1.0000+4.0000ia*b=5.0000+5.0000ia/b=-0.2000+1.4000i过程：x=-4.5*pi/180;y=7.6*pi/180;结果：sin<abs<x>+y>/sqrt<cos<abs<x+y>>>=0.2098心得：对于matlab中的角度计算应转为弧度。〔1过程：x=0:0.01:2*pi;y1=sin<x>;y2=cos<x>;y3=exp<x>;y4=log<x>;plot<x,y1,x,y2,x,y3,x,y4>plot<x,y1,x,y2,x,y3,x,y4>结果：〔2过程：>>subplot<2,2,1>>>plot<x,y1>>>subplot<2,2,2>>>plot<x,y2>>>subplot<2,2,3>>>plot<x,y3>>>subplot<>>>subplot<2,2,4>>>plot<x,y4>结果：心得：在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与holdon进行区别,后者为在同一坐标系中划出多条曲线。5、随机生成一个3x3矩阵A及3x2矩阵B,计算〔1AB,〔2对B中每个元素平方后得到的矩阵C,〔3sinB,〔4A的行列式,〔5判断A是否可逆,若可逆,计算A的逆矩阵,〔6解矩阵方程AX=B,〔7矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。过程：A=fix<rand<3,3>.*10>;B=fix<rand<3,3>.*10>;结果：<1>A*B=2228496476100<2>C=B.^2C=149162536<3>sin<B>ans=0.84150.90930.1411-0.7568-0.9589-0.2794<4>a=det<A>a=0<5>由det<A>存在,且det〔A=0,故A不可逆<6>X无解<7>l=A<2,1>,l1=l+1,m=A<2,2>,m1=m+1,n=A<2,3>,n1=n+1D=A,D<2,1>=l1,D<2,2>=m1,D<2,3>=n1D=123567789心得：熟练掌握matlab中对数组的操作方法尤为重要,其中,除了要牢记一些十分方便的函数外,应注意对数组的每一个元素进行操作时,应在运算符前加"."号,如.^.*等。6、设y=<x2+excosx+[x]>/x,分别计算x=1,3,5,7.4时y的值。其中[x]表示x的取整函数。过程：x=1,y=<x^2+exp<x>*cos<x>+fix<x>>/x结果：x=1y=3.4687x=3y=-2.6282x=5y=14.4198x=7.4y=105.29957、已知某地区1-12月份的平均气温为5,8,13,20,25,28,31,33,27,24,19,10,请绘图表示。过程:：T=[5,8,13,20,25,28,31,33,27,24,19,10];

bar<T>结果：第2次实验任务某大学本科生就业情况如下：458人考入研究生,60人签到外企或合资企业,184人签到国营大企业,87人签到私人企业,13人自主创业,画出饼图表示之；过程：x=[458601848713];pie3<x>答案：求极限过程：symsxa;limit<'<<x-a>/<x+a>>^x',x,inf>答案：ans=exp<-2*a>求极限过程：>>symsx;>>limit<'tan<x>^<1/log10<x>>',x,0,'right'>答案：ans=exp〔14.设过程：symsxqt;q=limit<'<1+1/x>^<2*x>',x,inf>;df=diff<'q^t',t>答案：df=q^t*log<q>5.展开多项式；过程：symsab;F=<a+b>^3+<b-1>^2+a+2*b;expand<F>答案：ans=a^3+3*a^2*b+3*a*b^2+b^3+b^2+1+a6.分解因式过程：symsx;F=x^5+3*x^4+2*x^2+x-1;factor<F>答案：ans=x^5+3*x^4+2*x^2+x-17.求方程的根过程：symsx;C=solve<‘x^3-2*x+1=0’>答案C=11/2*5^<1/2>-1/2-1/2-1/2*5^<1/2>8.过程：>>symsx;>>f=sqrt<x+sqrt<x+sqrt<x>>>;>>df=diff<f>答案：df=1/2/<x+<x+x^<1/2>>^<1/2>>^<1/2>*<1+1/2/<x+x^<1/2>>^<1/2>*<1+1/2/x^<1/2>>>9.求不定积分过程：symsx;F=1/<sqrt<2*x+3>+sqrt<2*x-1>>；Nf=int<F,x,>答案：Nf=-1/12*<2*x-1>^<3/2>+1/12*<2*x+3>^<3/2>过程：symsx;F=1/<1+cos<x>>;Nf=int<F>答案：Nf=tan<1/2*x>10求定积分过程：symsx;F=1/<x*sqrt<log<x>*<1-log<x>>>;Nf=int<F,x,exp<1/2>,exp<3/4>>答案：ans=og<2>-1>11.解方程组过程：symsxy;eq1=sym<‘2*x+y=8’>;eq2=sym<‘x-3*y=1’>;[x,y]=solve<eq1,eq2>答案：x=25/7y=6/712.求和过程：symsk;symsum<1/k^2,k,1,20>答案：ans=13.过程：symsx;tay答案：ans=1/2-3^<1/2>*<x-1/6*pi>-<x-1/6*pi>^2+2/3*3^<1/2>*<x-1/6*pi>^3+1/3*<x-1/6*pi>^4-2/15*3^<1/2>*<x-1/6*pi>^5-2/45*<x-1/6*pi>^6+4/315*3^<1/2>*<x-1/6*pi>^7+1/315*<x-1/6*pi>^8-2/2835*3^<1/2>*<x-1/6*pi>^9-2/14175*<x-1/6*pi>^10+4/155925*3^<1/2>*<x-1/6*pi>^11+2/467775*<x-1/6*pi>^12-4/6081075*3^<1/2>*<x-1/6*pi>^13-4/42567525*<x-1/6*pi>^1414.编写猜数游戏程序：首先由计算机随机产生一个[1,100]之间的一个整数,然后由用户猜测所产生的这个数。根据用户猜测的情况给出不同的提示,如果猜测的数大于产生的数,则显示"High",小于则显示"Low",等于则显示"Youwon！",同时退出游戏。用户最多有7次机会。过程：a=fix<rand<1>*100>;fori=1:1:7x=input<'pleaseinputthenumberyouguess'>;ifx>a fprintf<'high'>;continue;endifx<afprintf<'low'>;continue;endifx==afprintf<'youwon'>;break;endend结果：pleaseinputthenumberyouguess50lowpleaseinputthenumberyouguess75highpleaseinputthenumberyouguess63highpleaseinputthenumberyouguess56lowpleaseinputthenumberyouguess60youwon>>心得体会：〔1灵活运用符号方程和各类求解函数对于我们高数的学习和数学问题的验证将会产生很大帮助。〔2应熟练掌握生成某个范围内随机整数等常用技巧。〔3应正确区分循环中的continue与break的作用,即continue是跳过当次循环,而break是跳出整个循环。第三次Matlab上机实验报告第一题：每门课程考试阅卷完毕,任课教师都要对各班的考试成绩进行统计,统计内容包括:全班人数,总得分,平均得分,不及格的人数及90分<包括90分>以上的人数.请编制程序解决这一问题,并自给一组数据验证程序的正确性.要求:使用者在提示下通过键盘输入学生成绩,计算机自动处理后,显示需要的结果.问题分析及程序：1分析：运用循环来对数组中的每个数的值进行判定,从而达到统计各项人数的目的。2程序：a=input<'请输入数组学生成绩

a[n]='>;%输入学生成绩到数组中n=length<a>;

s=0;

p=0;

q=0;

for

i=1:n%在循环中进行统计

s=a<i>+s;

if

a<i>>=90

p=p+1;

elseif

a<i><60

q=q+1;

end

end

r=s/n;

fprintf<'全班人数%.0f

总得分%.0f

平均得分%.2f

不及格的人数%.0f

90分<包括90分>以上的人数%.0f\n',n,s,r,q,p>;%输出运行结果：请输入数组学生成绩a[n]=[8798679890706066]全班人数8总得分636平均得分79.50不及格的人数090分<包括90分>以上的人数3结果分析：灵活运用循环语句,选择语句等基本语句,能很方便地在程序中实现很多基本功能。第二题：1习题1：编程计算小于1000的水仙花数。问题分析及程序：分析：水仙花数指一个n位数<n≥3>,它的每个位上的数字的n次幂之和等于它本身。所以在外层设置1到1000的循环,内层对每一个数是否满足水仙花数的条件进行判断并输出。程序：i=10;g=0;s=0;b=0;fprintf<‘1到1000的水仙花数有：\n’>whilei<1000g=mod<i,10>;s=<mod<i,100>-g>/10;b=<i-10*s-g>/100;ifi==g^3+b^3+s^3fprintf<'%.0f\n',i>;endi=i+1;end运行结果:1到1000的水仙花数有：153370371407结果分析：〔1在1到1000内只有以上四个水仙花数；〔2其实,改程序的代码还可以进行进一步简化,用三个for循环语句的嵌套〔每一个for循环代表某一位上数字的变化来实现对一个三位数的检验。2习题3：求方程2x+y+4z=100的所有非负整数解。问题分析及程序：分析：对于该三元一次方程,无法用现成的数学函数直接求解,应用一个三层循环对x,y,z的取值进行判断。程序：fprintf<'方程的解有：\n'>;i=1;forx=1:50fory=1:100forz=1:25if2*x+y+4*z==100i=i+1;endendendendfprintf<‘%.0f个’,i>;结果：方程的解有：577个分析：总共有577组解〔因数量过多,故只在此显示解的数量。第三题：编写猜数游戏程序：首先由计算机随机产生一个[1,100]之间的一个整数,然后由用户猜测所产生的这个数。根据用户猜测的情况给出不同的提示,如果猜测的数大于产生的数,则显示"High",小于则显示"Low",等于则显示"Youwon！",同时退出游戏。用户最多有7次机会。问题分析及程序：分析：产生随机数,进行七次循环,每次对输入的值进行判定,若相同则提前跳出循环。程序：a=fix<rand<1>*100>;fori=1:1:7x=input<'pleaseinputthenumberyouguess'>;Ifx>a fprintf<'high'>;continue; endIfx<afprintf<'low'>;continue;endifx==afprintf<'youwon'>;break;endend结果：pleaseinputthenumberyouguess50lowpleaseinputthenumberyouguess75highpleaseinputthenumberyouguess63highpleaseinputthenumberyouguess56lowpleaseinputthenumberyouguess60youwon结果分析：用二分法能有较高概率猜中结果。第四题：Fibonacci数组元素满足Fibonacci规则：〔1用循环语句形成前20个分量的数组；〔2求该数组中第一个大于10000的数.问题分析及程序：分析：分别用for循环和while循环能解决以上两个问题。程序：〔1a<1>=1;a<2>=1fori=3:20a<i>=a<i-1>+a<i-2>;end结果：11235813213455891442333776109871597258441816765〔2a<1>=1;a<2>=1;i=3;whilea<i-1><10000a<i>=a<i-1>+a<i-2>;ifa<i>>10000fprintf<‘第一个大于10000的数是：%d’,a<i>>endi=i+1;end结果：第一个大于10000的数是：10946第五题：5.〔1编写一个function函数,调用该函数,可以求一组数的最大值及最小值.