24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wtisunday

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 1875.3
散金: 30
红花: 1
帖子: 936
在线: 568.7小时
虫号: 444484
注册: 2007-10-27
性别: GG
专业: 高分子材料的加工与成型

[求助] 非线性方程组matlab拟合

各位大侠，我想搞一下非线性最小二乘拟合，但我的变量x、y分别都是向量，我想用matlab的高斯牛顿法求解，我编的如下：
function ff=fun(vector)
x=vector(1)
y=vector(2)
for i=1:4
for j=1:4
      A(i,j)=0,B(i,j)=0,b(i)=0,c(i)=0
      for k=1:10
         A(i,j)=A(i,j)+(w(k)*x(j))^2*(w(k)*x(i))^2/(G1(k))^2/(1+(w(k)*...
            x(j))^2)/(1+(w(k)*x(i))^2);
         B(i,j)=B(i,j)+2*y(i)*w(k)^2*x(i)*(w(k)*x(j))^2/(G1(k)^2)/(1+...
            (w(k)*x(i))^2)^2/(1+(w(k)*x(j))^2);
         b(i)=b(i)+(w(k)*x(i))^2/G1(k)/(1+(w(k)*x(i))^2);
         c(i)=c(i)+2*y(i)*x(i)*(w(k))^2/G1(k)/(1+(w(k)*x(i))^2);
      end
end
end
b=b';
c=c';
ff(1)=A*x-b;
ff(2)=B*y-c;

在赋值x=[1 100 1000 10000],y=[100000,20000,900,20]后出现如下问题：
fsolve(@fun1,x0,y0)

x =

0.0100

y =

0.2000

A =

   0

B =

   0

b =

   0

c =

   0

??? Undefined function or method 'w' for input arguments of
type 'double'.

Error in ==> fun1 at 8
         A(i,j)=A(i,j)+(w(k)*x(j))^2*(w(k)*x(i))^2/(G1(k))^2/(1+(w(k)*...

Error in ==> fsolve at 254
         fuser = feval(funfcn{3},x,varargin{:});

Caused by:
Failure in initial user-supplied objective function
evaluation. FSOLVE cannot continue.

请大家帮帮忙了！万分感谢

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有220人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

在matlab中拟合曲面的表达式已经有29人回复
非线性拟合估计方程参数（matlab，spss，origin都可以）已经有4人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
matlab非线性参数拟合问题已经有7人回复
请教Matlab高手解多元非线性方程组中的参数已经有6人回复
非线性拟合用哪些软件？已经有8人回复
【求助】求MAtlab求解程序！！！已经有20人回复
【求助】求教matlab解非线性方程组已经有9人回复

1楼 2011-11-24 11:04:43

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

wtisunday(金币+20): 呵呵，整明白了，主要是fsolve只能解决未知量为向量的 2011-11-24 17:24:29

CODE:

function AA

x=[1 100 1000 10000 100000,20000,900,20];

fsolve(@fun,x)

function ff=fun(x)

w=[0.00100 0.02000 0.1300 0.910000 8.7000 76 1001 4301 20032 100000]';

G1=[1000 2034 3000 4700 6002 10032 50345 109820 620000 1200000]';

G2=[6000 7200 9000 11000 13000 23000 38000 70012 209903 632098]';

for i=5:8

    y(mod(i,5)+1)=x(i);

end

for i=1:4

    for j=1:4

        A(i,j)=0,B(i,j)=0,b(i)=0,c(i)=0

        for k=1:10

            A(i,j)=A(i,j)+(w(k)*x(j))^2*(w(k)*x(i))^2/(G1(k))^2/(1+(w(k)*...

                x(j))^2)/(1+(w(k)*x(i))^2);

            B(i,j)=B(i,j)+2*y(i)*w(k)^2*x(i)*(w(k)*x(j))^2/(G1(k)^2)/(1+...

                (w(k)*x(i))^2)^2/(1+(w(k)*x(j))^2);

            b(i)=b(i)+(w(k)*x(i))^2/G1(k)/(1+(w(k)*x(i))^2);

            c(i)=c(i)+2*y(i)*x(i)*(w(k))^2/G1(k)/(1+(w(k)*x(i))^2);

        end

    end

end

b=b';

c=c';

X=x(1:4)';

Y=x(5:8)';

ff1=A*X-b;

ff2=B*Y-c;

ff=[ff1;ff2];