²é¿´: 2185  |  »Ø¸´: 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µÄ»ØÌû

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µÄ»ØÌû
²é¿´È«²¿ 8 ¸ö»Ø´ð

csgt0

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

²ÊÉ«¹Òͼ

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

¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
pianfeng1986: ½ð±Ò+5 2013-05-16 15:07:20
ÓÃlsqcurvefitÊÔÊÔ
showmethemoney
2Â¥2013-05-14 15:54:36
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢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µÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 286Çóµ÷¼Á +6 lemonzzn 2026-03-16 9/450 2026-03-18 15:31 by vgtyfty
[¿¼ÑÐ] 344Çóµ÷¼Á +5 knight344 2026-03-16 6/300 2026-03-18 15:17 by knight344
[¿¼ÑÐ] 266Çóµ÷¼Á +5 ÑôÑôÍÛÈû 2026-03-14 9/450 2026-03-18 15:05 by stone_128
[¿¼ÑÐ] »¯Ñ§¹¤³Ì321·ÖÇóµ÷¼Á +15 ´óÃ×·¹£¡ 2026-03-15 18/900 2026-03-18 14:52 by haxia
[¿¼ÑÐ] Ò»Ö¾Ô¸ÎäÀí²ÄÁÏ305·ÖÇóµ÷¼Á +3 ÏëÉϰ¶µÄÀðÓã 2026-03-18 4/200 2026-03-18 14:33 by laoshidan
[¿¼ÑÐ] 296Çóµ÷¼Á +5 ´ó¿Ú³Ô·¹ ÉíÌ彡 2026-03-13 5/250 2026-03-17 21:05 by ²»»ó¿ÉÀÖ
[¿¼ÑÐ] 085601Çóµ÷¼Á +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[¿¼ÑÐ] ±¾ÈË¿¼085602 »¯Ñ§¹¤³Ì ר˶ +16 ²»ÖªµÀ½Ðʲô£¡ 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤Çóµ÷¼Á +5 Ð÷ÐÒÓë×Ó 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[¿¼ÑÐ] ²ÄÁϹ¤³Ìר˶274Ò»Ö¾Ô¸211Çóµ÷¼Á +6 Ñ¦ÔÆÅô 2026-03-15 6/300 2026-03-17 11:05 by ѧԱh26Tkc
[ÂÛÎÄͶ¸å] ÓÐûÓдóÀз¢Ð¡ÂÛÎÄÄÜ´øÎÒ¸ö¶þ×÷ +3 ÔöÈñ©ÈË 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[¿¼ÑÐ] »¯Ñ§µ÷¼Á0703 +8 °¡ÎÒÎÒµÄ 2026-03-11 8/400 2026-03-16 17:23 by ÎҵĴ¬Îҵĺ£
[¿¼ÑÐ] 304Çóµ÷¼Á +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by ÎҵĴ¬Îҵĺ£
[¿¼ÑÐ] 277²ÄÁÏ¿ÆÑ§Ó빤³Ì080500Çóµ÷¼Á +3 ×ÔÓɼå±ý¹û×Ó 2026-03-16 3/150 2026-03-16 14:10 by ÔËÆøyunqi
[¿¼ÑÐ] 0703 ÎïÀí»¯Ñ§µ÷¼Á +3 ÎÒ¿ÉÒÔÉϰ¶µÄ¶Ô 2026-03-13 5/250 2026-03-16 10:50 by ÎÒ¿ÉÒÔÉϰ¶µÄ¶ÔÂ
[¿¼ÑÐ] ѧ˶285Çóµ÷¼Á +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[¿¼ÑÐ] 330Çóµ÷¼Á +3 ?½´¸øµ÷¼Á¹òÁË 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[¿¼ÑÐ] £¨081700£©»¯Ñ§¹¤³ÌÓë¼¼Êõ-298·ÖÇóµ÷¼Á +12 11À²À²À² 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[¿¼ÑÐ] 295Çóµ÷¼Á +3 Сذ×ÐÖ­ 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[¿¼ÑÐ] 283Çóµ÷¼Á£¬²ÄÁÏ¡¢»¯¹¤½Ô¿É +8 ËÕ´òË®7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û