24小时热门版块排行榜    

查看: 1599  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yy-6133312

新虫 (初入文坛)

[求助] 用matlab绘制含参隐函数的问题

隐函数为
-30*h*a1*tr^2-12*h*a4*tr^5-15*h*a3*tr^4-20*h*a2*tr^3-60*h*a0*tr
+30*h*a1*tr0^2+12*h*a4*tr0^5+15*h*a3*tr0^4+20*h*a2*tr0^3+60*h*a0*tr0-15*q*r0*r0+15*q*r*r=0;
其中h=0.0003;r0=0.005;
a0=25.15;a1=0.16;a2=6.996e-004;a3=3.06e-006;a4=1.34e-008;
tr0=20;q=100;
求r和tr的关系曲线图
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dir51888

铁杆木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 谢谢参与! 2012-06-10 23:24:10
ezplot 是符号画图,你的符号定义在哪里?
在ezplot前边添加:
syms tr r

另外,检查公式合理性,画图后没有曲线,可以自己定义一下r tr的取值范围,分段看结果图。
3楼2012-06-10 07:01:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

yy-6133312

新虫 (初入文坛)

h=0.0003;r0=0.005;a0=25.15;a1=0.16;a2=6.996e-004;a3=3.06e-006;a4=1.34e-008;tr0=20;q=100;
f=@(r,tr)-30*h*a1*tr^2-12*h*a4*tr^5-15*h*a3*tr^4-20*h*a2*tr^3-60*h*a0*tr+30*h*a1*tr0^2+12*h*a4*tr0^5+15*h*a3*tr0^4+20*h*a2*tr0^3+60*h*a0*tr0-15*q*r0*r0+15*q*r*r;

ezplot(f)
为什么我这样编了以后没图形?
2楼2012-06-09 23:32:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dir51888

铁杆木虫 (小有名气)

【答案】应助回帖

引用回帖:
3楼: Originally posted by dir51888 at 2012-06-10 07:01:11
ezplot 是符号画图,你的符号定义在哪里?
在ezplot前边添加:
syms tr r

另外,检查公式合理性,画图后没有曲线,可以自己定义一下r tr的取值范围,分段看结果图。

哦,看错了,你用的@函数。

看不到图,可能数据有问题看,可以给r,tr赋值,带入实际数值分段看看。
4楼2012-06-10 07:03:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

captortie

金虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 谢谢参与! 2012-06-10 23:24:22
function exampleezplot
clear;clc
k=[0.0003,0.005,25.15,0.16,6.996*10^(-4),3.06*10^(-6),1.34*10^(-8),20,100];
ezplot(@(r,tr)symezplot(r,tr,k))
function f=symezplot(r,tr,k)
h=k(1);r0=k(2);a0=k(3);a1=k(4);a2=k(5);a3=k(6);a4=k(7);tr0=k(8);q=k(9);
f=-30*h*a1*tr.^2-12*h*a4*tr.^5-15*h*a3*tr.^4-20*h*a2*tr.^3-60*h*a0.*tr+30*h*a1*tr0^2+12*h*a4*tr0^5+15*h*a3*tr0^4+20*h*a2*tr0^3+60*h*a0*tr0-15*q*r0*r0+15*q.*r.*r;
结果显示不出来可能与你的公式或参数有关
5楼2012-06-10 07:18:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见