Znn3bq.jpeg
²é¿´: 2204  |  »Ø¸´: 7
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

pianfeng1986

Ìú³æ (³õÈëÎÄ̳)

[ÇóÖú] Çó½ÌmatlabÏßÐÔÓÅ»¯Çó×îСֵµÄÎÊÌâ

ÕâÁ½ÌìÓöµ½ÕâÑùÒ»¸öÎÊÌ⣬ ĿǰÓöµ½Ò»¸öº¯Êý F(a,b,c,x)£¬ Èç¹û¶¨Òåa,b,c£¬¿ÉÒÔ»­³öFÔÚxÒ»¸öÇø¼äÉϵÄÇúÏߣ¬µ«ÎÞ·¨ÇóµÃFµÄ¾ßÌå½âÎöʽ¡£  ÎÒµÄÄ¿µÄÊǰÑFÇúÏßÓëÊÖÍ·ÉϵÄʵÑéÇúÏß´ïµ½×îºÃÎǺϣ¬Í¬Ê±Çó³öÏàÓ¦µÄ²ÎÊýa,b,c¡£

      Ä¿Ç°ÎÒµÄ˼·ÊÇÕâÑùµÄ£ºÊÔÑéÇúÏßÓÐ10000×éÊý¾Ý £¨x,y)£¬ÎÒÓÃ×îС¶þ³ËµÄ˼Ï룬Çó³öÿһ¸öxÉ϶ÔÓ¦µÄFÓëyÖ®²îƽ·½ºÍ£¬¼´sum((y-F)^2),  È»ºó¸ø³öa, b , cµÄ·¶Î§£¬ÓÃfminconÀ´Çó½â£¬µ«ÊÇдÁ˳ÌÐò£¬·¢ÏÖÎÊÌ⣬³ÌÐòÎÞ·¨Çó³öÏàÓ¦½â¡£´íÎóΪ  Warning: Large-scale (trust region) method does not currently solve this type of problem,
switching to medium-scale (line search).
> In fmincon at 274
  In main at 19
Optimization terminated: magnitude of directional derivative in search
direction less than 2*options.TolFun and maximum constraint violation
  is less than options.TolCon.
No active inequalities.


ÇëÎÊ´óÉñÎÒ¸ÃÔõô°ì£¬ÊǺ¯ÊýÓôíÁËÂ𣿠»¹ÊÇÓñðµÄ·½·¨£¿

Ö÷³ÌÐòÈçÏ£º
clear all;clc;
global nr nt apt pidt vb freq nfrq

ui=sqrt(-1);
nr=50;
nt=51;
apt=100;
pidt=10d0;
vb=4212;

%% A£¬b¶¨ÒåÔ¼Êø     45 < x1<55;  0.22 < x2<0.28;  0.0005 < x3<0.001;
A=[1,0,0; 0,1,0;0,0,1;-1,0,0;0,-1,0;0,0,-1];
b=[55,0.28,0.001,-45,-0.22,-0.0005]';
x0= [49.13; 0.2592; 0.0007];    % Starting guess at the solution
[x,fval] = fmincon(@myfun,x0,A,b)
%%_______________________________________
function f = myfun(x);
global nr nt apt pidt vb freq nfrq
adf=load('y_para.dat');

fs=adf(:,1);
con=adf(:,2);
freq=fs;
nfrq=length(freq);

epson= 0.3302  ;
ita=0.6052  ;
c=1.0888  ;
ym=admres(x(1),x(2),epson,ita,c,x(3));
f=sum((real(ym)-con).^2);
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

pianfeng1986

Ìú³æ (³õÈëÎÄ̳)

ÒýÓûØÌû:
2Â¥: Originally posted by csgt0 at 2013-05-14 15:54:36
ÓÃlsqcurvefitÊÔÊÔ

лл ÎÒÑо¿ÏÂ
5Â¥2013-05-14 21:36:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 8 ¸ö»Ø´ð

csgt0

ÈÙÓþ°æÖ÷ (ÖøÃûдÊÖ)

²ÊÉ«¹Òͼ

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
pianfeng1986: ½ð±Ò+5 2013-05-16 15:07:20
ÓÃlsqcurvefitÊÔÊÔ
showmethemoney
2Â¥2013-05-14 15:54:36
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

xxz903

½ð³æ (СÓÐÃûÆø)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
csgt0: ½ð±Ò+1, лл 2013-05-14 16:19:41
pianfeng1986: ½ð±Ò+5, ¡ïÓаïÖú 2013-05-16 15:07:06
ÏßÐÔÎÊÌâÖ±½ÓÓÃorigin¾ÍÐа¡~matlabÓÅ»¯×îСÎó²îÒ²ÐÐ~~~²»¹ýÓÅ»¯ÄâºÏ×îºÃÓõϹÊÇ1stopt~ÓÐ1.5ÆÆ½âµÄ  ¾ø¶ÔÄܽâ¾öÄãµÄÎÊÌâ
3Â¥2013-05-14 16:01:07
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

