matlab解方程时返回[ empty sym ]f=3.125;w=8.2;syms omega;gp=(((100*(omega-5).*f)./(((omega-5).^2-w^2).^2+(2*(omega-5).*f).^2))+0.55)/1.6;a=(w^2-2*f^2)/2;b=sqrt(w^4+4*f^4-4*f^2)/2;m=sqrt(a+b);n=sqrt(b-a);c=125*f/(16*m*n*(n^2+m^2));syms T;d1=-n*(m

来源:学生作业帮助网 编辑:作业帮 时间:2024/04/28 14:43:49
matlab解方程时返回[ empty sym ]f=3.125;w=8.2;syms omega;gp=(((100*(omega-5).*f)./(((omega-5).^2-w^2).^2+(2*(omega-5).*f).^2))+0.55)/1.6;a=(w^2-2*f^2)/2;b=sqrt(w^4+4*f^4-4*f^2)/2;m=sqrt(a+b);n=sqrt(b-a);c=125*f/(16*m*n*(n^2+m^2));syms T;d1=-n*(m

matlab解方程时返回[ empty sym ]f=3.125;w=8.2;syms omega;gp=(((100*(omega-5).*f)./(((omega-5).^2-w^2).^2+(2*(omega-5).*f).^2))+0.55)/1.6;a=(w^2-2*f^2)/2;b=sqrt(w^4+4*f^4-4*f^2)/2;m=sqrt(a+b);n=sqrt(b-a);c=125*f/(16*m*n*(n^2+m^2));syms T;d1=-n*(m
matlab解方程时返回[ empty sym ]
f=3.125;
w=8.2;
syms omega;
gp=(((100*(omega-5).*f)./(((omega-5).^2-w^2).^2+(2*(omega-5).*f).^2))+0.55)/1.6;
a=(w^2-2*f^2)/2;
b=sqrt(w^4+4*f^4-4*f^2)/2;
m=sqrt(a+b);
n=sqrt(b-a);
c=125*f/(16*m*n*(n^2+m^2));
syms T;
d1=-n*(m^2+n^2)*exp(-n*T)*cos((5-m+omega)*T);
d2=(m^2+n^2)*(5-m+omega)*exp(-n*T)*sin((5-m+omega)*T);
d3=n*(m^2+n^2);
dm=2*(n^2+(5-m+omega)^2);
d=(d1+d2+d3)/dm;
e1=-n*(m^2+n^2)*exp(-n*T)*cos((5-m-omega)*T);
e2=(m^2+n^2)*(5-m-omega)*exp(-n*T)*sin((5-m-omega)*T);
e3=n*(m^2+n^2);
em=2*(n^2+(5-m-omega)^2);
e=(e1+e2+e3)/em;
g1=(n*(m^2-n^2)-2*m*n*(5+m+omega))*exp(-n*T)*cos((5+m+omega)*T);
g2=(-(m^2-n^2)*(5+m+omega)-2*m*n^2)*exp(-n*T)*sin((5+m+omega)*T);
g3=-n*(m^2-n^2)+2*m*n*(5+m+omega);
gm=2*(n^2+(5+m+omega)^2);
g=(g1+g2+g3)/gm;
k1=(n*(m^2-n^2)-2*m*n*(5+m-omega))*exp(-n*T)*cos((5+m-omega)*T);
k2=(-(m^2-n^2)*(5+m-omega)-2*m*n^2)*exp(-n*T)*sin((5+m-omega)*T);
k3=-n*(m^2-n^2)+2*m*n*(5+m-omega);
km=2*(n^2+(5+m-omega)^2);
k=(k1+k2+k3)/km;
gpt=d+e-g-k;
solve('gp-gpt',omega)
Warning:Explicit solution could not be found.
> In solve at 98
ans =
[ empty sym ]
上面是程序,gp是关于omega的函数,gpt是关于omega和T的函数,希望能从方程gp=gpt中接触omega和T的关系,可是总出错,

matlab解方程时返回[ empty sym ]f=3.125;w=8.2;syms omega;gp=(((100*(omega-5).*f)./(((omega-5).^2-w^2).^2+(2*(omega-5).*f).^2))+0.55)/1.6;a=(w^2-2*f^2)/2;b=sqrt(w^4+4*f^4-4*f^2)/2;m=sqrt(a+b);n=sqrt(b-a);c=125*f/(16*m*n*(n^2+m^2));syms T;d1=-n*(m
你的方程这么复杂,matlab在符号计算没有mathmatic等强(可能这个软件业解不出来),它求出来显示解(机器不是万能的),像这么复杂的方程,给你点建议:
手工加电脑简化方程,再试试.
用数值算法(fsolve)计算.不能求出omega与T的关系.(推荐)

  假设你用数值算法,可以把这些方程保存为一个函数,给omega和T一个任意初始值,求最后的方程eq=gp-gpt为0的解.

  例如:将下面的代码保存为testt.m函数.
  function eq=testt(OT)
  omega=OT(1);
  T=OT(2);
  f=3.125;
  w=8.2;
  % syms omega;
  gp=(((100*(omega-5).*f)./(((omega-5).^2-w^2).^2+(2*(omega-5).*f).^2))+0.55)/1.6;
  a=(w^2-2*f^2)/2;
  b=sqrt(w^4+4*f^4-4*f^2)/2;
  m=sqrt(a+b);
  n=sqrt(b-a);
  c=125*f/(16*m*n*(n^2+m^2));
  % syms T;
  d1=-n*(m^2+n^2)*exp(-n*T)*cos((5-m+omega)*T);
  d2=(m^2+n^2)*(5-m+omega)*exp(-n*T)*sin((5-m+omega)*T);
  d3=n*(m^2+n^2);
  dm=2*(n^2+(5-m+omega)^2);
  d=(d1+d2+d3)/dm;
  e1=-n*(m^2+n^2)*exp(-n*T)*cos((5-m-omega)*T);
  e2=(m^2+n^2)*(5-m-omega)*exp(-n*T)*sin((5-m-omega)*T);
  e3=n*(m^2+n^2);
  em=2*(n^2+(5-m-omega)^2);
  e=(e1+e2+e3)/em;
  g1=(n*(m^2-n^2)-2*m*n*(5+m+omega))*exp(-n*T)*cos((5+m+omega)*T);
  g2=(-(m^2-n^2)*(5+m+omega)-2*m*n^2)*exp(-n*T)*sin((5+m+omega)*T);
  g3=-n*(m^2-n^2)+2*m*n*(5+m+omega);
  gm=2*(n^2+(5+m+omega)^2);
  g=(g1+g2+g3)/gm;
  k1=(n*(m^2-n^2)-2*m*n*(5+m-omega))*exp(-n*T)*cos((5+m-omega)*T);
  k2=(-(m^2-n^2)*(5+m-omega)-2*m*n^2)*exp(-n*T)*sin((5+m-omega)*T);
  k3=-n*(m^2-n^2)+2*m*n*(5+m-omega);
  km=2*(n^2+(5+m-omega)^2);
  k=(k1+k2+k3)/km;
  gpt=d+e-g-k;
  eq=gp-gpt;


  然后在命令窗口中输入下面的代码:
  clear,clc
  x0=[2   2];
  [xj,fval]=fsolve(@testt,x0)
  求解结果如下:
  xj =  -56.2901   42.6107;
  fval = -2.5612e-010;
  fval 已经接近于0,说明求解已经收敛,
  这只是另外一种变通的方法,有时matlab无法求解符号问题,只能采用这种方法.这种方法对更加复杂的模型,也同样适应.

请教用matlab求方程符号解的问题我现在要用matlab求一个方程的符号解casein=solve('x^(0.12/0.88)+x+a','x')执行后程序返回Warning: Explicit solution could not be found.> In solve at 140 casein = [ empty sym ]但是如果把 matlab解方程时返回[ empty sym ]f=3.125;w=8.2;syms omega;gp=(((100*(omega-5).*f)./(((omega-5).^2-w^2).^2+(2*(omega-5).*f).^2))+0.55)/1.6;a=(w^2-2*f^2)/2;b=sqrt(w^4+4*f^4-4*f^2)/2;m=sqrt(a+b);n=sqrt(b-a);c=125*f/(16*m*n*(n^2+m^2));syms T;d1=-n*(m matlab怎么解方程 用matlab解方程求matlab代码. matlab解不定方程 2x+2y=xy+7x=-2:6;y=-2:6;[x,y]=meshgrid(x,y);z=2*x+2*y;k=find(z==x*y+7);x(k)',y(k)'运行结果为ans =Empty matrix:1-by-0ans =Empty matrix:1-by-0请问怎么改 用matlab怎么解方程 怎样用matlab解积分方程, 用matlab 解超越方程~ matlab解非线性方程问题 empty empty 用matlab解题出现y = [ empty sym ] 为什么 用matlab解方程时,怎么设置求解区间和结果类型?对于一个方程f(x)=0,用solve命令会求出它所有的理论解.怎样设置求解区间和解类型?比如,只需要在[0,1000]这个区间求解,而且解不是复数,只返回实 matlab里用solve函数解出来的方程[z]=solve('z-2=0')的解是sym类型的,怎么能把他变成图2的类型的变量呢?我是想用exlink把结果返回excel里,可发现sym这样的是无法返回的。而田字格图标的能返回去。 matlab求解积分总是出现empty sym的情况.解这个方程,假设其解为alfa.z=((x/(alfa*x^2+1-alfa))^2-1)^(-1/2),beta=(sin(theta)+(sin(theta)^2+4*alfa*(alfa-1))^(1/2))/(2*alfa).z和beta都是关于alfa的式子,H,R,theta都是已知.我用in matlab 解超越方程先解出方程,得到超越方程,再用MATLAB解,如图 在用matlab解矩阵方程时,什么是超静定方程 如何用matlab解一个二元方程?