请教大神 图中51定时器初值计算问题怎么做 再想请教一下 51定时器初值计算问题是不是有什么公式可以套 我该如何思考

小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
非线性方程求解,哪个大神能帮我编一下程序!
f(x)=a*ln(1+b*p*y/x)-c*ln(1+d*p*(1-y)/(1-x))=0;
其中,a=6.08565,b=0.00381;
& && &c=25.82587,d=0.00675;
& && &y=0.5;
& && &p为变量,从0、100、200、300、400、500···2000;
求相应的x的数值解,精度不高,小数点后4位就可以了,x的取值范围为0到1。
需要的是正值的实数解,这个方程在0到1之间应该是有解的,比如p=1000时,x=0.0036。
只需要x在0到1之间的解,而且在这个区间不应该存在多解的情况,请大神们帮帮忙。
你好,首先感谢你的回复,解的很不错,很用心,麻烦你了。确实像你说的,如果用matlab中的fsolve计算的话,初值的选择很重要,得一个一个试,很麻烦,而且回经常出现答案为NaN的情况。所以想请教一下你,这种情况的话matlab有没有方法可以进行更好的计算。还有,你说的lingo,以前没有接触过,它解非线性方程的时候不需要初值吗,是不是解的效果和精度都要比matlab高,你更加推荐它吗?
初值依赖是由于算法的局限性,MATLAB在这方面做得并不太突出,我平常用MATLAB也不多,太庞大而慢,不是我喜欢的类型,加上版权的缘故,:)
& &LINGO算是一款比较通用的专业优化软件吧,可以解决比较多的类型的问题,提不提供初值都可以,专业软件在算法上会额外用心的,也比较小巧,几十兆吧。
& &我觉得软件都是有适用性,而且看个人喜好吧。在你的这个问题上我是绝对推荐LINGO的。当然,MATLAB也有它的价值,我不否认,只是我不喜欢而已。:)
研究生必备与500万研究生在线互动!
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研求助大神,我是新手,不要笑话我啊,我想请教一下用sitemap x抓取后出现错误页,然后再怎么做?_百度知道小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
请教大神如何求解这个问题,最好能给个算法,把所有的金币都捐出来了。。。。
Ms=1.573e6;a=499,alpa=7.092e-6;k=1154.6;c=0.0198;dH/dt&0时,delt=1; dH/dt&0时,delt=-1;(H,M,B都是随时间t变化的),u0=48*pie-7;
Man=Ms(coth((H+alpa*M)/a)-(a/(H+alpa*M)));
B=u0*(M+H);
求H=上B的数值解,并画出B-H图。小弟不甚感激!
多谢指导,对,是求磁滞回线,最终是求B-H关系,我说B、H是时间的函数,这是为了说明dH/dt&0这个判断;经过查证后u0=4*pi*e-7(多谢提醒);Man的表达式的物理意义是修改后的朗之万函数,表示的是无磁滞时磁化曲线,由于之前参考的文献有无,应该修正为Man=Ms(coth((H+alpa*Man)/a)-(a/(H+alpa*Man)));望兄台不吝赐教,不胜感激,如此问题能圆满解决,将送上所有金币,多谢!
那这么说,Man的表达式是一个非线性方程了?
我做了一下,你给的公式似乎可行。其实就是数值积分一个常微分方程,解出M作为H的函数,然后求出B。初值始于M=0,H=0,然后多循环几圈,自然得到Hysteresis Loop。我用的是简单的predictor-corrector,没有稳定性问题。需要注意的是Langevin function L(x)=coth(x)-1/x 的求法,在0附近要用展开。结果如图,看起来数量级是对的。敢问这是什么材料的?
多谢,对,这是我参考文献上的一种含碳0.2%的钢的样本,能否把程序发一下,我后面自己的材料还需要修改参数,金币已经全部发给你了,多谢!
多谢!小弟电气学科的菜鸟,对fortran完全不了解,下午才装了一个vs2010和Intel Complier XE,现在还不知道怎么运行,还需要研究研究,不知道这个fortran能否转换成matlab代码?能麻烦大神把5楼那个Man函数形式的代码给出来吗?多谢(送了2朵鲜花以表感谢)!
刚刚运行了一下,怎么报错呢?能否加一下你的qq,向你请教一下。。。。
O5BV0VGK{BF_}RLZ4BROGON.jpg
大神,能否加一下QQ,,想请教一下你。。。。
我不用studio,你需要找人帮你把fortran 做好先。
研究生必备与500万研究生在线互动!
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研

我要回帖

更多关于 什么 请教大神 的文章

 

随机推荐