求救matlab绘图!!!! 15

我对matlab一窍不通,求大神帮忙写下代码。。。谢谢了万分感谢!!!后面写错了是v^2不是v... 我对matlab一窍不通,求大神帮忙写下代码。。。谢谢了

万分感谢!!!
后面写错了 是v^2 不是v
展开
 我来答
tianxiawulang
2013-12-22 · TA获得超过2.7万个赞
知道大有可为答主
回答量:4732
采纳率:89%
帮助的人:2641万
展开全部

分析

 

1、正如楼上已经试过的,这个方程无法求出解析解;

 

2、求解微分方程需要初始条件,但这个题并没给初始条件。经过试探发现,如果初值为正,则系统发散;如果初值为负,似乎不论取什么初值,最终都收敛到一个固定值(仿真结论,并未进行理论分析)。

 

3、收敛的稳态值可以按照下面的思路求出来:稳态意味着v不变,则dv/dt=0,也就是上述方程的右侧为0,解方程可以得到

>> solve('4687052/v+173703.8+14.197136*v^2') 
ans = 
[                                    -25.610148923197372792334831397939]
[ 12.805074461598686396167415698969-112.81418116868122677783012067587*i]
[ 12.805074461598686396167415698969+112.81418116868122677783012067587*i]

保留实数解,即v=25.61。

 

 

代码

V0 = [-0.001 -10 -30 -40];
lstr = cell(size(V0));
clf
for i = 1 :length(V0)
    v0 = V0(i);
    [t, v] = ode45(@(t,v)4687052/v+173703.8+14.197136*v^2,[0 0.001],v0);
    plot(t, v, 'color', rand(1,3));
    hold on
    lstr{i} = sprintf('v = %g', v0);
end
legend(lstr{:})
xlabel('t')
xlabel('v')

 

结果

东莞大凡
2024-11-19 广告
标定板认准大凡光学科技,专业生产研发厂家,专业从事光学影像测量仪,光学投影测量仪.光学三维测量仪,光学二维测量仪,光学二维测量仪,光学三维测量仪,光学二维测量仪.的研发生产销售。东莞市大凡光学科技有限公司创立于 2018 年,公司总部坐落于... 点击进入详情页
本回答由东莞大凡提供
贞皖0F
2013-12-22 · TA获得超过357个赞
知道小有建树答主
回答量:134
采纳率:0%
帮助的人:159万
展开全部
没有找到显函数,得到的是隐函数:
v=dsolve('Dv=4687052/v+173703.8+14.197136*v^2','t')
Warning: Explicit solution could not be found; implicit solution returned.
> In dsolve at 310

v =

t-Int(62500/(292940750000+10856487500*_a+887321*_a^3)*_a,_a = .. v)+C1 = 0
更多追问追答
追问
噢噢,结果不重要    我主要是想知道这个怎么画出图形!!
追答
你要画什么图形?
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
丸梓子儿
高粉答主

2020-05-08 · 每个回答都超有意思的
知道答主
回答量:9.1万
采纳率:2%
帮助的人:6805万
展开全部
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
收起 更多回答(1)
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式