Question

求助(matlab)：求助各位高手，我用fsolve函数解二元二次方程组，解100次for循环中，为何每次值都不一样？

求助(matlab)：求助各位高手，我用fsolve函数解二元二次方程组，解100次for循环中，为何每次值都不一样？而且相差而大。不胜感激。如有疑问，下附循环程序，谢谢！
for k=1:100
xinhao=sa*cos(w*t)+na*randn(size(t));%信号
wg=wgn(1,1901,1);
x0=[zeros(1,320) sa0.*cos(w.*(t)) zeros(1,1500)]+wg(k);
% % x0=sa1*cos(w*t)+n*randn(size(t));
x1=[zeros(1,640) sa1.*cos(w.*(t)) zeros(1,1180)]+wg(k);
x2=[zeros(1,960) sa2.*cos(w.*(t)) zeros(1,860)]+wg(k);
% figure(1);
% plot(x0);

% x1=sa2*cos(w*(t-t10))+n*randn(size(t));
% x2=sa3*cos(w*(t-t20))+n*randn(size(t));

% x1=[zeros(1,4)+n.*randn(4) sa3.*cos(w.*(t))+n.*randn(size(t)) zeros(1,20)+n.*randn(20)];
% x2=[zeros(1,2)+n.*randn(2) sa3*cos(w.*(t))+n.*randn(size(t)) zeros(1,22)+n.*randn(22)];
% wg=wgn(1,225,1);
% x1=[zeros(1,4) sa1.*cos(w.*(t)) zeros(1,140)]+wg;
% x2=[zeros(1,122) sa2.*cos(w.*(t)) zeros(1,22)]+wg;
% figure(2);
% plot(x1);
% figure(3);
% plot(x2);
%plot(x2);
c=xcorr(xinhao,x0);
c1=xcorr(xinhao,x1);
c2=xcorr(xinhao,x2);
% c=xcorr(x0,x1);
% c1=xcorr(x0,x2);
%x=[zeros(1,length(x1)) x2 zeros(1,length(x1))]
%for k=length(x1):length(x)
%y(k)=x(k-length(x1)+1:k)*x1';
%end
% figure(4);
% plot(c);
% figure(5);
% plot(c1);
%figure(6);
%plot(y);
% figure(6);
% plot(c2);
cdianshu=find(c==max(c)) %延迟
c1dianshu=find(c1==max(c1))
c2dianshu=find(c2==max(c2))
t10(k)=abs(cdianshu-c1dianshu)/fs
t20(k)=abs(cdianshu-c2dianshu)/fs

% syms target_x target_y
% eq1=abs(sqrt((-169.7056-target_x).^2+(-169.7056-target_y).^2)-sqrt((0-target_x).^2+(120-target_y).^2))-c*0.08;
% eq2=abs(sqrt((254.5584-target_x).^2+(-254.5584-target_y).^2)-sqrt((0-target_x).^2+(120-target_y).^2))-c*0.16;
% [target_x target_y]=solve(eq1,eq2,'target_x','target_y');
% function F = myfun(x)
% F = [abs(sqrt((zx1-x(1)).^2+(zy1-x(2)).^2)-sqrt((zx0-x(1)).^2+(zy0-x(2)).^2))-c*t10;
% abs(sqrt((zx2-x(1)).^2+(zy2-x(2)).^2)-sqrt((zx0-x(1)).^2+(zy0-x(2)).^2))-c*t20];
% % z=[target_x target_y];
% numeric(z);
% [target_x,target_y]=solve('abs(sqrt((x1-target_x)^2+(y1-target_y)^2))-sqrt((x0-target_x)^2+(y0-target_y)^2)))=c*t10','abs(sqrt((x2-target_x)^2+(y2-target_y)^2))=sqrt((x0-target_x)^2+(y0-target_y)^2)))=c*t20','target_x','target_y');
x0 = [-5; -5]; % Make a starting guess at the solution
options=optimset('Display','iter'); % Option to display output
[x,fval] = fsolve(@myfun,x0,options) % Call optimizer
target_x(1,k)=x(1);
target_y(1,k)=x(2);
% figure(1)
% plot(x(1),x(2),'.b','Markersize',12);
% hold on
end
我的QQ：541421241
邮箱：1hanjiangbaixue@163.com

我已经找到问题所在了，呵呵。

Answer 1

不好意思 看不懂啊 我是机械设计毕业的
帮不到你