| ²é¿´: 2269 | »Ø¸´: 6 | ||||
СÁõ8885ľ³æ (ÕýʽдÊÖ)
|
[ÇóÖú]
Matalab ³ÌÐòµ÷ÊÔ
|
|
ÏÂÊö³ÌÐòÔËÐгö´í£¬Çó¸ßÊÖÖ¸µã¡£ÕâÊÇÒ»¸öMÎļþ¡£ function f4=func_main(name,a,q1,q2,x2,flag) % flag=0 represent complex permittivity % flag=1 represent complex permeability data=load(name); xf=data(:,1); if (flag==0) y1=data(:,2)+data(:,3)*i; elseif (flag==1) y1=data(:,4)+data(:,5)*i; else disp('0 is complex perimittivity, 1 is complex permeability'); return; end x11=func1(a,q1,y1,x2); %µ÷Óú¯Êýfunc1£¬Çó½âx1 %Hanai formula µÄÇó½â½á¹û %Çó½âµÄÊÇx1£¬Îª¾ø¶Ô½éµç³£Êý x1=x11; a=real(x1); b=imag(x1); c=[a,b]; x0=zeros(length(data),1); F=fsolve(@(y2)func2(a,q2,x1,x2,y2),x0); %µ÷ÓÃǰÊöfunc1ËùÇóµÃµÄx1£¬´úÈ뺯Êýfunc2 %Çó½â³¬Ô½·½³ÌµÄy2Öµ¡£ f4=F; r = [real(f4), imag(f4)]; plotyy(xf,a,xf,b); save('data1.txt', 'c', '-ascii'); save('data2.txt', 'r', '-ascii'); function x1=func1(a,q1,y1,x2) x1=(q1.*x2-x2+y1.*(x2./y1).^a)./(q1+(x2./y1).^a-1); %Çó½âx1 function F=func2(a,q2,x1,x2,y2) F=(((1-q2)*(x2-x1))./(x2^a)).*(y2.^a)-y2+x1; %Çó½ây2 ÁíÍâµÄµ÷ÓôúÂëÈçÏ£º % Exop4 clear up; clc; a=1/3; q1=0.06; q2=0.107; x2=3.7+0.703*i; flag=0; func_main('3.txt',a,q1,q2,x2,flag); ¸´ÖƵ½command´°¿Ú¾Í¿ÉÒÔÔËÐС£ |
» ±¾Ìû¸½¼þ×ÊÔ´Áбí
-
»¶Ó¼à¶½ºÍ·´À¡£ºÐ¡Ä¾³æ½öÌṩ½»Á÷ƽ̨£¬²»¶Ô¸ÃÄÚÈݸºÔð¡£
±¾ÄÚÈÝÓÉÓû§×ÔÖ÷·¢²¼£¬Èç¹ûÆäÄÚÈÝÉæ¼°µ½ÖªÊ¶²úȨÎÊÌ⣬ÆäÔðÈÎÔÚÓÚÓû§±¾ÈË£¬Èç¶Ô°æÈ¨ÓÐÒìÒ飬ÇëÁªÏµÓÊÏ䣺xiaomuchong@tal.com - ¸½¼þ 1 : func_main.m
- ¸½¼þ 2 : 3.txt
- ¸½¼þ 3 : Exop4.m
2012-05-27 11:38:14, 945 bytes
2012-05-27 11:38:25, 3.1 K
2012-05-27 11:39:11, 119 bytes
» ²ÂÄãϲ»¶
070300»¯Ñ§319Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
08¹¤¿Æ 320×Ü·Ö Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´
Ò»Ö¾Ô¸Ìì½ò´óѧ»¯Ñ§¹¤ÒÕרҵ£¨081702£©315·ÖÇóµ÷¼Á
ÒѾÓÐ12È˻ظ´
307Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´
265Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´
301Çóµ÷¼Á
ÒѾÓÐ10È˻ظ´
303Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
²ÄÁϹ¤³Ì£¨×¨£©Ò»Ö¾Ô¸985 ³õÊÔ335Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
299Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´
083200ѧ˶321·ÖÒ»Ö¾Ô¸ôßÄÏ´óѧÇóµ÷¼Á
ÒѾÓÐ3È˻ظ´
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
ÎÞ·¨ÕÒµ½µ÷ÊÔÐÅÏ¢
ÒѾÓÐ3È˻ظ´
ÇóÖúһϳÌÐòµ÷ÊÔ¹ý³Ì³öÏֵĴíÎó
ÒѾÓÐ4È˻ظ´
ÂÝÐý²¨²¨Í· ³ÌÐòµ÷ÊÔ ³ýÁã´íÎó
ÒѾÓÐ12È˻ظ´
matlabµ÷ÊÔÎÊÌâ
ÒѾÓÐ3È˻ظ´
GDBµ¥²½µ÷ÊÔÎÊÌâ
ÒѾÓÐ8È˻ظ´
¶Ô"ÍÆ°åʽËíµÀµç×诡°µ÷ÊÔ¾«Í¨µÄÇë½ø£¡£¡
ÒѾÓÐ31È˻ظ´
¡¾ÇóÖú¡¿ÓÃC#±àдµÄWindowsÓ¦ÓóÌÐòÖУ¬Óõ½ÁËÊý×飬²»Äܳɹ¦µ÷ÊÔ
ÒѾÓÐ4È˻ظ´
¡¾ÇóÖú¡¿ÑáÑõµ÷ÊÔ£¿
ÒѾÓÐ8È˻ظ´
¡¾ÇóÖú¡¿VS2005ÊýÖµ¼ÆËã³ÌÐòµ÷ÊÔ¡¾¼±Ç󣡣¡£¡£¡¼Ó½ð±Ò¡¿
ÒѾÓÐ9È˻ظ´
¡¾ÇóÖú¡¿VC++ÊýÖµËã·¨¼¯Ó¦Óõ½MFCÖб¨´í£¬³ÌÐòµ÷ÊÔÇóÖú
ÒѾÓÐ8È˻ظ´
¡¾ÇóÖú¡¿VC++³ÌÐòËã·¨¼¯µ÷ÊÔ±¨´í
ÒѾÓÐ8È˻ظ´
¡¾ÇóÖú¡¿ÇëÎÊÔõôÓÃMATALAB±àдһ¸ö¿Õ¼ä¼¸ºÎÌåµÄÊÜÁ¦±äÐγÌÐò£¿£¿
ÒѾÓÐ3È˻ظ´
¡¾ÇóÖú¡¿KRUSS K10±íÃæÕÅÁ¦Òǵ÷ÊÔ·½·¨¡¾ÓÐЧÆÚÖÁ2010Äê7ÔÂ11ÈÕ¡¿
ÒѾÓÐ3È˻ظ´
¡¾ÌÖÂÛ¡¿¹ØÓÚÕë¶ÔCASTEPÊÕÁ²ÎÊÌâµÄ²ÎÊýµ÷ÊÔ×ܽᡣ
ÒѾÓÐ79È˻ظ´
csgt0
ÈÙÓþ°æÖ÷ (ÖøÃûдÊÖ)
²ÊÉ«¹Òͼ
-

