matlab怎么求求数值解

结果：代码：clearallclcf=@(x,y)([y(2);   0.357*y(1)-0.1905*y(1)*y(2)]);[x,Y]=ode45(f,[0100]

h=solve('1+a*x+b*x^2+c*x^3=0');a=1;b=2;c=3;eval(h)ans=-0.78390.0586+0.6495i0.0586-0.6495i或者>>subs(h)

symsa1a2c1c2%以下的值你自己改一下aa1=1;cc1=2;aa2=3;cc2=4;[x5,y5]=solve('y5=a1*x5+c1','y5=a2*x5+c2','x5','y5');

你没有给出参数的具体值,我只好编了.functionhhfsolve(@myfun,[100,10])functionf=myfun(x)T=x(1);t=x(2);c1=1;c2=9;a=8;b=3

去问你的老师啥给他打电话不久知道了

B1的值是一直不变吗?最终是要求解出变量B?参考如下程序：f = @(a,b,c,B)A1-B1   % 这里把A1,B1具体表达式代入,B

你可能误解了解析解与数值解的来源.solve给你的是解析解,你得到解析解以后,给参数赋值,然后再使用eval()函数,即可得到结果（但不是你说的数值解）.如果需要数值解,用积分来说吧,我们有辛普森(s

symx;v=int(e^(x^2/(-2)),0,1)vpa(v)

程序：clear[xy]=solve('A1=x*(y-B1)','A2=x*(y-B2)','x','y')x=-(A1-A2)/(B1-B2)y=(A1*B2-A2*B1)/(A1-A2)clea

你这个属于边值问题.得用打靶算法.如果想简单,你令一个变量为x=y-2001,就变成正常的微分方程.

首先建立vdp1.m文件functiondx=vdp1(t,x)dx=zeros(2,1);dx(1)=x(2);dx(2)=(10000*x+(96.321*exp(-t/378.47108)+28

a=1;b=.1;[t,x]=ode45(@(t,x)[ab;-ba]*x,[0,10*pi],[1;1]);plot(x(:,1),x(:,2));xlabel('x1');ylabel('x2')

clear;clc[T,Dxyz]=ode45('sdfd',[012],[0;1;1]);plot(T,Dxyz)--------------------------functiondxyz=sdf

改动：functionS=history(t)S=ones(3,1);改为S=ones(4,1);OVER!PS.多翻翻help!

只说a,b,c,d,h,N为已知量,又不给出具体值我随便取了一组值,你看看吧：clearall;clc;symsabcdhN[x,y]=solve('1+y*d*(x-1)/(6.283*b*a*(2

A = [0 1 2 0 4 5];B=  [-1 2 0 1 3&n

fun=@(x,y)y.*log10(y)-x创建关于x,y的匿名函数log10(x)由于Matlab中自然对数的符号是log(与我们数学上的不同),以10为底的对数符号是log10,同理有log2,

这个很简单,一个命令就够了：ode45(inline('y-2*y^3*sin(t)','t','y'),3,1)先用inline创建一个内联函数描述微

在Matlab下输入：edit,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%functiony=zhid

Ihopethisiswhatyouwant:figure;set(gca,'ytick',[0.1.3.9]);set(gca,'yticklabel',{'a','b','c','.9'});gr