nbcj.net
当前位置:首页 >> mAtlAB解方程+画图 >>

mAtlAB解方程+画图

x=10:0.01:10; y=x.*sin(x)-1/2; plot(x,y) z=@(x)x.*sin(x)-1/2; a=fzero(z,0) a=-0.740840955095491 a就是函数在0附近的零点。求其他值就把0换成其他数 这是个偶函数,所以知道半个轴就知道另一半了。

问题出在符号量和数值量的的关系上。 建议修改 syms kx ky k3x k5x k2y k4y S=solve('k3x^2=(2*pi/1.55)^2*(1.4949^2-1.465^2)-kx^2',... 'k5x^2=(2*pi/1.55)^2*(1.4949^2-1.465^2)-kx^2',... 'kx*6=pi+atan(k3x/kx)+atan(k5x/kx)',... 'k2y^2=(...

代码供参考: n = 1.5;% x1 = φ, x2 = Rdx = @(z,x)[1/(x(2)*(n+(n+1)*(2*cos(x(1))+cos(x(1))*cos(x(1))))); -1/tan(x(1))];[t,x]=ode45(dx,[0 1],[0.5 0.5]);plot(t,x(:,2)) 你现在的问题有很多不清楚的地方: 1、第一个表达式很乱,我不确定...

dsolve()求解微分方程后,使用ezplot()作图,参考代码: >> y = dsolve('Dy=y*cos(x)/(1+2*y^2)','y(0)=1','x');>> ezplot(y,[-10,10]) % [-10,10]为x区间

syms t v=dsolve('Dv=(190.708-90.64*v^2)/47.27','v(0)=0','t'); t=0:0.00001:0.002; v=eval(v); plot(t,v) 使用这样的方法求解,但从结果看好像你的方程有问题!

syms x f=(x-1)/x-2*x/(x-1)-1; ezplot(f) x=solve(f)

说明 求解常微分方程一般使用ode系列函数(其中最常用的是ode45),基本做法就是把微分方程用函数表示出来,然后调用ode函数求解即可。对于本题而言,使用匿名函数描述微分方程更为方便一些。 参考代码 k0=8; k1=1; K=1; k2=5; n=2; % 常数定义d...

你好!给你变了一段代码: 主程序 clc;clear allglobal tA=1:1:50;x=zeros(size(A));y=x;for i=1:length(A) x0=1; t=A(i); y(i)=fsolve(@(x)f(x),x0); x(i)=216.77-(0.4*(428.15-y(i)))/(0.9*(exp(-t*51.84/904.5)-1));endfigure(1)plot(A,x,'r-...

由于方程比较复杂,解析解不能用初等函数表示 只是要获得图像的话,用数值计算的办法可能更方便 fun=@(t,y) 1.44*(10^9)*(1-y).*exp(-109170./(8.314*t)); [T,Y] = ode23t(fun,500:600,0); DY=fun(T,Y); AX=plotyy(T,Y,T,DY); set(get(AX(1),'Xl...

该二阶微分方程可以ode()函数求得,其方程的数值解,然后用plot()函数绘制其x—y、x—dy/dx函数图形,以及y—dy/dx函数的相平面图形。 主要执行代码

网站首页 | 网站地图
All rights reserved Powered by www.nbcj.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com