数学实验与数学建模实验报告及承诺 书.docx
- 文档编号:10567698
- 上传时间:2023-05-26
- 格式:DOCX
- 页数:36
- 大小:334.05KB
数学实验与数学建模实验报告及承诺 书.docx
《数学实验与数学建模实验报告及承诺 书.docx》由会员分享,可在线阅读,更多相关《数学实验与数学建模实验报告及承诺 书.docx(36页珍藏版)》请在冰点文库上搜索。
数学实验与数学建模实验报告及承诺书
数学实验与数学建模实验报告
学院:
物理科学与技术
专业班级:
物理升华
姓名:
席运泽
学号:
1104090119
完成时间:
2009年12月日
承诺书
本人承诺所呈交的数学实验与数学建模作业都是本人通过学习自行进行编程独立完成,所有结果都通过上机验证,无转载或抄袭他人,也未经他人转载或抄袭。
若承诺不实,本人愿意承担一切责任。
承诺人:
席运泽
2009年12月26日
注:
所有结果交电子打印稿,并将电子文档发送到邮箱:
xuanyunqin@
(word文档命名:
姓名+学号+数学实验作业)
数学实验学习体会
(每个人必须要写,占总成绩的20%)
通过几周的学习,我对MATLAB数学实验与建模有了更深的认识。
我感受到MATLAB强大的运算能力和实用性,为以后的数学建模竞赛打下了坚实的基础。
但最深刻的感受就是:
要不断的用它。
Matlab是个好工具,但如若不用它来解决问题,只知道一点语法,那是连皮毛都没有学到。
还有就是程序的设计,对于程序的运行效率非常有帮助。
有的时候,编出来的程序,能够运行,但是耗时太长,程序虽然没有错,但是不适合实际。
这就需要对程序的结构和算法问题进行改进,要时刻思考多动脑,找到适宜的解题途径。
还有就是学习MATLAB要多动手,找一个习题实际操作一下或者找一个实际的程序来动手编一下,才能更好的对MATLAB有所了解,进一步巩固知识。
你要在编程的过程中学习,程序需要什么知识再去补充,编程是一点一点积累的,需做一些随手笔记,我就是在这个时候有所懈怠才发懵的。
当然,除了要去用它以外,辅导书可以很大程度上提高我们的知识与技能,通过模仿别人编写的程序,可以大大加快我们掌握它的进度,并且学到一些课堂中所没有的知识。
实验1MATLAB运算基础
1.先求下列表达式的值,然后显示MATLAB工作空间的使用情况并保存全部变量。
(1)
>>x=exp
(2)
定义变量x=
7.3891
>z=2*sin(85*pi/180)/(1+x)
z=
0.2375
(2)
,其中
>>
则x=[2,1+2i;-0.45,5];>>z=1/2*log(x+sqrt(1+x^2))
运算z=
0.7114-0.0253i0.8968+0.3658i
0.2139+0.9343i1.1541-0.0044i
(3)
,
则a=-3.0:
0.1:
3.0;>>
z=(exp(0.3*a)-exp(-0.3*a)).*sin(a+0.3)/2
z=
Columns1through13
0.43870.50720.56400.60890.64240.66480.67660.67840.67090.65490.63130.60110.5652
Columns14through26
0.52470.48050.43370.38540.33660.28810.24100.19620.15430.11620.08240.05350.0299
Columns27through39
0.0120-0.0000-0.0060-0.006000.01170.02880.05090.07750.10800.14180.17800.2159
Columns40through52
0.25470.29340.33110.36690.39980.42890.45320.47190.48410.48910.48620.47480.4543
Columns53through61
0.42430.38460.33500.27550.20610.12720.0392-0.0574-0.1619
(4)
其中
程序:
>>fort=0:
0.5:
2.5;if0<=t<=1
z=t^2
elseif1 z=t^2-1 elsez=t.^2-2*t+1 end end z= 0 z= 0.2500 z= 1 z= 2.2500 z= 4 z= 6.2500 2.已知 求下列表达式的值: (1)A+6*B和A-B+I(I为单位矩阵)。 (2)A*B和A.*B (3)A^3和A.^3 (4)A/B和B\A (5)[A,B]和[A([1,3],: );B^2] 1. >>A=[1234-4;34787;3657]; >>B=[13-1;203;3-27]; >>C=A+6*B C= 1852-10 467105 215349 >>D=A*B D= 684462 309-72596 154-5241 >>E=A.*B E= 121024 680261 9-13049 >>F=A^3 F= 3722623382448604 247370149188600766 78688454142118820 >>G=A.^3 G= 172839304-64 39304343658503 27274625343 >>H=A/B H= 16.4000-13.60007.6000 35.8000-76.200050.2000 67.0000-134.000068.0000 >>I=B/A I= 0.1027-0.0062-0.0069 0.06170.0403-0.0366 0.02050.0855-0.0507 >>J=[A,B] J= 1234-413-1 34787203 36573-27 >>K=[A([1,3],: );B^2] K= 1234-4 3657 451 11019 20-540 3.设有矩阵A和B , (1)求他们的乘积C。 (2)求矩阵C的右下角3x2子矩阵赋给D。 (3)查看MATLAB工作空间的使用情况。 >>C=A*B C= 9315077 258335237 423520397 588705557 753890717 >>D=C(3: 5,2: 3) D= 520397 705557 890717 4.完成下列操作: (1)求[100,999]之间能被21整除的数的个数。 (2)建立一个字符串向量,删除其中的大写字母。 实验2选择结构程序设计 1.求下列分段函数的值 编程: >>x=input(': ') if-3 y=x^2+x-6 elseifx<-3 y=x^2+x-6 elseif0<=x<2 y=x^2-5*x+6 elseif2 y=x^2-5*x+6 elseif3 y=x^2-5*x+6 else y=x^2-x-1 end : -1 x= -1 y= 12 2.输入一个百分制成绩,要求输出成绩等级A,B,C,D,E。 其中90~100分为‘A’,80~89分为‘B’,70~79分为‘C’,60~69分为‘D’,60分以下为‘E’。 要求: (1)分别用if语句和switch语句实现。 (2)输入百分制成绩后要判断该成绩的合理性,对不合理的成绩应输出错误成绩。 (1)x=input(': ') if90<=x<=100 A=x elseif80<=x<=89 B=x elseif70<=x<=79 C=x elseif60<=x<=69 D=x elseE=x end : 98 x= 98 A= 98 Switch语句: >x=input(': ') switch case90<=x<=100 A=x case80<=x<=89 B=x case70<=x<=79 C=x case60<=x<=69 D=x otherwiseE=x end 3.假定某地区电话收费标准为: 通话时间3分钟以下,收费0.5元;3分钟以上,则每超过分钟加收0.15元;通话时间在7: 00~22: 00之间者,按上述标准全价收费,在其他通话时间者,按上述标准半价收费。 计算某人在t1到t2通话时间,应缴多少电话费 程序如下: 4.建立一个5×6矩阵,要求输出矩阵第n行元素。 当n值超过矩阵的行数时,自动转为输出矩阵最后一行元素,并给出错误信息。 实验三循环结构程序设计 1.根据 ,求π的近似值。 当n分别取100、1000、1000时,结果是多少? i=input(': ') sum=0;n=1000000 fori=1: n sum=sum+i*i end : 100 sum= 3.3234e+017 2.根据 求: (1)y<3时的最大值。 (2)与 (1)的n值对应的y值 3一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。 试输出全部的水仙花数。 要求 (1)用循环结构实现。 (2)用向量运算实现。 4.已知 求f1~f100中: (1)最大值、最小值、各数之和。 (2)正数、零、负数的个数。 试验四函数与文件 1、一个自然数是素数,且它的各位数字位置经过任意对换后仍为素数,则称为绝对素数。 试求所有的两位绝对素数。 实验五高层绘图操作 1.画出空间曲线 在 范围内的图形,并画出相应的等高线。 [x,y]=meshgrid(-30: 0.5: 30); z=10*sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+1); meshc(x,y,z); 2.根据给定的参数方程,绘制下列曲面的图形。 a)椭球面 >>symsuv u=0: 0.2: 2*pi; v=0.001: 0.1: 2; [u,v]=meshgrid(u,v); z=sin(u); y=2*cos(u).*cos(v); >>x=3*cos(u).*sin(v); >>meshc(x,y,z) b)椭圆抛物面 symsuv u=0: 0.2: 2*pi; v=0.001: 0.1: 2; [u,v]=meshgrid(u,v); z=4*u.^2; y=2*u.*cos(v); x=3*u.*sin(v); mesh(x,y,z) c)单叶双曲面 symsuv u=0: 0.2: 2*pi; v=0.001: 0.1: 2; [u,v]=meshgrid(u,v); z=4*tan(u); y=2*cos(v).*sec(u); x=3*sec(u).*sin(v); mesh(x,y,z) d) 双曲抛物面 symsuv u=-5: 0.1: 5 v=0.001: 0.1: 2; [u,v]=meshgrid(u,v); z=(u.^2-v.^2); y=v; x=u; mesh(x,y,z) e)旋转面 symsuv; u=0: 0.2: 4*pi; v=0.001: 0.1: 2; [u,v]=meshgrid(u,v); z=u; y=log(u).*cos(v); x=log(u).*(sin(v)); mesh(x,y,z) f)圆锥面 symsuv; u=0: 0.2: 4*pi; v=0.001: 0.1: 2; [u,v]=meshgrid(u,v); z=u; y=u.*cos(v); x=u.*(sin(v)); mesh(x,y,z); view(125,-9) g)环面 symsuv; u=0: 0.2: 4*pi; v=0.001: 0.1: 2; [u,v]=meshgrid(u,v); z=0.4*sin(v); x=(3+0.4*cos(u)).*cos(v); y=(3+0.4*cos(u)).*(sin(v)); mesh(x,y,z); view(125,-9) h)正螺面 symsuv; u=0: 0.2: 4*pi; v=0.001: 0.1: 2; [u,v]=meshgrid(u,v); z=4*v; y=u.*cos(v); x=u.*(sin(v)); mesh(x,y,z); view(125,-9) 3.已知 ,完成下列操作: (1)在同一坐标系下用不同颜色和线形绘制3条曲线。 (2)以子图形式绘制3条曲线。 (3)分别用条形图、阶梯图、杆图和填充图绘制3条曲线。 1.x=(0: pi/100: 2*pi)'; y1=x.^2; y2=cos(2*x); y3=y1.*y2; plot(x,y1,'g: ',x,y2,'b--',x,y3,'rp'); 2. x=-2*pi : 0.1 : 2*pi; y1=x.^2;y2=cos(2*x); y3=y1.*y2 subplot(2,2,1); plot(x,y1);title('x.^2');axis([-2*pi,2*pi,-30,30]); subplot(2,2,2); plot(x,y2);title('cos(2*x)');axis([-2*pi,2*pi,-30,30]); subplot(2,2,3); plot(x,y3);title('y3');axis([-2*pi,2*pi,-30,30]); 3.x=0: 0.2: 7; y=x.^2; subplot(2,2,1);bar(x,y,'g'); title('bar(x,y,''g'')');axis([0,7,0,50]); subplot(2,2,2);fill(x,y,'r'); title('fill(x,y,''r'')');axis([0,7,0,50]); subplot(2,2,3);stairs(x,y,'b'); title('stairs(x,y,''b'')');axis([0,7,0,50]); subplot(2,2,4);stem(x,y,'k'); title('stem(x,y,''k'')');axis([0,7,0,50]); x=0: 0.1: 7; y=cos(2*x); subplot(2,2,1);bar(x,y,'g'); title('bar(x,y,''g'')');axis([0,7,0,2]); subplot(2,2,2);fill(x,y,'r'); title('fill(x,y,''r'')');axis([0,7,0,2]); subplot(2,2,3);stairs(x,y,'b'); title('stairs(x,y,''b'')');axis([0,7,0,2]); subplot(2,2,4);stem(x,y,'k'); title('stem(x,y,''k'')');axis([0,7,0,2]); x=0: 0.1: 7; y=cos(2*x).*(x.^2); subplot(2,2,1);bar(x,y,'g'); title('bar(x,y,''g'')');axis([0,7,0,10]); subplot(2,2,2);fill(x,y,'r'); title('fill(x,y,''r'')');axis([0,7,0,10]); subplot(2,2,3);stairs(x,y,'b'); title('stairs(x,y,''b'')');axis([0,7,0,10]); subplot(2,2,4);stem(x,y,'k'); title('stem(x,y,''k'')');axis([0,7,0,10]); 4.绘制极坐标曲线 ,并分析参数 对曲线形状的影响。 symsxpabn x=-2*pi: 0.1: 2*pi;p=a*sin(b+n*x);polar(x,p) 5.绘制函数的曲面图和等高线。 (1) x=-10: 0.1: 10; [x,y]=meshgrid(x); z=(x.^2-2.*x).*exp(-x.^2-y.^2-x.*y); mesh(x,y,z) [x,y,z]=peaks; contour3(x,y,z,12,'k'); (2) 程序: x=-5: 0.2: 5; [x,y]=meshgrid(x); f=1./sqrt((x-1).^2+y.^2)-1./sqrt((x+1).^2+y.^2)+eps; mesh(x,y,f); contour3(x,y,f,12,'k'); 实验六数据处理与多项式计算 1、将100个学生5门功课的成绩存入矩阵P中,进行如下处理: (1)分别求5门课的最高分、最低分和相应学生的学号。 (2)分别求5门课的平均分和标准方差。 (3)5门课总分的最高分、最低分及相应学生学号。 (4)将5门课总分按从大到小的顺序存入zcj中,相应学生学号存入xsxh。 2、某气象观测站测得某日6: 00~18: 00之间每隔2小时室外温度如下表所示: 时间h 6 8 10 12 14 16 18 室外温度t1 18.0 20.0 22.0 25.0 30.0 28.0 24.0 室外温度t2 15.0 19.0 24.0 28.0 34.0 32.0 30.0 试用三次样条差值分别求出该日室内外6: 00~17: 30时之间每隔两小时各点的近似温度。 3、已知lg(x)在[1,101]区间11个整数采样点的函数值如下表所示。 x 1 11 21 31 41 51 61 71 81 91 101 Lg(x) 0 1.0414 1.3222 1.4914 1.6128 1.7076 1.7853 1.8513 1.9085 1.9590 2.0043 试求lg(x)的5次拟合多项式p(x),并分别绘制出lg(x)和p(x)在[1,101]区间上的函数曲线。 4、有三个多项式 , , ,试进行下列操作: (1)求 。 P1=[12405];p2=[12];p3=[123];p=conv(p2,p3)+p1;r=roots(p) (2)求 的根。 (3)当x取矩阵A的每一元素时,求的值。 其中 (4)当以矩阵A为自变量时,求 的值。 实验七符号计算基础和符号微积分 1、已知 ,利用符号表达式求 2、已知 求: (1)B=P1*P2*A。 (2)B的逆矩阵并验证结果。 (3)包括B矩阵主对角线元素的下三角阵。 (4)B的行列式值。 3、用符号法求下列极限和导数。 (1) (2)已知 ,求 4、用符号法求下列积分。 (1) x=sym(‘x’);f=1/(1+x^4+x^8);int(f) (2) x=sym(‘x’);f=(1+x^2)/(1+x^4);int(f,x,0,inf) 实验八级数和符号方程求解 1、级数符号求和。 (1)计算 symsx;s1=symsum(1/(2*n+1),n,1,10) s1= 17181176/14549535 (2)求级数 之和函数,并求 之和 symsn;symsxreal; s2=symsum(x^(n-1)*n^2,n,1,inf) s2= piecewise([abs(x)<1,-(x^2+x)/(x*(x-1)^3)]) symsx;s3=symsum(n^2/5^n,n,1,inf) s3= 15/32 2、将ln(x)在x=1处按5次多项式展开为泰勒级数。 x=sym(‘x’);y=log(x);seriestaylor(y,x,6) 3、求微分方程的符号解。 3、求微分方程的符号解。 >>symsyxab [xy]=solve(‘(y.^2)/(x.^2)+(k.^2)*y=0’,’y(0)=0’,‘y’’=o,‘x,y’) 4、求下列方程和方程组的符号解。 (1) x=solve('3*x*exp(x)+5*sin(x)-78.5=0','x') x= 2.3599419584772910151699327715486 (2) [xy]=solve('sqrt(x^2+y^2)=100','3*x+5*y-8=0','x,y') x= (10*21246^(1/2))/17+12/17 12/17-(10*21246^(1/2))/17 y= 20/17-(6*21246^(1/2))/17 (6*21246^(1/2))/17+20/17 实验六数据处理与多项式计算 2、将100个学生5门功课的成绩存入矩阵P中,进行如下处理: (5)分别求5门课的最高分、最低分和相应学生的学号。 (6)分别求5门课的平均分和标准方差。 (7)5门课总分的最高分、最低分及相应学生学号。 (8)将5门课总分按从大到小的顺序存入zcj中,相应学生学号存入xsxh。 已知矩阵P(100*5) 最高分及学生学号: [Y,U]=max(P) 最低分及学生学号: [M,N]=min(P) 平均分[A,B]=mean(P) 标准方差: std(A,0,1) 总分最高: [C,D]=max(consum(A(1: 100,;)) 总分最低: [E,F]=min(consum(A(1: 100,;)) 排列: [zcj,xsxh]=sort(consum(A(1: 100,;)) 2、某气象观测站测得某日6: 00~18: 00之间每隔2小时室外温度如下表所示: 时间h 6 8 10 12 14 16 18 室外温度t1 18.0 20.0 22.0 25.0 30.0 28.0 24.0 室外温度t2 15.0 19.0 24.0 28.0 34.0 32.0 30.0 试用三次样条差值分别求出该日室内外6: 00~17: 30时之间每隔两小时各点的近似温度。 >>h=6: 2: 18 t1=[18.0,20.0,22.0,25.0,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学实验与数学建模实验报告及承诺 数学 实验 建模 报告 承诺