查看: 7931|回复: 3|关注: 0

[已答复] Matlab solve 如何解带符号变量的一元二次方程

[复制链接]

新手

5 麦片

财富积分


050


1

主题

8

帖子

0

最佳答案
发表于 2013-2-19 23:41:38 | 显示全部楼层 |阅读模式
本帖最后由 afterwrite 于 2013-2-20 11:55 编辑

需要解一个一元二次方程,未知数为q,但是方程中有很多符号变量,但之前已经进行了赋值,该方程如下

  1. eq4= 'q=mm*(Tc-Tg+(-0.000000000476692 *(((Qc - q)/Mgn*1000)^2 + 0.001329528253416*( (Qc - q)/Mgn*1000 )+ 9.16131382407062)-(-0.000000000179015* ((Qg + q)/vgm)^2 + 0.001051881416728*( (Qg +q)/vgm ) - 0.624598027137816))'
复制代码

使用solve(eq4)命令报错,(错误内容大概是这个方程不合法?)
如果把符号变量的数值用手工输进去,是一个非常简单的一元二次方程,求解没问题,但是加入了变量就无法求解。
请问我该如何用matlab解这个方程?

R14版 XP系统
刚学三四天,还请各位指点

完整代码如下:

  1. Tg=303;%空气温度
  2. Pg=111325;%空气压力
  3. Vg=3.117284;%空气流速  (标况下?)
  4. L=3  %篦冷机长度
  5. W=4  %篦冷机宽度
  6. Nx=1000  %x向网格
  7. pp=0.029* Pg /8.314/ Tg;%空气密度ρ=pp
  8. vgm=pp*Tg*Vg*L*W/273.15*101325/Pg*10^(-6) %空气质量流速
  9. Cvg=955+0.14387* Tg +3.8525*10^-5* Tg ^2+2.1036*10^-10* Tg ^3+1.2052*10^-13* Tg ^4;%空气比热
  10. Tc=1721.35; %熟料入口温度
  11. Lgc=0.5; %熟料层厚
  12. Nz=5000; %z向网格
  13. Mgc=5000 %熟料产量
  14. Mgn= Mgc/24 *1000/3600* Lgc/Nz; %熟料流速
  15. vvc=0.156;%熟料导热系数 λc=vvc
  16. dp=0.04;%熟料颗粒直径
  17. bb=0.4; %bb=ε孔隙率
  18. Cvc=699.5+0.31812* Tc -6.2308*10^-5* Tc ^2-1.3753*10^-10* Tc ^3-5.1388*10^(-14)* Tc ^4; %熟料比热
  19. Par=0.703;%普朗特常数
  20. aa=15.06*10^-6; %aa=α空气粘度系数
  21. v=19.6*10^-6;  %v=ν 运动粘度系数
  22. vv=2.59*10^(-2);%vv=λ导热系数
  23. Re=Vg * dp /aa; %雷诺数
  24. B= Par ^(1/3);
  25. C= Re ^0.5;
  26. K=1/(dp/vv/(2+1.8*B*C)+dp/2*0.25/vvc);%传热系数
  27. Qc= Cvc* Mgn* Tc/1000 %熟料总热量
  28. Qg= Cvg*pp*Vg * Tg*W/ Nx* L/ Nx; %入口废气热量
  29. mm=K / 2 * 6 * (1-bb) / dp * L / Nx * W / Nx* Lgc/ Nz;

  30. eq4= 'q=mm*(Tc-Tg+(-0.000000000476692 *(((Qc - q)/Mgn*1000)^2 + 0.001329528253416*( (Qc - q)/Mgn*1000 )+ 9.16131382407062)-(-0.000000000179015* ((Qg + q)/vgm)^2 + 0.001051881416728*( (Qg +q)/vgm ) - 0.624598027137816))'
  31. s=solve(eq4)
  32. q1=subs(s(1))
  33. q2=subs(s(2))
复制代码
回复主题 已获打赏: 0 积分

举报

新手

5 麦片

财富积分


050


1

主题

8

帖子

0

最佳答案
 楼主| 发表于 2013-2-27 11:44:16 | 显示全部楼层
沉了?沉了!沉了。。。。
回复此楼 已获打赏: 0 积分

举报

新手

10 麦片

财富积分


050


0

主题

2

帖子

0

最佳答案
发表于 6 天前 | 显示全部楼层
问题解决了吗?怎么解决的?
回复此楼 已获打赏: 0 积分

举报

新手

10 麦片

财富积分


050


0

主题

2

帖子

0

最佳答案
发表于 6 天前 | 显示全部楼层
请问怎么解决的?
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /5 下一条

快速回复 返回顶部 返回列表