24小时热门版块排行榜    

查看: 1971  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +37 医学老男孩 2026-08-13 87/4350 2026-08-18 06:35 by 医学老男孩
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
请填处理意见