EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB中常微分方程的求解、利用MATLAB求极值。5 c7 V4 l2 J+ o u! ^" `5 A: \4 \( @7 ?8 K
例题如下:
8 H$ M7 F7 [) r
程序如下: 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,:)) 5 l( j& p$ d, ` s5 H7 g
|