pianfeng1986

Ìú³æ (³õÈëÎÄ̳)

ÒýÓûØÌû:
3Â¥: Originally posted by xxz903 at 2013-05-14 16:01:07
ÏßÐÔÎÊÌâÖ±½ÓÓÃorigin¾ÍÐа¡~matlabÓÅ»¯×îСÎó²îÒ²ÐÐ~~~²»¹ýÓÅ»¯ÄâºÏ×îºÃÓõϹÊÇ1stopt~ÓÐ1.5ÆÆ½âµÄ  ¾ø¶ÔÄܽâ¾öÄãµÄÎÊÌâ

ÎÒËùÇóµÄº¯ÊýľÓоßÌåµÄ½âÎöʽ£¬1stopt²»ºÃÓð¡
4Â¥2013-05-14 21:35:05
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 0702ÎïÀíѧѧ˶299Çóµ÷¼Á +5 ÆîÆâÁ¬ 2026-04-06 5/250 2026-04-09 09:27 by wj165256
[¿¼ÑÐ] µç×ÓÐÅÏ¢270Çóµ÷¼Á +8 terminal469 2026-04-07 8/400 2026-04-09 08:56 by ¿ÆÑÐëë³æ»¯µû
[¿¼ÑÐ] ÉúÎïѧ328·ÖÇóµ÷¼Á +9 ÉÁµçkkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[¿¼ÑÐ] 347²ÄÁÏר˶Çóµ÷¼Á +18 zj8215216 2026-04-06 18/900 2026-04-08 16:27 by luoyongfeng
[¿¼ÑÐ] Ò»Ö¾Ô¸¼ª´ó»¯Ñ§327Çóµ÷¼Á +12 ÍõÍõ°×ʯ 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[¿¼ÑÐ] 283Çóµ÷¼Á +19 A child 2026-04-04 19/950 2026-04-08 14:26 by xingguangj
[¿¼²©] É격 +8 IQwQl 2026-04-04 8/400 2026-04-08 09:43 by 0608104024
[¿¼ÑÐ] 305·ÖÇóµ÷¼Á +3 ¹þ_¹þ_¹þ_¹þ_¹þ 2026-04-04 5/250 2026-04-07 14:49 by ¹þ_¹þ_¹þ_¹þ_¹þ
[¿¼ÑÐ] µ÷¼Á +4 mcbbc 2026-04-06 5/250 2026-04-07 12:33 by upczlm1989
[¿¼ÑÐ] ÉúÎïÓëÒ½Ò©086000µ÷¼ÁÒ»Ö¾Ô¸Î÷±±Å©ÁÖ320·Ö +3 ÃÀÃÀŮʿ 2026-04-03 3/150 2026-04-05 21:55 by ѧԱ8dgXkO
[¿¼ÑÐ] 284Çóµ÷¼Á +7 Ðìͬѧ_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[¿¼ÑÐ] 085600µ÷¼Á +9 ¶«ÕÕÕÕÕÕ 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +7 dxyµ÷¼Á 2026-04-04 7/350 2026-04-05 09:15 by İÇï26
[¿¼ÑÐ] 320Çóµ÷¼Á +3 Ò»ÑùÔ² 2026-04-04 3/150 2026-04-04 22:29 by à£à£à£0119
[¿¼ÑÐ] 341Çóµ÷¼Á +3 Âå¶àÂÞ 2026-04-02 4/200 2026-04-04 21:36 by ÖÇÄÜÖÇ»Û
[¿¼ÑÐ] 359Çóµ÷¼Á +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[¿¼ÑÐ] ±¾¿Æ211£¬×¨Òµ085404£¬293·ÖÇëÇóµ÷¼Á +5 Á«²Ë¾ÍÊÇź°É 2026-04-04 5/250 2026-04-04 14:08 by ÕâÊÇÒ»¸öÎÞÁĵÄê
[¿¼ÑÐ] Çóµ÷¼Á +4 ѹÁ¦??´ó 2026-04-03 4/200 2026-04-03 21:36 by à£à£à£0119
[¿¼ÑÐ] 320Çóµ÷¼Á +3 ũҵ¹¤³ÌÓëÐÅÏ¢¼ 2026-04-03 3/150 2026-04-03 11:40 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] 22408 266Çóµ÷¼Á +3 masss11222 2026-04-02 3/150 2026-04-02 18:11 by ±ÊÂä½õÖÝ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û