| 查看: 476 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
matlab编程
|
||
|
T=1200:10:1400; a=exp(-(2.98.*log(T)+4.*T/100000+12500./T.^2+47787./T-52.01)./1.987); b=exp(-(5.79.*log(T)+1.75.*T/10000+48150./T.^2+92542./T-116.93)./1.987); c=exp(-(-0.42.*log(T)+1.1.*T/10000-89000./T.^2+114925./T-76.16)./1.987); p=[1 0 -b -2.*(a.*b-c.^2)]; A=roots(p) 求出根后，我想画出其中任意一个根随T变化的图像，该怎么编程 |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
3楼2013-03-02 13:21:21
johnliu1983
至尊木虫 (著名写手)
- 应助: 43 (小学生)
- 金币: 15243.3
- 红花: 8
- 帖子: 2468
- 在线: 150.7小时
- 虫号: 370817
- 注册: 2007-05-14
- 性别: GG
- 专业: 物理学II
【答案】应助回帖
★
感谢参与，应助指数 +1
jjdg: 金币+1, 辛苦了 2013-03-02 14:12:10
感谢参与，应助指数 +1
jjdg: 金币+1, 辛苦了 2013-03-02 14:12:10
|
Tvalue=1200:10:1400 for k=1:length(Tvalue) T=Tvalue(k); a=exp(-(2.98.*log(T)+4.*T/100000+12500./T.^2+47787./T-52.01)./1.987); b=exp(-(5.79.*log(T)+1.75.*T/10000+48150./T.^2+92542./T-116.93)./1.987); c=exp(-(-0.42.*log(T)+1.1.*T/10000-89000./T.^2+114925./T-76.16)./1.987); p=[1 0 -b -2.*(a.*b-c.^2)]; A=roots(p); target(k,1)=A(1); target(k,2)=A(2); target(k,3)=A(3); end figure(1) plot(Tvalue',target(:,1),'r') figure(2) plot(Tvalue',target(:,2),'g') figure(3) plot(Tvalue',target(:,3),'k') 图一到三是三个解随T的变化图像 |
2楼2013-03-02 01:20:21
20