function  f= solveMagneticDipoleEquation(x)
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
%ÒÑÖªÊý
a=-150;
b=150;
Bx=-48030;
By=29370;
Bz=-19077.3;
%δ֪Êý
Bt=x(1);
c=x(2);
m=x(3);
n=x(4);
p=x(5);
Rl=sqrt(a.^2+b.^2+c.^2);
f(1)=Bt*(3*(m*a+n*b+p*c)*a/Rl.^5-m/Rl.^3)/Bx-1;
f(2)=Bt*(3*(m*a+n*b+p*c)*b/Rl.^5-n/Rl.^3)/By-1;
f(3)=Bt*(3*(m*a+n*b+p*c)*c/Rl.^5-p/Rl.^3)/Bz-1;
f(4)=m.^2+n.^2+p.^2-1;
end

Parameters Bt[10^10,10^12];
Parameters c[160,200];
Parameters m[-0.55,0.54];
Parameters n[-0.52,0.54];
Parameters p[0.65,1];
Parameters Rl;
Constant a=-150;
Constant b=150;
Constant Bx=-48030; //Bx,By,BzÊDZä²ÎÊý£¬´ÓÎļþÖжÁ³öÒ»ÐÐÔÙ½âÒ»´Î·½³Ì£¬¶ÔÓÚÿһ×éBx,By,BzÐèÒªÇóÒ»´Î·½³Ì×é
Constant By=29370;
Constant Bz=-19077.3;
//BatchFileModel;
//VarConstant Bx="Sheet1[A1:A240]";
//VarConstant By="Sheet1[B1:B240]";
//VarConstant Bz="Sheet1[C1:C240]";
//DataFile "D:\Mag1230.xls[Sheet1[A1:A240]]";
//RowDataSet;
     //Bx="Sheet1[A1:A240]";
     //By="Sheet1[B1:B240]";
     //Bz="Sheet1[C1:C240]";
//EndRowDataSet;
Function Bt*(3*(m*a+n*b+p*c)*a/Rl^5-m/Rl^3)/Bx-1=0;
Bt*(3*(m*a+n*b+p*c)*b/Rl^5-n/Rl^3)/By-1=0;
Bt*(3*(m*a+n*b+p*c)*c/Rl^5-p/Rl^3)/Bz-1=0;
m^2+n^2+p^2-1=0;
Rl-sqrt(a^2+b^2+c^2)=0;

x0=[999999999999 161.9 0.401 -0.0465 0.9148 266.856];
options=optimset('Display','iter','MaxFunEvals',10^5,'TolFun',1e-18,'TolX',1e-18,'MaxIter',1e4,'Algorithm','Levenberg-Marquardt');
[x val]= fsolve(@solveMagneticDipoleEquation,x0,options);

3Â¥: Originally posted by mjjackey at 2015-03-03 11:11:16
ÎÒÒ²ÏëÍêÈ«ÓÃ1stopt£¬¹Ø¼üÊÇÎÒ²»Ö¹Òª½â³ö·½³Ì×éµÄ½â£¬»¹Òª£º¡°Òª¼ÌÐøÒ»²½×÷·ÖÎö£¬Bx£¬ By£¬ BzÊDzâÊÔÊý¾Ý£¬Òª´Ó±í¸ñÖжÁ³öµÄ£¬¶ÔӦÿһ×éÊý¾ÝÇóÒ»´Î·½³Ì×飬Ï൱Óڱ䳣Á¿£¬°´ÕÕÉÏÃæ×¢ÊÍÖвÉÓõÄÁ½ÖÖ·½·¨£¬¶¼Ã»Óдï ...

Constant a=-150, b=150, Bx=-48030, By=29370,Bz=-19077.3 ;
Parameters Bt[10^10,10^12],c[160,200],m[-0.55,0.54],n[-0.52,0.54],p[0.65,1]   ;

Function Bt*(3*(m*a+n*b+p*c)*a/(sqrt(a^2+b^2+c^2))^5-m/(sqrt(a^2+b^2+c^2))^3)/Bx-1=0;
Bt*(3*(m*a+n*b+p*c)*b/(sqrt(a^2+b^2+c^2))^5-n/(sqrt(a^2+b^2+c^2))^3)/By-1=0;
Bt*(3*(m*a+n*b+p*c)*c/(sqrt(a^2+b^2+c^2))^5-p/(sqrt(a^2+b^2+c^2))^3)/Bz-1=0;
m^2+n^2+p^2-1=0;

function solve_nonl_eqs
clear all;clc
format long

x0=[1e12 160 0.4 -0.04  0.9];
options=optimset('Display','iter','Algorithm','Levenberg-Marquardt','TolFun',1e-21,'TolX',1e-21);
[x fval]=fsolve(@objfun,x0,options')


function  f=objfun(x)

%ÒÑÖªÊý
a=-150;
b=150;
Bx=-48030;
By=29370;
Bz=-19077.3;
%δ֪Êý
Bt=x(1);
c=x(2);
m=x(3);
n=x(4);
p=x(5);
Rl=sqrt(a.^2+b.^2+c.^2);



f(1)=Bt*(3*(m*a+n*b+p*c)*a/Rl.^5-m/Rl.^3)/Bx-1;
f(2)=Bt*(3*(m*a+n*b+p*c)*b/Rl.^5-n/Rl.^3)/By-1;
f(3)=Bt*(3*(m*a+n*b+p*c)*c/Rl.^5-p/Rl.^3)/Bz-1;
f(4)=m.^2+n.^2+p.^2-1;

2Â¥: Originally posted by ÔÂÖ»À¶ at 2015-03-02 17:15:39
Ô­·½³Ì×éδ֪ÊýÖ®¼äµÄÊýÁ¿¼¶Ïà²î¹ý´ó£¬ÔÚ²»¶ÔÔ­·½³Ì½øÐд¦Àí£¨ÖîÈçÈ¡¶ÔÊý½µÊýÁ¿¼¶Ö®ÀࣩµÄÇé¿öÏ£¬Çó½âµÄ½á¹û²»Îȶ¨ÊÇÕý³£µÄ£¬¶øÇÒ´Ëʱ·½³ÌÓÃMALTABÇó½â±È½ÏÀ§ÄÑ£¬¾ÍÊdzõÖµÄÑÒÔÑ¡È¡¡£½¨ÒéÓÃ1stoptÈí¼þÇó½â¼´¿É¡£
Â¥ ...

4Â¥: Originally posted by ÔÂÖ»À¶ at 2015-03-03 11:13:29
1stoptÈí¼þ¼°Æä½Ì³Ì£º
http://muchong.com/bbs/viewthread.php?tid=4586018&fpage=1&target=blank...