Matlab简单实例学习

1、Matlab程序代码绘制的函数图象。 fvclear;t=0:0.02:10;f1=10/sqrt(7.75).*exp(-1.5*t);f2=sin(sqrt(7.75).*t);y=f1.*f2;plot(t,y,-k,t,y,ok);xlabel(t);ylabel(y(t) );title(函数图像)axis(-2 10 -0.5 2)拉氏变换clear;clc;syms s t fs1 fs2 fs3 ft1 ft2 ft3;L=1,C=0.1,R=1.5 3 5;h1=1/(L*C*s2+R(1)*C*s+1);h2=1/(L*C*s2+R(2)*C*s+1);h3=1/(L*C*

2、s2+R(3)*C*s+1);fs1=h1*(1/s);fs2=h2*(1/s);fs3=h3*(1/s);ft1=ilaplace(fs1,s,t);ft2=ilaplace(fs2,s,t);ft3=ilaplace(fs3,s,t);ezplot(t,ft1);hold on;ezplot(t,ft2);hold on;ezplot(t,ft3);信号编码对1 1 0 1 1 1 0 1 0 0 1进行编码。clear;clc;c=1 1 0 1 1 1 0 1 0 0 1for i=1:length(c) if i=1 d1(i)=0;d2(i)=0; elseif i=2 d1(i)

3、=c(i-1);d2(i)=c(i-1); elseif i=3 d1(i)=mod(c(i-1)+c(i-2),2); d2(i)=c(i-1); else d1(i)=mod(c(i-1)+c(i-2),2); d2(i)=mod(c(i-1)+c(i-3),2); endendd1d2迭代法使用一般迭代法求解方程的解。第一根：clear;clc;x0=10;err=1;while (err10-6) F1=x0; F2=log(3*x02); err=abs(F1-F2); x0=F2;endx0第二根：clear;clc;x=1;err=1;while (err10-6) F1=x;

