大专数学建模必须会的知识概要1.doc

1、微分方程的求解 Matlab上机学习指导1.doc 四 符号计算 5 解微分方程 例1：解微分方程 y=dsolve(D2y+2*Dy+2*y=0) 结果： y = C1*exp(-t)*sin(t)+C2*exp(-t)*cos(t) %默认自变量为t 例2：解微分方程，， y=dsolve(D2y+2*Dy+2*y=0,y(0)=1,Dy(0)=0) 结果： y = exp(-t)*sin(t)+exp(-t)*cos(t) 对于上面的符号函数y，可以用下面的命令画出函数图像 ezplot(y) 2、线性规划模型的求解 min=x+y; 4.5*x+10*y100; 51*(x+y)1000 x-60; y40; @free(x);@free(y); (2) model: max=x+y; 4.5*x+10*y100; 51*x+51*y1000; x-60; y40; @GIN( X); @GIN( y); End （3）用LINGO求解矩阵 model: sets: jh/1,2,3,4/:; link(jh,jh):a,x; endsets data: a=15,18,21,24 19,23,22,18 26,17,16,19 19,21,23,17; enddata min=@sum(link(i,j):a(i,j)*x(i,j)); @for(jh(i):@sum(jh(j):x(i,j))=1); @for(jh(j):@sum(jh(i):x(i,j))=1); @for(link(i,j): @BIN( x(i,j))); （4）LINGO求解集合最小值 model: sets: jh2/1,2,3,4,5,6/:x; endsets data: x=5,6,3,4,1,6; enddata y=@min(jh2(j)|j#LE#5:x); 注释： LE:小于等于 GT:大于 LINGO中的函数（edit）： （5）LINGO中矩阵的运算 model: sets: jh1/1..5/:h,q; jh2/1,2,3/:l,c; link(jh1,jh2):x,y,z; endsets data: x=5,6,3 1,2 3 9,7,5 5,6,7 4,5,6; y=1,2,3 5,6,7 7,6,9 8,5,6 9,5,3; q=1 2 3 4 5; enddata @for(link(i,j):z(i,j)=x(i,j)+y(i,j)); @for(jh2(j):l(j)=@sum(jh1(i):x(i,j))); @for(jh1(i):h(i)=@sum(jh2(j):x(i,j))); @for(jh2(j):c(j)=@sum(jh1(i):q(i)*x(i,j))); (1)定义两个5行3列的矩阵，求两个矩阵的和。 (2)求其中一个矩阵前三行数的最大值。 (3)对其中一个矩阵的行求和及列求和。 (4)用一个1行5列的矩阵乘以其矩阵。 3、matlab制图 实验三 MATLAB的图形功能.doc Matlab上机学习指导1.doc 三 函数作图 4、matlab程序设计 （1）选择结构 例5：定义函数， 计算，，，， 建立函数文件fff.m function y=fff(x) if x2 y=x+1; elseif x=2 x=8 y=3*x; elseif x8 x=20 y=4*x-5; else y=cos(x)+sin(x); end 在命令窗口中计算 y=[fff(0.1),fff(1),fff(9),fff(22),fff(2*pi)] 结果y = 1.1000 2.0000 31.0000 -1.0088 18.8496 （2）循环结构 例6：一个简单的for循环示例 for i=1:10; % i依次取1,2,…10,. x(i)=2*i; % 对每个i值，重复执行该指令 end; % 表示循环结束，每一个for要对应一个end x % 要求显示运行后数组x的值。 输入后观察结果，体会for语句的作用。 while循环语句 例7: Fibonacci 数列：1，1，2，3，5，8，… 即：，（ 1，2，3…）现要求该数列中第一个大于10000 的元素。 a(1)=1;a(2)=1;i=2; while a(i)=10000 % 当a(i)=10000时执行后面的语句，否则跳出循环 a(i+1)=a(i-