| 4 | 1/1 | 返回列表 |
| 查看: 1560 | 回復(fù): 3 | |||
zxlzqz木蟲(chóng) (正式寫(xiě)手)
|
[求助]
關(guān)于dM/dH退磁曲線 已有1人參與
|
| 大神們有誰(shuí)知道如何得到dM/dH的圖像,是M對(duì)H的整個(gè)求導(dǎo),還是部分?M,和H都是磁滯回線就是上M和H 嗎如何說(shuō)明交換耦合 |
新蟲(chóng) (初入文壇)
|
論文J-A磁化建模理論的質(zhì)疑和修正方法里面有相關(guān)的公式推導(dǎo)。根據(jù)他的論文寫(xiě)的Matlab代碼如下: function [x,y]=lx17() u0=4*pi*(1e-7); alpha=7.092e-4; x12=linspace(0,2e4,100); % x12=linspace(3000,4000,100); x13=linspace(2e4,-2e4,200); x14=linspace(-2e4,2e4,100); M0=[0.0001,0]; % 樣本1 [x,y11]=ode45(@ode_fun11,x12,M0); [x,y12]=ode45(@ode_fun12,x13,[y11(100,1),y11(100,2)]); [x,y13]=ode45(@ode_fun11,x14,[y12(200,1),y12(200,2)]); figure; plot(x13',y12(:,2),x14',y13(:,2)); xlabel('H'); ylabel('M'); title('樣本1'); figure; plot(x13',u0*(x13'+y12(:,2)),x14',u0*(x14'+y13(:,2))); xlabel('H'); ylabel('B'); title('樣本1'); % 樣本2 [x,y21]=ode45(@ode_fun21,x12,M0); [x,y22]=ode45(@ode_fun22,x13,[y21(100,1),y21(100,2)]); [x,y23]=ode45(@ode_fun21,x14,[y22(200,1),y22(200,2)]); figure; plot(x13',y22(:,2),x14',y23(:,2)); xlabel('H'); ylabel('M'); title('樣本2'); figure; plot(x13',u0*(x13'+y22(:,2)),x14',u0*(x14'+y23(:,2))); xlabel('H'); ylabel('B'); title('樣本2'); % 樣本3; [x,y31]=ode45(@ode_fun31,x12,M0); [x,y32]=ode45(@ode_fun32,x13,[y31(100,1),y31(100,2)]); [x,y33]=ode45(@ode_fun31,x14,[y32(200,1),y32(200,2)]); figure; plot(x13',y32(:,2),x14',y33(:,2)); xlabel('H'); ylabel('M'); title('樣本3'); figure; plot(x13',u0*(x13'+y32(:,2)),x14',u0*(x14'+y33(:,2))); xlabel('H'); ylabel('B'); title('樣本3'); figure plot(x13',u0*(x13'+y12(:,2)),'g',x14',u0*(x14'+y13(:,2)),'g',x13',u0*(x13'+y22(:,2)),'r',x14',u0*(x14'+y23(:,2)),'r',x13',u0*(x13'+y32(:,2)),'b',x14',u0*(x14'+y33(:,2)),'b'); grid on % [x,Man]=ode45(@ode_fun,x12,0.0001); % % figure; % plot(x,Man(:,1)) end %樣本1磁場(chǎng)增大求解方程 function dy=ode_fun11(x,y) ms=1.5743e6; a=499; alpha=7.092e-4; k=1154.6; c=0.0198; deta=1; km=k*(1-0.96*(y(2)/ms)^2); dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)); dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c)); end %樣本1磁場(chǎng)減小求解方程 function dy=ode_fun12(x,y) ms=1.5743e6; a=499; alpha=7.092e-4; k=1154.6; c=0.0198; deta=-1; km=k*(1-0.96*(y(2)/ms)^2); dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)); dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c)); end %樣本2求磁場(chǎng)增大解方程 function dy=ode_fun21(x,y) ms=1.5755e6; a=1408.1; alpha=2.4e-3; k=2356.5; c=0.0382; deta=1; km=k*(1-0.96*(y(2)/ms)^2); dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)); dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c)); end %樣本2求磁場(chǎng)減小求解方程 function dy=ode_fun22(x,y) ms=1.5755e6; a=1408.1; alpha=2.4e-3; k=2356.5; c=0.0382; deta=-1; km=k*(1-0.96*(y(2)/ms)^2); dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)); dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c)); end %樣本3求增大解方程 function dy=ode_fun31(x,y) ms=1.5827e6; a=1017.9; alpha=1.2e-3; k=2735.8; c=0.1051; deta=1; km=k*(1-0.96*(y(2)/ms)^2); dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)); dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c)); end %樣本3求減小解方程 function dy=ode_fun32(x,y) ms=1.5827e6; a=1017.9; alpha=1.2e-3; k=2735.8; c=0.1051; deta=-1; km=k*(1-0.96*(y(2)/ms)^2); dy(1,1)=ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)/(1-alpha*ms*(-(csch((x+alpha*y(1))/a))^2/a+a/(x+alpha*y(1))^2)); dy(2,1)=(-1/km/deta*(y(1)-y(2))-c/(1-c)*dy(1,1))/(alpha/km/deta*(y(1)-y(2))-1/(1-c)); end |
新蟲(chóng) (小有名氣)
銀蟲(chóng) (小有名氣)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +6 | Tyoumou 2026-03-18 | 9/450 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +5 | 張vvvv 2026-03-15 | 7/350 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 機(jī)械專(zhuān)碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|