求助大侠指点一下关于fsolve的求解问题
function yueshu01global L n
L=input('请输入杆件长度矩阵L='); n=length(L); q0=rand(n,1);
F=fsolve(@myfun2,q0)
function F=myfun2(q)
global L n
for i=1:n, XX(i)=X(i)-L(i)*cos(q(i)); YY(i)=Y(i)-L(i)*sin(q(i)) ; end
F=
出错信息是X，Y没定义。这里的X，Y该怎样定义或可以用其他的什么办法解决呀？我用syms X Y 这样定义了一下但是长度不够。还恳请各位大哥大侠帮帮忙啊！！X(i),Y(i)都是关于q(i)的变量也是一个要求的未知量 q(i)也是要求的未知量
[ 本帖最后由 ChaChing 于 2009-5-2 10:41 编辑 ] 没人回复我自己顶一个。呵呵把XX(i)=X(i)-L(i)*cos(q(i)) YY(i)=Y(i)-L(i)*sin(q(i))换成 X(i)=(i)*cos(q(i)) -0Y(i)=(i)*sin(q(i))-0；把F=换成F=就可以了 呵呵
求助：关于fsolve的循环求解的问题
各位大侠，小弟现在有一关于fsolve的循环求解的问题求助，还请各位大家多多指教，小弟先行拜谢了！现给出如下代码（注：代码是可以运行的，只是无法循环计算）。我要求能够分别输出关于时间t的各时刻的结果，但是我这样做了后，只能多次循环输出时间t的最后一组结果。还请各位指点一二，小弟感激不尽！！function F=yueshufangcheng(q,t)
L=load('ganchang.txt');
n=length(L);
syms t
T=load('shijian.txt');
N=length(T);
load yizhijiao.mat;
o=length(U);%U为独立坐标变量%
%===========================================
u=subs(U,t,T);
for i=1:N
MD=q(2*n+1:2*n+o)-u(i);%附加驱动约束%
end
%u是随时间t变化的量，因此，MD也会随着时间而变化
%==============================================
if n==1
MK(1)=q(1)-1/2*L(1)*cos(q(2*n+1));
MK(2)=q(n+1)-1/2*L(1)*sin(q(2*n+1));
elseif n>=2
MK(1)=q(1)-1/2*L(1)*cos(q(2*n+1));
MK(2)=q(n+1)-1/2*L(1)*sin(q(2*n+1));
for i=2:n
MK(2*i-1)=q(i)-1/2*L(i)*cos(q(2*n+i))-q(i-1)-1/2*L(i-1)*cos(q(2*n+i-1));
MK(2*i)=q(n+i)-1/2*L(i)*sin(q(2*n+i))-q(n+i-1)-1/2*L(i-1)*sin(q(2*n+i-1));
end
MK(2*n+1)=0.4-q(n)-1/2*L(n)*cos(q(3*n));
MK(2*n+2)=0.1-q(2*n)-1/2*L(n)*sin(q(3*n));
end
%=====================================
F=;
%MD会随着时间变化
%=====================================
end
function weixingjisuan01
xr0=;
yr0=;
m0=;
q0=;
options=optimset('Display','iter');
T=load('shijian.txt');
N=length(T);
for i=1:N
q=fsolve(@yueshufangcheng01,q0,options)
end
save zuobiao.txt q -ASCII
相关数据如下：
杆长L=
时间T=
已知角U=
再次拜谢各位啦！！ 很着急结果呀！求求那位高手帮忙一下，小弟感激不尽啦！！
页:
[1]