最小二乘法求多元目标函数的极值
本帖最后由 jgn824 于 2015-5-31 18:42 编辑想给定一组初始值，然后求使中间的式子等于0的一组函数值。不知道真么求。目的是求解使式子一成立的方程组的解。这是一个6-sps求数值正解的过程。谁给指点一下。
主程序
tic
clear
pack
clc
X=fsolve(myfun,)%给定初始位置(X0,Y0,Z0,alpha,beta,gamma)
调用函数
function f=myfun(X0,Y0,Z0,alpha,beta,gamma)
B1=[-1145,-1145,0];B2=;B3=;B4=;B5=[-1145,1145,0];B6=[-2055,0,0]; %定平台6个铰链点坐标
P1=[-907,-1045,0];P2=;P3=;P4=;P5=[-907,1045,0];P6=[-1703,0,0]; %定平台6个铰链点坐标
L1=950;L2=950;L3=1026;L4=950;L5=950;L6=1026;%指定的杆长
T1=;%转换矩阵【T1;T2;T3】
T2=;
T3=[-sin(beta),cos(beta)*sin(alpha),cos(beta)*cos(alpha)];
%运用公式P定=【T1;T2;T3】*P动+P移 将动坐标系中的P坐标转化到定坐标系之中
Po1=;%P1点在定坐标系中的坐标
Po2=;%P2点在定坐标系中的坐标
Po3=;%P3点在定坐标系中的坐标
Po4=;%P4点在定坐标系中的坐标
Po5=;%P5点在定坐标系中的坐标
Po6=;%P6点在定坐标系中的坐标
f=(L1-sqrt(dot((Po1-B1),(Po1-B1))))^2+(L2-sqrt(dot((Po2-B2),(Po2-B2))))^2+(L3-sqrt(dot((Po3-B3),(Po3-B3))))^2+(L4-sqrt(dot((Po4-B4),(Po4-B4))))^2+(L5-sqrt(dot((Po5-B5),(Po5-B5))))^2+(L6-sqrt(dot((Po6-B6),(Po6-B6))))^2;
提示的错误是gamma未定义！！！！
页:
[1]