&& 查看话题
请教大神如何求解这个问题，最好能给个算法，把所有的金币都捐出来了。。。。
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图。小弟不甚感激！
u0=48*pi*e-7; 你这是求磁滞回线吧？ 不就是求出B-H曲线吗， 为什么是时间的函数呢？ 而且应该是 \mu_0=4\pi 10^{-7} 是吗？你用的是SI单位？ 我在编程时，碰到一些问题，能否检查 你的Man 表达式？ 在M=0，H=0 时，当dt=1时，应有dM/dH &0, 可是，按你的Man 表达式，似乎dM/dH&0。再者，Man的物理意义是什么？ : Originally posted by pippi6 at
我在编程时，碰到一些问题，能否检查 你的Man 表达式？ 在M=0，H=0 时，当dt=1时，应有dM/dH &0, 可是，按你的Man 表达式，似乎dM/dH&0。再者，Man的物理意义是什么？ 多谢指导，对，是求磁滞回线，最终是求B-H关系，我说B、H是时间的函数，这是为了说明dH/dt&0这个判断；经过查证后u0=4*pi*e-7（多谢提醒）；Man的表达式的物理意义是修改后的朗之万函数，表示的是无磁滞时磁化曲线，由于之前参考的文献有无，应该修正为Man=Ms(coth((H+alpa*Man)/a)-(a/(H+alpa*Man)));望兄台不吝赐教，不胜感激，如此问题能圆满解决，将送上所有金币，多谢！ : Originally posted by zoumi89 at
多谢指导，对，是求磁滞回线，最终是求B-H关系，我说B、H是时间的函数，这是为了说明dH/dt&0这个判断；经过查证后u0=4*pi*e-7（多谢提醒）；Man的表达式的物理意义是修改后的朗之万函数，表示的是无磁滞时磁化曲 ... 那这么说，Man的表达式是一个非线性方程了？ : Originally posted by zoumi89 at
多谢指导，对，是求磁滞回线，最终是求B-H关系，我说B、H是时间的函数，这是为了说明dH/dt&0这个判断；经过查证后u0=4*pi*e-7（多谢提醒）；Man的表达式的物理意义是修改后的朗之万函数，表示的是无磁滞时磁化曲 ... 我做了一下，你给的公式似乎可行。其实就是数值积分一个常微分方程，解出M作为H的函数，然后求出B。初值始于M=0，H=0，然后多循环几圈，自然得到Hysteresis Loop。我用的是简单的predictor-corrector，没有稳定性问题。需要注意的是Langevin function L(x)=coth(x)-1/x 的求法，在0附近要用展开。结果如图，看起来数量级是对的。敢问这是什么材料的？
快照4.jpg : Originally posted by pippi6 at
我做了一下，你给的公式似乎可行。其实就是数值积分一个常微分方程，解出M作为H的函数，然后求出B。初值始于M=0，H=0，然后多循环几圈，自然得到Hysteresis Loop。我用的是简单的predictor-corrector，没有稳定性问 ... 多谢，对，这是我参考文献上的一种含碳0.2%的钢的样本，能否把程序发一下，我后面自己的材料还需要修改参数，金币已经全部发给你了，多谢！ 这是fortran程序。用的是你第一回给的公式。你如果想用5楼改进的公式，需要加一个Newton solver。你如想做，我可以帮你
subroutine&&HysteresisLoop
&&implicit none
&&integer, parameter :: nn=200,nt=8*nn
&&integer :: i
&&real*8, parameter :: pi=3.793d0 &
& && & ,cmu0=4.d-7*3.793d0 &
& && & ,ms=1.573d6,a=499d0,alpha=7.092d-6,k=,c=0.0198d0,hmax=20000d0 &
& && & ,eps=1.d-16
&&real*8 :: man,h,m,dmdh,dmandh,dt,dh,m1,lx,x,l1x
&&real*8, dimension(-nn:nn) :: mm,hh
&&real*8, dimension(0:nt) :: mt,ht,bt
&&dh=hmax/nn
&&do i=1,nt
& &&&! dt and next h
& &&&call get_dt
& &&&ht(i)=ht(i-1)+dh*dt
& &&&! predictor
& &&&m=mt(i-1)
& &&&h=ht(i-1)
& &&&call get_man
& &&&call get_dmdh
& &&&m1=mt(i-1)+(ht(i)-ht(i-1))*dmdh
& &&&! corrector
& &&&m=(mt(i-1)+m1)/2
& &&&h=(ht(i-1)+ht(i))/2
& &&&call get_man
& &&&call get_dmdh
& &&&mt(i)=mt(i-1)+(ht(i)-ht(i-1))*dmdh
& &&&bt(i)=cmu0*(mt(i)+ht(i))
& &&&if(mod(i,nt/100)==0) then
& && &&&print&(1(a,i7),5(a,f15.6))&,' i=',i,' m=',mt(i),' h=',ht(i),' b=',bt(i) &
& && && && & ,' dmdh=',dmdh,' dt=',dt
& &&&end if
&&open(20,file='output.data/HysteresisLoop.dat',status='unknown',form='unformatted')
&&write(20) nt
&&write(20) (ht(i),i=0,nt)
&&write(20) (bt(i),i=0,nt)
&&close(20)
&&subroutine get_dt
& & if(dt==1 .and. abs(ht(i-1)-hmax) & eps ) then
& && & dt=-1
& & end if
& & if(dt==-1 .and. abs(ht(i-1)+hmax) & eps ) then
& && & dt=1
& & end if
&&end subroutine get_dt
&&subroutine get_dmdh
& & implicit none
& & dmdh=((1-c)*(man-m)+c*dt*k*dmandh)/(dt*k-alpha*(1-c)*(man-m))
&&end subroutine get_dmdh
&&subroutine get_man
& & implicit none
& & x=(h+alpha*m)/a
& & call Langevin_function
& & man=ms*lx
& & dmandh=ms*l1x/a
&&end subroutine get_man
&&subroutine Langevin_function
& & if(abs(x) & 1.d-3) then
& && & lx=x/3-x**3/45+2*x**5/945d0
& && & l1x=1/3d0-x**2/15+2*x**4/189d0
& && & lx=1/tanh(x)-1/x
& && & l1x=1-1/tanh(x)**2+1/x**2
& & end if
&&end subroutine Langevin_function
end subroutine HysteresisLoop : Originally posted by pippi6 at
这是fortran程序。用的是你第一回给的公式。你如果想用5楼改进的公式，需要加一个Newton solver。你如想做，我可以帮你
subroutine&&HysteresisLoop
&&implicit none
&&integer, parameter :: nn=200,nt=8*nn
... 多谢！小弟电气学科的菜鸟，对fortran完全不了解，下午才装了一个vs2010和Intel Complier XE，现在还不知道怎么运行，还需要研究研究，不知道这个fortran能否转换成matlab代码？能麻烦大神把5楼那个Man函数形式的代码给出来吗？多谢（送了2朵鲜花以表感谢）！ : Originally posted by pippi6 at
这是fortran程序。用的是你第一回给的公式。你如果想用5楼改进的公式，需要加一个Newton solver。你如想做，我可以帮你
subroutine&&HysteresisLoop
&&implicit none
&&integer, parameter :: nn=200,nt=8*nn
... 刚刚运行了一下，怎么报错呢？能否加一下你的qq，向你请教一下。。。。
O5BV0VGK{BF_}RLZ4BROGON.jpg : Originally posted by pippi6 at
这是fortran程序。用的是你第一回给的公式。你如果想用5楼改进的公式，需要加一个Newton solver。你如想做，我可以帮你
subroutine&&HysteresisLoop
&&implicit none
&&integer, parameter :: nn=200,nt=8*nn
... 大神，能否加一下QQ，，想请教一下你。。。。 : Originally posted by zoumi89 at
大神，能否加一下QQ，，想请教一下你。。。。... 我不用studio，你需要找人帮你把fortran 做好先。这两道计算题我怎么算都算不出结果求大神解答第一道的³不是3的，是立方根那个3 后面的小3同_百度作业帮
这两道计算题我怎么算都算不出结果求大神解答第一道的³不是3的，是立方根那个3 后面的小3同
这两道计算题我怎么算都算不出结果求大神解答第一道的³不是3的，是立方根那个3 后面的小3同
1、4-3*(2/3)=22、-2+2/3+7/8=-11/48
³√1-19/27
我算出来是4-3
³√6/27
好吧我全错了T.T
第一道：4-立方根（8*27/27）=4-立方根8=2第二道：-2+2/3-立方根（343/512）=-2+2/3-7/8求大神算出这两道题 _百度作业帮
求大神算出这两道题
求大神算出这两道题&回复：求大神算一下这道题_李毅吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0可签7级以上的吧50个
本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：20,195,569贴子：
回复：求大神算一下这道题
第一!绝对不意气用事!
第二!绝对不漏判任何一件坏事!
第三!绝对裁判的公正漂亮!
裁判前来进见!
这场争夺战由我来做裁判!
我的灯呢？
卧槽我的灯呢？
有没有谁看见我的灯？
很多时候换位思考一下结果就出了，顾客赚了多少，王师傅就亏了多少，鞋子20，王师傅又退了30，也就是说王师傅共亏了50
生命本是一次旅行，在每一次停泊时都要清理自己的口袋，把更多的位置空出来，让自己活得更轻松、更自在。心无杂物，才会有满心室的暖暖阳光，才会有从容生活轻松涉世的自信和勇气。只有拥有一颗空灵的心，才能注入快乐，注入幸福。
当你老了，头白了，睡意昏沉， 炉火旁打盹，请取下这部诗歌， 慢慢读，回想你过去眼神的柔和， 回想它们昔日浓重的阴影; 多少人爱你青春欢畅的时辰， 爱慕你的美丽，假意或真心， 只有一个人爱你那朝圣者的灵魂， 爱你衰老了的脸上痛苦的皱纹; 垂下头来，在红光闪耀的炉子旁， 凄然地轻轻诉说那爱情的消逝， 在头顶的山上它缓缓踱着步子， 在一群星星中间隐藏着脸庞。
王用假的50去邻居家换了真的50，然后给了客户30。手里还有20。赔给邻居50，就是赔了30，再加上鞋的成本30。加起来60。
如果钱是真的，亏了10块，后来发现钱是假的，又少了50身上，所以不是60是多少？不服来辩
向邻居换钱，等于向邻居借了50又还了可以不计算。进鞋花了30，又找出30，故损失60.
哭泣的刀前来顶帖哭泣的刀前来顶帖（伪面君纸）
60　　　--来自助手版贴吧客户端 看而不言，笑而不语，我是来送小尾巴的。。。
反正我算是亏70，楼主你就没有答案么
答案是亏出翔
留名“最锋利的剑的刀锋也无法与一颗和平的心的平静相匹敌。”——无极剑圣，易大师
假钱-50 进价-30 卖+20。结论60。此贴完结。
♚♚♚♚♚♚♚♚♚我种下一个蘑菇终于长成大蘑菇今天是个伟大日子把它种在草丛里把它种在河道里让它为我指明道路变身斥候牺牲自己只为探探路把我一切都献给你团战不能输我让每次carry都会变得有意义团战虽短决定成败不能放弃你是我的小呀小提莫怎么爱你都不嫌多萌萌的身影走遍每个角落点亮那地图的火 火火火火你是我的小呀小提莫怎么卖萌都不嫌多你的蘑菇种在我心的角落种下希望就会收获从来不觉得厌倦每天都在第一线每天游戏都很新鲜你的蘑菇开视线有你的Q致盲点你混助攻我拿人头中路carry卖掉自己换点助攻钱上路隐身蹲点逼掉敌方的闪现下路和我一起双杀打爆了对面小龙大龙埋伏一波 基情温暖你是我的小呀小提莫怎么爱你都不嫌多萌萌的身影走遍每个角落点亮那地图的火 火火火火你是我的小呀小提莫怎么卖萌都不嫌多你的蘑菇种在我心的角落种下希望就会收获你是我的小呀小提莫怎么爱你都不嫌多萌萌的身影走遍每个角落点亮那地图的火 火火火火你是我的小呀小提莫怎么卖萌都不嫌多你的蘑菇种在我心的角落种下希望就会收获
如果顾客买1双鞋是亏60元，如果买2双是亏70元钱！没错了这就是正确答案！
没有算的意义，不符合生活常识，而且问题还钱时说的很模糊
八十不多说，一双鞋子30，假钱50。
内&&容:使用签名档&&
保存至快速回贴这道题，我算的是90，而出题人说得正确答案是60，我不解，求大神赐教 _百度作业帮
这道题，我算的是90，而出题人说得正确答案是60，我不解，求大神赐教
这道题，我算的是90，而出题人说得正确答案是60，我不解，求大神赐教&
王师傅赔了一双鞋，30元，赔了邻居50元，手里剩了20元，所以一共赔了60元