《数学模型实验报告65992》由会员分享，可在线阅读，更多相关《数学模型实验报告65992（10页珍藏版）》请在金锄头文库上搜索。
1、数学模型实验报告姓名专业：数学与应用数学学号：实验一 熟悉Matlab软件基本功能一、实验目的和要求： 熟悉Matlab软件的基本功能，数组的输入及运算。二、实验内容和步骤：1、 熟悉Matlab软件的窗口的功能。2、 熟悉Matlab软件工作区以及基本命令。3、 熟悉Matlab软件帮助系统。4、 熟悉数组运算。5、 (1) a=1 2 3,b=5 7 8 c=4,计算：a+c; a*c; a+b; a.*b; a ./b ; a.b ; a.c ;c.a ; a.b .(2) A=1 2 3; 2 3 4; 3 4 5,B=1 1 1; 2 2 2; 3 3 3, c=3,计算A+c; A
2、+B; A*B; A.*B; A.B; A.c; Ac; A.c. 6、在同一坐标系中作出 的图形。三、实验现象及结果：解：5】（1） a=1 2 3; b=5 7 8; c=4; a+cans =5 6 7 a*cans = 4 8 12 a+bans = 6 9 11 a.*bans =5 14 24 a./bans =0.2000 0.2857 0.3750 a.bans =5.0000 3.5000 2.6667 a.cans = 1 16 81 c.aans = 4 16 64 a.bans = 1 128 6561（2） A=1 2 3; 2 3 4; 3 4 5; ,B=1 1
3、1; 2 2 2; 3 3 3; B=1 1 1; 2 2 2; 3 3 3; c=3; A+cans =4 5 6 5 6 76 7 8 A+Bans = 2 3 4 4 5 6 6 7 8 A*Bans = 14 14 14 20 20 20 26 26 26 A.*Bans = 1 2 3 4 6 8 9 12 15 A.Bans = 1.0000 0.5000 0.3333 1.0000 0.6667 0.5000 1.0000 0.7500 0.6000 A.cans =3.0000 1.5000 1.0000 1.5000 1.0000 0.75001.0000 0.7500 0.
4、6000 Acans =132 192 252 192 279 366 252 366 480 A.cans = 1 8 27 8 27 64 27 64 1256】 x=0:2*pi/30:2*pi; y1=sin(x);y2=cos(x); plot(x,y1,x,y2,:) 实验二 熟悉Matlab软件M文件一、实验目的和要求： 熟悉Matlab软件的命令式文件和函数式文件，熟悉循环和选择控制结构。二、实验内容和步骤：1、熟悉Matlab软件的程序编器。2、熟悉命令式文件。3、熟悉函数式文件。4、(1) 编制循环和选择控制结构程序并运行 (2） 编制一个命令式文件，计算5！。（3）编制一
5、个命令式文件，计算100以内质数之和。 （4） 编制一个函数式文件,参数n，返回n!。三、实验现象及结果：4】（1）循环结构：计算（1/n/n）（n t=0; n=1; while(n tt = 1.6350选择控制结构：计算y=cos(x+1)+sqrt(x2+1) x=10y=x*sqrt(x+sqrt(x); x=102 x=10; if x=10y=cos(x+1)+sqrt(x2+1);elsey=x*sqrt(x+sqrt(x);end yy = 10.0543 x=5; if x=10y=cos(x+1)+sqrt(x2+1);elsey=x*sqrt(x+sqrt(x);end
6、 yy = 13.4500(2)M文件：f1.mn=5;a=1;for i=1:n a=a*i;endaf1a=120(3)M文件：f2.ma=0;for i=1:100 if isprime(i) a=a+i; endendaf2a=1060(4)M文件：fun.mfunction y=fun(n)y=1;for i=1:n y=y*i;end n=10; fun(n)ans = 3628800 n=5; fun(n)ans = 120实验三 席位分配一、实验目的和要求： 编制程序，解决席位分配问题。二、实验内容和步骤：1、掌握课本席位分配的几个方法。2、编制程序求解课本席位分配问题。3、分
7、析实验结果。三、实验现象及结果：P283 M文件：f1.mclear;clc;a=9061 7179 5259 3319 1182N=26L=length(a);b=a/sum(a)*N;GR=floor(b);jj=N-sum(GR); bb,xi=sort(b-GR,descend);for i=1:jj; GR(xi(i)=GR(xi(i)+1;endGRHu=zeros(L,5); for n=0:N-1; D(n+1,:)=n+1 n+1/2 sqrt(n*(n+1) 2*n*(n+1)/(2*n+1) n; endfor k=1:2 for i=1:N B(i,:)=a/D(i,k
8、); end bb=B(:); bbb,x=sort(bb,descend); xx=x(1:N); for i=1:L Hu(k,i)=sum(xx-(i-1)*N0)-sum(xx-i*N0); endend N=N-L; for k=3:5 for i=1:N B(i,:)=a/D(i+1,k); end bb=B(:); bbb,x=sort(bb,descend); xx=x(1:N); for i=1:L Hu(k,i)=sum(xx-(i-1)*N0)-sum(xx-i*N0); end end Hu(3:5,:)=Hu(3:5,:)+ones(3,L); GD=Hu(1,:)
9、MF=Hu(2,:) EP=Hu(3,:) HM=Hu(4,:) SD=Hu(5,:)f1a =9061 7179 5259 3319 1182N =26GR =9 7 5 4 1GD =10 7 5 3 1MF =9 8 5 3 1EP =9 7 6 3 1HM = 9 7 5 4 1SD = 9 7 5 3 2P285M文件：f2.mclear;clc;a=91490 1660 1460 1450 1440 1400 1100N=100L=length(a);b=a/sum(a)*N;GR=floor(b);jj=N-sum(GR); bb,xi=sort(b-GR,descend);fo
10、r i=1:jj; GR(xi(i)=GR(xi(i)+1;endGRHu=zeros(L,5); for n=0:N-1; D(n+1,:)=n+1 n+1/2 sqrt(n*(n+1) 2*n*(n+1)/(2*n+1) n; endfor k=1:2 for i=1:N B(i,:)=a/D(i,k); end bb=B(:); bbb,x=sort(bb,descend); xx=x(1:N); for i=1:L Hu(k,i)=sum(xx-(i-1)*N0)-sum(xx-i*N0); end end N=N-L; for k=3:5 for i=1:N B(i,:)=a/D(i
11、+1,k); end bb=B(:); bbb,x=sort(bb,descend); xx=x(1:N); for i=1:L Hu(k,i)=sum(xx-(i-1)*N0)-sum(xx-i*N0); end end Hu(3:5,:)=Hu(3:5,:)+ones(3,L); GD=Hu(1,:) MF=Hu(2,:) EP=Hu(3,:) HM=Hu(4,:) SD=Hu(5,:)f2a = Columns 1 through 5 91490 1660 1460 1450 1440 Columns 6 through 7 1400 1100N =100GR = 92 2 2 1 1 1 1GD =94 1 1 1 1 1 1MF =93 2 1 1 1 1 1EP =90 2 2 2 2 1 1HM