4、F2=sqrt(exp(x)/3); err=abs(F1-F2); x=F2;endx第三根：clear;clc;x=10;err=1;while (err10-6) F1=x; F2=-sqrt(exp(x)/3); err=abs(F1-F2); x=F2;endx牛顿迭代法使用牛顿迭代法求解方程的解。syms x fx fx1;fx=2*x3-4*x2+3*x-6;fx1=diff(fx)err=1;k=1;x0=1.5;while (err10-6 & k10-6 & k10-6) & (k10-6) & (k10-6) & (err210-6) & (err310-6) & (kL

5、egend可加标注。最小二乘拟合function y=f(x) %f文件定义t=0.25 0.5 1 1.5 2 3 4 6 8;c=19.21 18.15 15.36 14.10 12.98 9.32 7.45 5.24 3.01;y=c-x(1).*exp(-x(2).*t);clear;clc; %主程序t=0.25 0.5 1 1.5 2 3 4 6 8;c=19.21 18.15 15.36 14.10 12.98 9.32 7.45 5.24 3.01;p,s=polyfit(t,c,7); %求七阶多项式拟合y1=polyval(p,t);r1=corrcoef(c,y1)x0=

6、1 2;x=lsqnonlin(f,x0) %最小二乘拟合xy2=x(1).*exp(-x(2).*t);r2=corrcoef(c,y2)v=0:0.01:8;y11=polyval(p,v);plot(t,c,ok,v,y11,b,t,y2,r);legend(原函数,多项拟合,最小二乘拟合);最小二乘拟合function y=f1(x) %定义f1t=0.5 1.0 1.5 2.0 2.5 3.0;c=1.75 2.45 3.81 4.80 8.00 8.60;y=c-x(1).*exp(-x(2).*t);clear;clc; %主程序t=0.5 1.0 1.5 2.0 2.5 3.0

7、;c=1.75 2.45 3.81 4.80 8.00 8.60;p,s=polyfit(t,c,7);y1=polyval(p,t);r1=corrcoef(c,y1)x0=1 2;x=lsqnonlin(f1,x0)y=x(1).*exp(-x(2).*t);r2=corrcoef(c,y)v=0.5:0.01:3;y11=polyval(p,v);y2=x(1).*exp(-x(2).*v);plot(t,c,o,v,y11,r,v,y2,g);legend(数据,七次多项式拟合,最小二乘拟合);线性和非线性规划有约束优化求解函数在约束条件下最值。function f=f(x)f=100

8、*(x(2)-x(1)2).2+(1-x(1)2;function f1=f(x)f1=-(100*(x(2)-x(1)2).2+(1-x(1)2);clear;clc;x0=1,2;A=1 1;-1 -1;b=1.5;0;x=fmincon(f,x0,A,b)f=100*(x(2)-x(1)2).2+(1-x(1)2x=fmincon(f1,x0,A,b)f=100*(x(2)-x(1)2).2+(1-x(1)2function f=f(x)f=-2*x(1)-x(2);function g,ceq=nonlcon(x)g(1)=x(1)2+x(2)2-25;g(2)=x(1)2-x(2)2

9、-7;ceq=;clear;clc;VLB=0,0;VUB=5,10;x0=1 2;x=fmincon(f,x0,VLB,VUB,nonlcon,)min=-2*x(1)-x(2)function f=f(x)f=exp(x(1).*(4*x(1)2+2*x(2)2+4*x(1)*x(2)+2*x(2)+1);function g,ceq=nonlcon(x)g(1)=1.5+x(1)*x(2)-x(1)-x(2);g(2)=-x(1)*x(2)-10;ceq=x(1)+x(2);clear;clc;x0=1 2;x=fmincon(f,x0,nonlcon,)min=exp(x(1).*(4

10、*x(1)2+2*x(2)2+4*x(1)*x(2)+2*x(2)+1)function f=f(x)f=exp(x(1)+x(2)+x(3)+x(4)+x(5);function g,ceq=nonlcon(x)g(1)=x(3)+x(4)-3;g(2)=x(1)+x(5)-6;ceq(1)=x(1)2+x(2)2+x(3)2+x(4)2+x(5)2-10;ceq(2)=x(2)*x(3)-5*x(4)*x(5);ceq(3)=x(1)+x(2)+1;clear;clc;x0=1 2 3 4 5;VLB=-2.3 -2.3 -2.3 -2.3 -2.3;VUB=2.3 2.3 3.2 3.2

11、 3.2;x=fmincon(f,x0,VLB,VUB,nonlcon,)min=exp(x(1)+x(2)+x(3)+x(4)+x(5)求解clear;clc;f=-1 -3 -3 0 0 0;Aeq=3 1 2 0 1 0;1 0 1 0 2 1;1 0 2 1 2 0;Beq=5;2;6;LB=1e-6 1e-6 1e-6 1e-6 1e-6 1e-6;UB=Inf Inf Inf Inf Inf Inf;format long;x=linprog(f,Aeq,Beq,LB,UB)min=-x(1)-3*x(2)-3*x(3)现有三块地，分别为I=100hm2、II=300hm2、III

12、=200hm2，需要种植水稻、大豆和玉米，要求最低收获水稻190000kg、大豆130000kg、玉米350000kg，三种作物的但产量如下图所示，市价水稻￥1.2/kg、大豆￥1.5/kg、玉米￥0.8/kg。问：1、获得最大产量的种植方案。2、获得最大产值的种植方案。Kg/hm2IIIIII水稻1100095009000大豆800068006000玉米140001200010000clear;clc;f1=-11000 8000 14000 9500 6800 12000 9000 6000 10000;f2=-11000*1.2 8000*1.5 0.8*14000 1.2*9500 1

13、.5*6800 0.8*12000 1.2*9000 1.5*6000 0.8*10000;A=-11000 0 0 9500 0 0 9000 0 0;0 8000 0 0 6800 0 0 6000 0;0 0 14000 0 0 12000 0 0 10000;B=-190000;130000;350000;Aeq=1 1 1 0 0 0 0 0 0;0 0 0 1 1 1 0 0 0;0 0 0 0 0 0 1 1 1;Beq=100;300;200;LB=zeros(9,1);UB=100 100 100 300 300 300 200 200 200;x,fval=linprog(

14、f1,A,B,Aeq,Beq,LB,UB);xmax1=-fvalx,fval=linprog(f2,A,B,Aeq,Beq,LB,UB);xmax2=-fval某产品1900到2010年每个10年的产量为：75.995,91.972,105.711,123.203,131.699,150.696,179.323,203.212,226.505,249.633,256.344,267.893，使用差值法求解1995年的产量，并画出相关图像。clear;clc;x=1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 2000 2010;y=75.9

15、95,91.972,105.711,123.203,131.699,150.696,179.323,203.212,226.505,249.633,256.344,267.893;Y=interp1(x,y,1995,spline)t=1900:0.1:2010;y1=interp1(x,y,t,spline);plot(t,y1,r,x,y,ok);legend(差值函数,原始数据);1950到1990年工龄10、20、30年工作人员劳动报酬如下：服务年限年份1020301950150.697169.592187.6521960179.323195.072250.2871970203.212

16、239.092322.7671980226.505273.706426.7301990249.633370.281598.243计算1975年时15年工龄的工作人员平均工资。clear;clc;X=1950 1960 1970 1980 1990;1950 1960 1970 1980 1990;1950 1960 1970 1980 1990;Y=10 10 10 10 10;20 20 20 20 20;30 30 30 30 30;Z=150.697 179.323 203.212 226.505 249.633;169.592 195.072 239.092 273.706 370.2

17、81;187.652 250.287 322.767 426.730 598.243;x=1975;y=15;z=interp2(X,Y,Z,x,y,linear)plot3(X,Y,Z)数据如下表，求通过这些点的差值曲面（提示：使用mesh（x,y,z）、surf(x,y,z)画图，在其后加shading inerp可是曲面光滑）。XY1234561121011111315216222835272031821263228254202530333220clear;clc;X=1:1:6;Y=(1:1:4);Z=12 10 11 11 13 15;16 22 28 35 27 20;18 21 26 32 28 25;20 25 30 33 32 20;x=1:0.01:6;y=(1:0.01:4);z=interp2(X,Y,Z,x,y,spline);subplot(2,1,1)mesh(x,y,z);colorbar;subplot(2,1,2)contour(x,y,z);colorbar;matlab求解微分方程解析