±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1586  |  »Ø¸´: 8

quan_zhang

Òø³æ (СÓÐÃûÆø)

[½»Á÷] ¡¾ÇóÖú¡¿matlabÓÃlsqnonlin·ÇÏßÐÔÄâºÏµÄÎÊÌ⡾Òѽâ¾ö¡¿ ÒÑÓÐ1È˲ÎÓë

´ó¼ÒºÃ°¡£¬ÎÒÓиöÒ»¸öÀûÓÃmatlabµÄ×îС¶þ³Ë·¨À´Çó½â²ÎÊýµÄÎÊÌ⣬À§ÈÅÁËÎҺþÃÁË¡£
¾ßÌåÇé¿öÈçÏ£ºÎÒÓÐÒ»¸öº¯Êýf£¨x£©=346/15975*sin(4*pi*x/3.195)*pi-3196/3195*sin(2*pi*x/3.195)*piÇÒx=[0,1],ÁíÍâÒ»¸ö´ø²ÎÊýµÄº¯ÊýΪ


            
ÇÒx=[0,1],  a1+a2+a3=1£¬ÇÒa1,a2,a3=(0,1),ÎÒÏëÓÃlsqnonlinÀ´ÄâºÏ¡£
ÎÒµÄ×ö·¨ÈçÏÂ:
1.ÏȽ¨Á¢º¯ÊýÎļþ£¬²¢±£´æÎªmyfunlin.m
function y=myfunlin(x,xdata)
xdata=[0    0.3195  0.639   0.9585  1.278   1.43775 1.5975  1.75725 1.917   2.2365  2.556   2.8755  3.195];
ydata=[0    -1.782446806    -2.948772499    -3.028762147    -1.911872775    -1.011104195    0   1.011104195 1.911872775 3.028762147 2.948772499 1.782446806 1E-15];
y=4.3688*((x(1)*(xdata-x(2)))./((xdata-x(2)).^2+x(3).^2)+(x(4))*(xdata-x(5))./((xdata-x(5)).^2+x(6).^2)+(1-x(1)-x(4))*(xdata-x(7))./((xdata-x(7)).^2+x(8).^2));

2.È»ºóµ÷ÓÃÓÅ»¯³ÌÐò£º
x0=[0.12  3.32  0.16  0.12  -0.13 0.17  0.16  1.2];
x=lsqnonlin('myfunlin',x0)
µÃµ½½á¹û
  -1.4874    2.8738    0.9386   -1.7023    0.3343    0.9831    1.5173    2.0017´ËʱaiµÄ·¶Î§²»ÄÜ¿ØÖÆ£¬ÄÇôӦ¸ÃÈçºÎ¿ØÖÆaiµÄ·¶Î§ÄØ£¿

[ Last edited by nono2009 on 2010-10-31 at 07:24 ]
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

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

change0618

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

·½ÕÉ´óʦ

quan_zhang(½ð±Ò+1): 2010-10-26 15:38:46
help lsqnonlin
2Â¥2010-10-26 13:52:48
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

quan_zhang

Òø³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
Originally posted by change0618 at 2010-10-26 13:52:48:
help lsqnonlin

ʲôÒâ˼°¡£¿ÄܾßÌåµãÂð£¿
3Â¥2010-10-26 15:39:08
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

change0618

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

·½ÕÉ´óʦ

quan_zhang(½ð±Ò+1): 2010-10-27 09:15:46
quan_zhang(½ð±Ò+1): 2010-10-27 09:16:08
X = LSQNONLIN(FUN,X0,LB,UB) defines a set of lower and upper bounds on
    the design variables, X, so that the solution is in the range LB <= X
    <= UB. Use empty matrices for LB and UB if no bounds exist. Set LB(i)
    = -Inf if X(i) is unbounded below; set UB(i) = Inf if X(i) is
    unbounded above.
4Â¥2010-10-26 17:40:52
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lidaxue

ľ³æ (ÕýʽдÊÖ)

Ö®ºõÕßÒ²

