EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB中常微分方程的求解、利用MATLAB求极值。
4 S# Q; i2 | h# w6 i) P 例题如下:
9 w! L; g& p% R O) u2 d程序如下: 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,:))
2 U$ C m# O; s" u4 ?& J4 |' ? |