ר¼Ò¾Ñé: +2 - Ó¦Öú: 367 (˶ʿ)
- ¹ó±ö: 0.457
- ½ð±Ò: 7807
- É¢½ð: 8897
- ºì»¨: 20
- ɳ·¢: 7
- Ìû×Ó: 2809
- ÔÚÏß: 515.8Сʱ
- ³æºÅ: 449046
- ×¢²á: 2007-11-02
- רҵ: ½ðÊô·Ç¾§Ì¬¡¢×¼¾§ºÍÄÉÃ×¾§
- ¹ÜϽ: ¼ÆËãÄ£Äâ

2Â¥2012-05-28 09:10:18
lijie169
ͳæ (ÖøÃûдÊÖ)
- Ó¦Öú: 159 (¸ßÖÐÉú)
- ½ð±Ò: 161.5
- É¢½ð: 2018
- ºì»¨: 7
- Ìû×Ó: 1338
- ÔÚÏß: 491.5Сʱ
- ³æºÅ: 1180408
- ×¢²á: 2010-12-30
- רҵ: ·ÖÀë¹ý³Ì
3Â¥2012-05-28 13:11:52
СÁõ8885
ľ³æ (ÕýʽдÊÖ)
- Ó¦Öú: 3 (Ó×¶ùÔ°)
- ½ð±Ò: 10042.8
- É¢½ð: 200
- Ìû×Ó: 468
- ÔÚÏß: 160.5Сʱ
- ³æºÅ: 735175
- ×¢²á: 2009-03-30
- ÐÔ±ð: MM
- רҵ: ²ÄÁÏÎïÀí»¯Ñ§
|
Äú¿ÉÒÔ°ïÎÒ¿´¿´£¬¸ÄÕýµã³ËÖ®ºó£¬ÔËÐÐÖ®ºóÏÔʾÏÂÃæµÄ´íÎó£º Error using trustnleqn (line 28) Objective function is returning undefined values at initial point. FSOLVE cannot continue. Error in fsolve (line 399) [x,FVAL,JACOB,EXITFLAG,OUTPUT,msgData]=... Error in func_main (line 27) F=fsolve(@(y2)func2(a,q2,x1,x2,y2),x0); |
4Â¥2012-05-28 18:43:31
СÁõ8885
ľ³æ (ÕýʽдÊÖ)
- Ó¦Öú: 3 (Ó×¶ùÔ°)
- ½ð±Ò: 10042.8
- É¢½ð: 200
- Ìû×Ó: 468
- ÔÚÏß: 160.5Сʱ
- ³æºÅ: 735175
- ×¢²á: 2009-03-30
- ÐÔ±ð: MM
- רҵ: ²ÄÁÏÎïÀí»¯Ñ§
|
Õâ¾äÓï¾ä¾ÍÊǻһ¸öË«yÖáµÄµÄͼ£¬ºÃÏñûÓÐ´í°¡¡£ function f4=func_main(name,a,q1,q2,X,flag) % flag=0 represent complex permittivity % flag=1 represent complex permeability data=load(name); xf=data(:,1); x2=X*ones(length(data),1); if (flag==0) y1=data(:,2)+data(:,3)*i; elseif (flag==1) y1=data(:,4)+data(:,5)*i; else disp('0 is complex perimittivity, 1 is complex permeability'); return; end x11=func1(a,q1,y1,x2); %µ÷Óú¯Êýfunc1£¬Çó½âx1 %Hanai formula µÄÇó½â½á¹û %Çó½âµÄÊÇx1£¬Îª¾ø¶Ô½éµç³£Êý x1=x11; a=real(x1); b=imag(x1); c=[a,b]; plotyy(xf,a,xf,b); x0=zeros(length(data),1); F=fsolve(@(y2)func2(a,q2,x1,x2,y2),x0); %µ÷ÓÃǰÊöfunc1ËùÇóµÃµÄx1£¬´úÈ뺯Êýfunc2 %Çó½â³¬Ô½·½³ÌµÄy2Öµ¡£ f4=F; r = [real(f4), imag(f4)]; save('data1.txt', 'c', '-ascii'); save('data2.txt', 'r', '-ascii'); function x1=func1(a,q1,y1,x2) x1=(q1.*x2-x2+y1.*(x2./y1).^a)./(q1+(x2./y1).^a-1); %Çó½âx1 function F=func2(a,q2,x1,x2,y2) F=(((1-q2)*(x2-x1))./(x2.^a)).*(y2.^a)-y2+x1; %Çó½ây2 ÎÒ°Ñ plotyy(xf,a,xf,b); ÕâÒ»Óï¾äÌáµ½Ç°ÃæÈ¥ÊÇ¿ÉÒÔÖ´Ðе쬻¹»³öÁËͼ¡£ÔËÐÐÖ®ºóÌáʾÎÒÎÊÌâ³öÔÚ fsolve ÄǸöÓï¾ä¡£ |
5Â¥2012-05-28 18:46:10
csgt0
ÈÙÓþ°æÖ÷ (ÖøÃûдÊÖ)
²ÊÉ«¹Òͼ
-