quan_zhang(½ð±Ò+2): 2010-10-27 09:16:03
quan_zhang(½ð±Ò+2): 2010-10-27 21:15:11
function seqcurvefit11
clear all;
load seqdata;
beta0=[0.005 0.001];

lb=[0 0];ub=[inf inf];
[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...
    lsqnonlin(@seqfun,beta0,lb,ub,[],t,c);


function y = seqfun(beta,t,c)
ca=exp(-beta(1)*t);
cb=beta(1)/(beta(2)-beta(1))*(exp(-beta(1)*t)-exp(-beta(2)*t));
y=[ca-c(:,1) cb-c(:,2)];
Â¥Ö÷ ÒÀºù«»­Æ°°É~
Comeon£¡
5Â¥2010-10-26 18:48:56
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

quan_zhang

Òø³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
Originally posted by lidaxue at 2010-10-26 18:48:56:
function seqcurvefit11
clear all;
load seqdata;
beta0=[0.005 0.001];

lb=[0 0];ub=[inf inf];
[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...
    lsqnonlin(@seqfun,beta0,lb,ub, ...

residual,exitflag,output,lambda,jacobian
¸ÃÔõô¶¨Òå°¡
6Â¥2010-10-28 09:25:45
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

quan_zhang

Òø³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
Originally posted by lidaxue at 2010-10-26 18:48:56:
function seqcurvefit11
clear all;
load seqdata;
beta0=[0.005 0.001];

lb=[0 0];ub=[inf inf];
[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...
    lsqnonlin(@seqfun,beta0,lb,ub, ...

ÊDz»ÊÇÒª¶Ôt£¬c¸³Öµ°¡£¬ÎÒÊäÈë³öÏÖ Undefined function or variable 't'.
7Â¥2010-10-28 09:43:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lidaxue

ľ³æ (ÕýʽдÊÖ)

Ö®ºõÕßÒ²

¡ï
nono2009(½ð±Ò+1):¹ÄÀøÓ¦Öú 2010-10-28 14:02:08
quan_zhang(½ð±Ò+3): 2010-10-30 15:01:21
ÒýÓûØÌû:
Originally posted by quan_zhang at 2010-10-28 09:43:51:

ÊDz»ÊÇÒª¶Ôt£¬c¸³Öµ°¡£¬ÎÒÊäÈë³öÏÖ Undefined function or variable 't'.

ÔΠ Â¥Ö÷ûÓÐÀí½âÎÒµÄÒâ˼°¡£¬[]ºóÃæµÄt£¬cÊÇ´«µÝµÄ²ÎÊý£¬ÕâÀïÖ»¸øÄãÒ»¸öÀý×Ó£¬ÈÃÄã¸ù¾Ý×Ô¼ºµÄÐèÒªÐÞ¸ÄһϵÄ
Comeon£¡
8Â¥2010-10-28 13:12:05
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

quan_zhang

Òø³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
Originally posted by lidaxue at 2010-10-28 13:12:05:

ÔΠ Â¥Ö÷ûÓÐÀí½âÎÒµÄÒâ˼°¡£¬[]ºóÃæµÄt£¬cÊÇ´«µÝµÄ²ÎÊý£¬ÕâÀïÖ»¸øÄãÒ»¸öÀý×Ó£¬ÈÃÄã¸ù¾Ý×Ô¼ºµÄÐèÒªÐÞ¸ÄһϵÄ

ллÄã°¡£¬ÎÒÒѾ­¸ã¶¨ÁË¡£
9Â¥2010-10-30 15:01:54
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ quan_zhang µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Ò»Ö¾Ô¸Çà¿Æ085500£¬³õÊÔ295·Ö£¬¹«¹²¿Î213·Ö +3 Óöµ½µÄÈËÔ¸Íû¶¼Ä 2026-04-05 3/150 2026-04-05 18:45 by À¶ÔÆË¼Óê
[Óлú½»Á÷] ¼×´¼/¶þÂÈ 1:15¹ýÖù×Ó 5+3 a°¥yßÏwι 2026-03-31 3/150 2026-04-05 10:42 by 88817753
[¿¼ÑÐ] 323·Ö£¨¼ÆËã»úÊÓ¾õºÍ´óÄ£ÐÍÏîÄ¿£©ÄÜÖ±½ÓÉÏÊÖ +3 chaoxiicy 2026-04-01 3/150 2026-04-05 00:50 by chongya
[¿¼ÑÐ] Çóµ÷¼Á +7 xzghyuj 2026-04-04 7/350 2026-04-04 22:25 by oooqiao
[¿¼ÑÐ] Ò»Ö¾Ô¸ÉϺ£´óѧÉúÎïѧ346 +3 ÉϺ£´óѧ346µ÷¼Á 2026-04-03 3/150 2026-04-04 20:20 by dongzh2009
[¿¼ÑÐ] Ò»Ö¾Ô¸±±½»´ó²ÄÁϹ¤³Ì×Ü·Ö358 +6 cs0106 2026-04-03 6/300 2026-04-04 11:20 by w_xuqing
[¿¼ÑÐ] һ־Ը˫·Ç085502£¬267·Ö£¬¹ýËļ¶Çóµ÷¼Á +3 ÔÙæҲҪ³Ô·¹°¡ 2026-04-03 3/150 2026-04-04 05:03 by gswylq
[¿¼ÑÐ] 334Çóµ÷¼Á +9 Trying] 2026-03-31 9/450 2026-04-03 15:18 by ×ÁçíØ¼
[¿¼ÑÐ] 085600 295·ÖÇóµ÷¼Á +19 W55j 2026-03-30 23/1150 2026-04-03 09:53 by ǧǧÔËÆø
[¿¼ÑÐ] Çóµ÷¼Á£¡ÉúÎïÓëҽҩר˶ +4 Äæ×ªÂ½ÏÈÉú 2026-04-01 5/250 2026-04-03 08:33 by Jaylen.
[¿¼ÑÐ] Ò»Ö¾Ô¸aÇø211£¬085601-307·ÖÇóµ÷¼Á +13 µ³¼ÎºÀ 2026-03-31 26/1300 2026-04-03 08:33 by 495374996
[¿¼ÑÐ] ѧ˶»¯Ñ§¹¤³ÌÓë¼¼Êõ£¬Ò»Ö¾Ô¸Öйúº£Ñó´óѧ320+Çóµ÷¼Á +8 ÅûÐÇºÓ 2026-04-02 8/400 2026-04-02 14:12 by oooqiao
[¿¼ÑÐ] ²ÄÁÏ»¯¹¤340Çóµ÷¼Á +5 jhx777 2026-03-30 5/250 2026-04-02 12:45 by smileboy2006
[¿¼ÑÐ] ²ÄÁÏר˶322·Ö +11 ¹þ¹þ¹þºðºðºð¹þ 2026-04-01 11/550 2026-04-02 10:52 by lnilvy
[¿¼ÑÐ] °²È«¹¤³Ì 285 Çóµ÷¼Á +3 Xinyu56 2026-04-01 4/200 2026-04-01 21:50 by ¾²¾²¾²¾²¾²¾²¾²¾
[¿¼ÑÐ] ²ÄÁÏרҵµ÷¼Á +5 À²À²À²¿Þ 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[¿¼ÑÐ] 318Çóµ÷¼Á +8 ÆßÒä77 2026-04-01 8/400 2026-04-01 10:37 by Jaylen.
[¿¼ÑÐ] 339Çóµ÷¼Á +5 zjjkt 2026-03-31 5/250 2026-04-01 09:18 by JourneyLucky
[¿¼ÑÐ] ²ÄÁÏר˶ 085600Çóµ÷¼Á +7 BBQ233 2026-03-30 7/350 2026-03-30 17:44 by oooqiao
[¿¼ÑÐ] 293Çóµ÷¼Á +3 ĩδmm 2026-03-30 5/250 2026-03-30 17:23 by Íõ±£½Ü33
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û