24小时热门版块排行榜    

查看: 1595  |  回复: 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的回帖

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

铁杆木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +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的回帖

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的回帖

yy-6133312

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by captortie at 2012-06-10 07:18:10
function exampleezplot
clear;clc
k=;
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* ...

我把x的区间取在0,0.005,y取在0,100,出来的是一条直线,但结果应该是一条向下的曲线啊,是公式不对吗?
6楼2012-06-10 08:19:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yy-6133312

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by dir51888 at 2012-06-10 07:03:51
哦,看错了,你用的@函数。

看不到图,可能数据有问题看,可以给r,tr赋值,带入实际数值分段看看。...

我把x的区间取在0,0.005,y取在0,100,出来的是一条直线,但结果应该是一条向下的曲线啊,是公式不对吗?
7楼2012-06-10 08:19:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yy-6133312

新虫 (初入文坛)

出来的图是这样


8楼2012-06-10 08:27:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

captortie

金虫 (小有名气)

【答案】应助回帖


xiegangmai: 金币+1, 谢谢参与! 2012-06-10 23:24:35
再小点试试。比如0.00001
9楼2012-06-10 21:01:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yy-6133312 的主题更新
信息提示
请填处理意见