ר¼Ò¾Ñé: +2 - Ó¦Öú: 367 (˶ʿ)
- ¹ó±ö: 0.457
- ½ð±Ò: 7807
- É¢½ð: 8897
- ºì»¨: 20
- ɳ·¢: 7
- Ìû×Ó: 2809
- ÔÚÏß: 515.8Сʱ
- ³æºÅ: 449046
- ×¢²á: 2007-11-02
- רҵ: ½ðÊô·Ç¾§Ì¬¡¢×¼¾§ºÍÄÉÃ×¾§
- ¹ÜϽ: ¼ÆËãÄ£Äâ
¡¾´ð°¸¡¿Ó¦Öú»ØÌû
¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
СÁõ8885: ½ð±Ò+29, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸, ¹ûÈ»ÊÇÅ£ÈË£¬´ó²Å¡£ 2012-05-29 10:09:20
СÁõ8885: ½ð±Ò+29, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸, ¹ûÈ»ÊÇÅ£ÈË£¬´ó²Å¡£ 2012-05-29 10:09:20
|
x0=zeros(length(data),1); ¸Ä³Éx0=ones(length(data),1); Èç¹ûx0=0£¬y2=0,F2¾ÍÎÞÇî´óÁËËùÒÔ±¨´í |
» ±¾ÌûÒÑ»ñµÃµÄºì»¨£¨×îÐÂ10¶ä£©

6Â¥2012-05-29 09:21:27
СÁõ8885
ľ³æ (ÕýʽдÊÖ)
- Ó¦Öú: 3 (Ó×¶ùÔ°)
- ½ð±Ò: 10042.8
- É¢½ð: 200
- Ìû×Ó: 468
- ÔÚÏß: 160.5Сʱ
- ³æºÅ: 735175
- ×¢²á: 2009-03-30
- ÐÔ±ð: MM
- רҵ: ²ÄÁÏÎïÀí»¯Ñ§
7Â¥2012-05-29 10:05:46













»Ø¸´´ËÂ¥
СÁõ8885