| 查看: 472 | 回复: 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
已经有13人回复
-
filecode
已经有11人回复
-
听说今天filecode变了
已经有26人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
固定端突然变了，今天
已经有6人回复
-
基金中了
已经有4人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
johnliu1983
至尊木虫 (著名写手)
- 应助: 43 (小学生)
- 金币: 15215.3
- 红花: 8
- 帖子: 2464
- 在线: 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
3楼2013-03-02 13:21:21
50