±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 2197  |  »Ø¸´: 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µÄ»ØÌû

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µÄ»ØÌû

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µÄ»ØÌû

xxz903

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

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

1stopt¿ÉÒÔ×Ô¶¯Ñ°ÕÒ¹«Ê½µÄ~~¶øÇÒÄâºÏ¶È»¹Âù¸ßµÄ~
6Â¥2013-05-14 21:41:15
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

pianfeng1986

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

ÒýÓûØÌû:
6Â¥: Originally posted by xxz903 at 2013-05-14 21:41:15
1stopt¿ÉÒÔ×Ô¶¯Ñ°ÕÒ¹«Ê½µÄ~~¶øÇÒÄâºÏ¶È»¹Âù¸ßµÄ~

ÎÒÕâ¸öÆäʵ·ÇÄâºÏÎÊÌâÁË£¬ÓÅ»¯Çó×îСֵµÄÎÊÌâ~~ ²»ÖªµÀÄúÏë·¨ÊÇÈçºÎµÄ£¿
7Â¥2013-05-16 15:06:56
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

xxz903

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

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

¡ï
csgt0: ½ð±Ò+1, лл 2013-05-16 17:14:33
×îСֵÎÊÌâ¿ÉÒÔʹÓÃmatlabµÄÓÅ»¯¹¤¾ßÏ䣬Ҳ¿ÉÒÔʹÓÃ1stoptµÄÓÅ»¯¹¦ÄÜ£¨min  »òÕßmax£© ¶¼¿ÉÒԵĠ  ÎÊÌâ²»´ó
8Â¥2013-05-16 16:10:41
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ pianfeng1986 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ÉúÎïѧ308·ÖÇóµ÷¼Á£¨Ò»Ö¾Ô¸»ª¶«Ê¦´ó£© +6 ÏàÐűػá¹ââÍòÕ 2026-04-06 7/350 2026-04-06 11:51 by lijunpoly
[¿¼ÑÐ] 085410È˹¤ÖÇÄÜ ³õÊÔ316·Ö Çóµ÷¼Á +7 ²ÐÐÇ·÷Êï 2026-03-31 7/350 2026-04-06 10:09 by À¶ÔÆË¼Óê
[¿¼ÑÐ] Çóµ÷¼Á +7 ÕÅ.1 2026-04-05 7/350 2026-04-05 20:40 by à£à£à£0119
[¿¼ÑÐ] 282Çóµ÷¼Á +7 aaa³µÁ¾ 2026-04-02 11/550 2026-04-05 17:24 by yulian1987
[¿¼ÑÐ] µ÷¼Á +3 ºÃºÃ¶ÁÊé¡£ 2026-04-02 3/150 2026-04-05 13:02 by arrow8852
[¿¼ÑÐ] Çóµ÷¼Á +3 С›g 2026-04-03 3/150 2026-04-05 09:10 by sihailian3
[¿¼ÑÐ] 290Çóµ÷¼Á +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
[¿¼ÑÐ] 344²ÄÁÏÓ뻯¹¤µ÷¼Á +9 µ÷¼ÁÉϰ¶«^ 2026-04-03 9/450 2026-04-04 23:10 by happyddm
[¿¼ÑÐ] ÇóÉúÎïѧѧ˶µ÷¼Á¡ª¡ª364·Ö +7 ÔÆ¶äåÞÍäÖ¸ÄÏ 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
[¿¼ÑÐ] 359Çóµ÷¼Á +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[¿¼ÑÐ] 321Çóµ÷¼Á +13 ÈÏÕæÇóÉÏѧ 2026-04-02 13/650 2026-04-04 18:23 by macy2011
[¿¼ÑÐ] 311Çóµ÷¼Á +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[¿¼ÑÐ] 295Çóµ÷¼Á +3 ÉÐټѽ 2026-04-03 4/200 2026-04-03 21:23 by zhq0425
[˶²©¼ÒÔ°] ÇóÀÏʦÊÕÁô +9 lllq123 2026-04-03 9/450 2026-04-03 13:48 by ºôÎü¶¼ÊǼõ·Ê
[¿¼ÑÐ] ½¨»·£¬ÄÜÔ´£¬ÍÁľÀÏʦ·¹ý¿´Ò»¿´£¡£¡£¡ +5 ºÙºÙuu 2026-04-01 5/250 2026-04-03 11:47 by znian
[¿¼ÑÐ] 325·Ö»¯Ñ§µ÷¼Á +5 15771691647 2026-04-02 5/250 2026-04-03 09:58 by ChemPharm
[¿¼ÑÐ] 285Çóµ÷¼Á +11 AZMK 2026-04-01 11/550 2026-04-01 22:40 by peike
[¿¼ÑÐ] 085600£¬321·ÖÇóµ÷¼Á +13 ´ó²öС×Ó 2026-03-31 13/650 2026-04-01 12:35 by chemdavid
[¿¼ÑÐ] 323·Ö ʳƷÓëÓªÑøµ÷¼Á +3 ºÙooo 2026-03-31 3/150 2026-03-31 09:38 by longlotian
[¿¼ÑÐ] ±¾¿Æ211×Ü·Ö289£¬08¹¤Ñ§ÕæÐÄÇóµ÷¼Á +3 utopiaE 2026-03-30 3/150 2026-03-30 23:42 by ms629
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û