数学建模中的数据处理方法.pptVIP

数值微分MATLAB用命令diff按两点公式计算差分;此题自编程序用三点公式计算相关变化率.编程如下(diff3.m):fori=1:length(x)ifi==1r(1)=(-3*x(1)+4*x(1+1)-x(1+2))/(20*x(1));elseifi~=length(x)r(i)=(x(i+1)-x(i-1))/(20*x(i));elser(length(x))=(x(length(x)-2)-4*x(length(x)-1)+3*x(length(x)))/(20*x(length(x)));endendr=r;第29页，讲稿共79页，2023年5月2日，星期三数值微分保存为diff3.m文件听候调用.再在命令窗内键入X=[1900,1910,1920,1930,1940,1950,1960,1970,1980,1990];x=[76.0,92.0,106.5,123.2,131.7,150.7,179.3,204.0,226.5,251.4];diff3;由于r以离散数据给出,所以要用数值积分计算.键入x(1,1)*exp(trapz(X(1,1:9),r(1:9)))数值积分命令：trapz(x),trapz(x,y),quad(‘fun’,a,b)等.第30页，讲稿共79页，2023年5月2日，星期三微分方程数值解(单摆问题)单摆问题的数学模型是在初始角度不大时,问题可以得到很好地解决,但如果初始角较大,此方程无法求出解析解.现问题是当初始角为100和300时,求出其解,画出解的图形进行比较。第31页，讲稿共79页，2023年5月2日，星期三微分方程数值解(单摆问题)解:若θ0较小,则原方程可用来近似.其解析解为θ(t)=θ0cosωt,.?若不用线性方程来近似,那么有两个模型:第32页，讲稿共79页，2023年5月2日，星期三微分方程数值解(单摆问题)取g=9.8,l=25,100=0.1745,300=0.5236.用MATLAB求这两个模型的数值解,先要作如下的处理:令x1=θ,x2=θ’,则模型变为第33页，讲稿共79页，2023年5月2日，星期三微分方程数值解(单摆问题)再编函数文件(danbai.m)functionxdot=danbai(t,x)xdot=zeros(2,1);xdot(1)=x(2);xdot(2)=-9.8/25*sin(x(1));第34页，讲稿共79页，2023年5月2日，星期三微分方程数值解(单摆问题)在命令窗口键入()[t,x]=ode45(‘danbai’,[0:0.1:20],[0.1745,0]);[t,y]=ode45(‘danbai’,[0:0.1:20],[0.5236,0]);plot(t,x(:,1),’r’,t,y(:,1),’k’);第35页，讲稿共79页，2023年5月2日，星期三优化问题线性规划有约束极小问题?非线性规划有约束极小问题?非线性无约束极小问题?非线性最小二乘问题二次规划第36页，讲稿共79页，2023年5月2日，星期三线性规划有约束极小问题模型?????用命令[x,fval]=linprog(f,A,b,A1,b1,lb,ub)第37页，讲稿共79页，2023年5月2日，星期三线性规划有约束极小问题Findxthatminimizesf(x)=-5x1-4x2-6x3subjecttox1-x2+x3≦203x1+2x2+4x3≦423x1+2x2≦300≦x1,0≦x2,0≦x3第38页，讲稿共79页，2023年5月2日，星期三线性规划有约束极小问题第39页，讲稿共79页，2023年5月2日，星期三线性规划有约束极小问题解问题把问题极小化并将约束标准化第40页，讲稿共79页，2023年5月2日，星期三线性规划有约束极小问题键入c=[-2,-3,5];a=[-2,5,-1];b=-10;a1=[1,1,1];b1=7;LB=[0,0,0];[x,y]=linprog(c,a,b,a1,b1,LB)得当X=(6.4286,0.5714,0.0000)时,z=-14.5714最大.第41页，讲稿共79页，2023年5月2日，星期三线性规划有约束极小问题解问题第42页，讲稿共79页，2023年5月2日，星期三线性规划有约束极小问题解:键入c=[-