Znn3bq.jpeg
²é¿´: 2205  |  »Ø¸´: 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µÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 297Çóµ÷¼Á +25 GENJIOW 2026-04-07 28/1400 2026-04-09 18:47 by laoshidan
[¿¼ÑÐ] 305Çóµ÷¼Á +4 77Qi 2026-04-07 4/200 2026-04-09 17:27 by wp06
[¿¼ÑÐ] 314Çóµ÷¼Á +20 wakeluofu 2026-04-09 21/1050 2026-04-09 16:14 by zhuimr
[¿¼ÑÐ] »¯Ñ§¹¤³ÌÓë¼¼Êõרҵһ־Ը¹þ¹¤³Ì 291·ÖBÇø ¹ú¼Ò¼¶´ó´´¸ºÔðÈË ÓÐÒ»×÷ÂÛÎÄ +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by onlyÖÜ
[¿¼ÑÐ] Ò»Ö¾Ô¸0807 ÊýÒ»Ó¢Ò» 313 ÓÐûÓжþÂÖµ÷¼Á +11 emokidd 2026-04-08 12/600 2026-04-09 09:24 by wyf236
[¿¼ÑÐ] 0703µ÷¼Á +17 ʰ¾ÁÒ¼ 2026-04-04 19/950 2026-04-09 09:18 by may_ÐÂÓî
[¿¼ÑÐ] ¿¼ÑÐÇóµ÷¼Á +4 ö©??? 2026-04-08 4/200 2026-04-08 21:44 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] 298Çóµ÷¼Á +6 ¶¤¶£ß˶¬¹Ï 2026-04-07 8/400 2026-04-08 10:51 by ÖзÉÔº¿Õ¹ÜѧԺÑ
[¿¼ÑÐ] һ־Ը̫ԭÀí¹¤´óѧ¼ÆËã»ú¼¼Êõר˶348£¬Çóµ÷¼ÁÖ¸µ¼ +3 nexious 2026-04-05 3/150 2026-04-07 08:19 by jp9609
[¿¼ÑÐ] Çóµ÷¼Á +4 µçÆøÐ¡Éñͯ 2026-04-04 6/300 2026-04-07 00:14 by guanxin1001
[¿¼ÑÐ] 308Çóµ÷¼Á +4 maverick^_^ 2026-04-03 4/200 2026-04-05 19:08 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 295Çóµ÷¼Á +8 FZAC123 2026-04-03 8/400 2026-04-05 17:46 by À¶ÔÆË¼Óê
[¿¼ÑÐ] ¿¼Ñе÷¼Á +5 ËÄ´¨ÍõÌÎ 2026-04-04 5/250 2026-04-04 22:18 by à£à£à£0119
[¿¼ÑÐ] 319Çóµ÷¼Á +4 ÐÇÐDz»Õ£ÑÛà¶ 2026-04-03 4/200 2026-04-04 16:25 by ÖзÉÔº¿Õ¹ÜѧԺÑ
[¿¼ÑÐ] ÔõôɾÌû×Ó°¡ +3 ·ìêØ1000 2026-04-04 3/150 2026-04-04 14:20 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] 26µ÷¼Á 086003 +6 ʧ»îµÄϸ°û 2026-04-04 6/300 2026-04-04 09:50 by zhangdingwa
[¿¼ÑÐ] ²ÄÁÏ¿ÆÑ§Ó빤³Ì¿¼ÑÐ +10 Õü¾ÈÆ¤ÌØÍÐÏÈÉú 2026-04-02 10/500 2026-04-03 23:57 by userper
[¿¼ÑÐ] 320µ÷¼Á +4 ũҵ¹¤³ÌÓëÐÅÏ¢¼ 2026-04-03 4/200 2026-04-03 21:40 by lbsjt
[¿¼ÑÐ] µ÷¼Á +5 asdasdassda 2026-04-03 6/300 2026-04-03 20:27 by °¶ÉϵÄÒ»ÌõÓã
[¿¼ÑÐ] Ò»Ö¾Ô¸±±¾©½»Í¨´óѧ²ÄÁϹ¤³Ì×Ü·Ö358 +4 cs0106 2026-04-03 4/200 2026-04-03 13:41 by °ÙÁéͯ888
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û