²é¿´: 2177  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 268Çóµ÷¼Á +6 Ò»¶¨ÓÐѧÉÏ- 2026-03-14 7/350 2026-03-17 10:14 by Sammy2
[¿¼ÑÐ] [¸´ÊÔµ÷¼Á]Î÷ÄϿƼ¼´óѧ¹ú·À/²ÄÁϵ¼Ê¦ÍƼö +3 ¼â½ÇСºÉ 2026-03-16 6/300 2026-03-16 23:21 by ¼â½ÇСºÉ
[¿¼ÑÐ] ҩѧ383 Çóµ÷¼Á +3 ҩѧchy 2026-03-15 4/200 2026-03-16 20:51 by Ôª×Ó^0^
[¿¼ÑÐ] 085601Çóµ÷¼Á +3 Du.11 2026-03-16 3/150 2026-03-16 20:42 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 328Çóµ÷¼Á£¬Ó¢ÓïÁù¼¶551£¬ÓпÆÑо­Àú +3 ÉúÎ﹤³Ìµ÷¼Á 2026-03-16 4/200 2026-03-16 20:13 by Wangjingyue
[¿¼ÑÐ] 311Çóµ÷¼Á +5 26ÑÐ0 2026-03-15 5/250 2026-03-16 16:21 by a²»Ò×
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤Çóµ÷¼Á +3 Ϊѧ666 2026-03-16 3/150 2026-03-16 15:09 by ¼ÓºÅ+
[¿¼ÑÐ] 0703 ÎïÀí»¯Ñ§µ÷¼Á +3 ÎÒ¿ÉÒÔÉϰ¶µÄ¶Ô 2026-03-13 5/250 2026-03-16 10:50 by ÎÒ¿ÉÒÔÉϰ¶µÄ¶ÔÂ
[¿¼ÑÐ] 326Çóµ÷¼Á +4 Éϰ¶µÄСÆÏ 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[¿¼ÑÐ] ¡¾0703»¯Ñ§µ÷¼Á¡¿-Ò»Ö¾Ô¸»ªÖÐʦ·¶´óѧ-Áù¼¶475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by ŶŶ123
[¿¼ÑÐ] 321Çóµ÷¼Á +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[¿¼ÑÐ] 304Çóµ÷¼Á +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] ËÄ´¨´óѧ085601²ÄÁϹ¤³Ìר˶ ³õÊÔ294Çóµ÷¼Á +4 ×£ÎÒÃǺÃÔÚ¶¬Ìì 2026-03-11 4/200 2026-03-13 21:39 by peike
[¿¼ÑÐ] ¡¾0856¡¿»¯Ñ§¹¤³Ì£¨085602£©313 ·Ö£¬±¾¿ÆÑ§¿ÆÆÀ¹ÀAÀàԺУ»¯Ñ§¹¤³ÌÓ빤ÒÕ£¬³ÏÇóµ÷¼Á +7 СÁõ¿ì¿ìÉϰ¶ 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[¿¼ÑÐ] ÍÁľµÚÒ»Ö¾Ô¸276Çóµ÷¼Á£¬¿ÆÑкͼ¼ÄÜÊ®·Ö·á¸»£¬ÇóÐÂÐË·½ÏòµÄµ¼Ê¦ÊÕÁô +3 ÍÁľСÌì²Å 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[¿¼ÑÐ] 26¿¼ÑÐÇóµ÷¼Á +5 ؼºêSir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[¿¼ÑÐ] 290Çóµ÷¼Á +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[¿¼ÑÐ] ¹¤¿Æ0856ר˶»¯Ñ§¹¤³Ì269Äܵ÷¼ÁÂð +10 ÎÒÏë¶ÁÑÐ11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[¿¼ÑÐ] 333Çóµ÷¼Á +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[¿¼ÑÐ] 298Çóµ÷¼Á +3 Vvѽ£¡ 2026-03-10 3/150 2026-03-10 22:40 by ½£Ê«¶Å¿µ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û