EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB中常微分方程的求解、利用MATLAB求极值。
+ m. F7 t* |' M8 `& z9 V5 c 例题如下:
% X a% D9 o) ^8 h
程序如下: 1.syms x syms y dsolve('(x^2-1)*Dy+2*x*y-sin(x)') 2.(1) P=0:0.01:10000; R=P.*(200./(P+4)-18); [Rmax,Pmax]=max(R) (2) X=0:0.01:10000; R=10.*X.*exp(-(X.^2)); [Rmax,Xmax]=max(R) C=(108+5.*X+X.^(3./2))./100; [Cmin,Xmin]=min(C) W=10.*X.*exp(-(X.^2))-(108+5.*X+X.^(3./2))./100; [Wmax,Xmax]=max(W) 3. syms x disp('Ô-oˉêyf') f=sin(x).*exp(-(x.^2)) disp('ÏóoˉêyF') F=fourier(f) 4. syms x f=sin(x).*exp(-(x.^2)); Fs=laplace(f,x) 5. syms x f=sin(x.*sin(x)) p50=taylor(f,'Order',5,'ExpansionPoint',0) ezplot(f) 6. x=linspace(0,4); solinit =bvpinit(linspace(0,4,5),[0 1]); sol=bvp4c(@twoode,@twobc,solinit); y=deval(sol,x); y2=sin(pi*x); plot(x,y(1,:),x,y2(1,:))
. P& x" n" X& G$ R; `$ h4 o |