| 查看: 219 | 回复: 0 | ||||
| 【悬赏金币】回答本帖问题，作者linkaiglk将赠送您 20 个金币 | ||||
[求助]
看看这段程序 我不知道从哪里开始粘贴
|
||||
|
这是给的程序 我是个matlab小白 不知道那一段是要粘贴到matlab里的？ function xmc02 clear all clc % deta global A r=[0.0051,0.0065,0.0073,0.008,0.0086,0.0043,0.0059,0.0067,0.0073,0.0079,0.0065,0.0039,0.0024,0.0019,0.0013]'; X1=[0.1 0.15 0.2 0.25 0.3 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2]'; X2=[0.8 0.8 0.8 0.8 0.8 0.2 0.4 0.6 0.8 1 0.8 0.8 0.8 0.8 0.8]'; X3=[0 0 0 0 0 0 0 0 0 0 0.02 0.04 0.06 0.08 0.1]'; A=[X1,X2,X3]; % objfun r=a*b*c*X1*X2/(1+aX1+bX2+dX3)^2 % 非线性拟合 beta0=[1 2 1 1]; % 依次为 a b c d [beta,resnorm,resid,exitflag] = ... lsqnonlin(@OptObjFunc,beta0,[],[],[],r,A) % ------------------------------------------------------------------ function f = OptObjFunc(beta,r,A) f = r-beta(1)*beta(2)*beta(3).*A(:,1).*A(:,2)./(1+beta(1).*A(:,1)+beta(2).*A(:,2)+beta(4).*A(:,3)).^2; 运行结果： Optimization terminated: relative function value changing by less than OPTIONS.TolFun. beta = 5.006082883296007 1.799732413072861 0.061596203399673 35.660702878551518 resnorm = 2.798689414484294e-006 resid = -0.000034980338998 -0.000041384319011 -0.000199113436192 -0.000145723489928 0.000027180441694 0.000318305366339 -0.000095947962104 -0.000315230960303 -0.000199113436192 0.000217442282819 0.001354807417095 0.000152165876446 -0.000451091514209 -0.000341537035504 -0.000508450568973 exitflag = 3  QQ截图20140701204831.png |
回复此楼» 收录本帖的淘贴专辑推荐
 文献 |
» 猜你喜欢
杭电、天津科技、青农和宁波工程学院如何选？
已经有8人回复
-
关于2023的收入
已经有32人回复
-
关于博后基金的bug问题
已经有6人回复
-
物理化学一对一辅导
已经有5人回复
-
车间生产，真空度很高，温度很高，但减压蒸馏速度很慢。
已经有32人回复
-
为什么我的博后基金还在流动站审核中？不会是学院给我卡了吧？
已经有22人回复
-
BO4的YQ答辩通知发布了吗？
已经有4人回复
-
Nature 11日发文，中国著名学者们称造假迫不得已
已经有3人回复
-
博后基金，以往的结果点不开，怎么回事呢？最后一次机会了，两次都没中前面。
已经有10人回复
-
请问事业编制和年薪制冲突吗？
已经有4人回复
点击这里搜索更多相关资源
5