24小时热门版块排行榜

返回列表

kjxywg

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 707.5
帖子: 21
在线: 31.3小时
虫号: 2061642
注册: 2012-10-15
性别: MM
专业: 化学信息学

[求助] 用matlab求解符号微分方程组时报错，哪位仁兄能帮忙解决下，小弟不胜感激

r=dsolve('Da=-k(1)*a','Db=k(2)*(((k(5)^2*k(7)^2/(k(6)+k(7))^2*b^2+4*k(2)*(k(1)*a+k(3)*b))^0.5-k(5)*k(6)/(k(6)+k(7)))/2/k(2))^2-(k(3)+k(4))*b-k(5)*b*((k(5)^2*k(7)^2/(k(6)+k(7))^2*b^2+4*k(2)*(k(1)*a+k(3)*b))^0.5-k(5)*k(6)/(k(6)+k(7)))/2/k(2)+k(6)*k(5)/(k(6)+k(7))*b*((k(5)^2*k(7)^2/(k(6)+k(7))^2*b^2+4*k(2)*(k(1)*a+k(3)*b))^0.5-k(5)*k(6)/(k(6)+k(7)))/2/k(2)-(k(9)+k(11)+k(13)+k(16))*b^2+k(10)*d+k(12)*e+k(14)*f','Dc=k(7)*k(5)/(k(6)+k(7))*b*((k(5)^2*k(7)^2/(k(6)+k(7))^2*b^2+4*k(2)*(k(1)*a+k(3)*b))^0.5-k(5)*k(6)/(k(6)+k(7)))/2/k(2)-k(8)*c','Dd=k(9)*b^2-k(10)*d','De=k(11)*b^2-k(12)*e','Df=k(13)*b^2-(k(14)+k(15))*f','Dg=k(15)*f','Dh=k(16)*b^2','Di=k(4)*b+k(8)*c','a(0)=100','b(0)=0','c(0)=0','d(0)=0','e(0)=0','f(0)=0','g(0)=0','h(0)=0','i(0)=0','t')
??? Error using ==> mupadmex
Error in MuPAD command: Unknowns must be functions [ode::new]

Error in ==> sym.sym>sym.mupadmexnout at 2018
out = mupadmex(fcn,args{:});

Error in ==> dsolve>mupadDsolve at 190
[var_list,R] = mupadmexnout('symobj::dsolve',sys,x,ignoreConstraints);

Error in ==> dsolve at 97
[R,vars] = mupadDsolve(ignoreConstraints,varargin{1:narg});

回复此楼