[求助] sinx要做分母，怎么解决0点数值无穷大的问题

Fx=load('Fxx.txt');  %Fxx.txt 在附件中
Fy=load('Fyy.txt');  %Fyy.txt 在附件中

n=length(Fx);
theta=zeros(n,1);
Ft=zeros(n,1);
Kt=zeros(n,1);
Ft_fit=zeros(n,1);

for k=1:n
    theta(k)=pi/n*k;
end
for k=1:n
    Ft(k)=-cos(theta(k))*Fx(k)+sin(theta(k))*Fy(k);
end
% Calculate Kt,Kr
ft=0.1;ap=3;
for k=1+30:n-30
    Kt(k)=Ft(k)/(ap*ft*sin(theta(k)));  
end
E=5;                                              %设定拟合次数
a=polyfit(theta,Kt,E);                            %拟合多项式
Kt_fit=polyval(a,theta);
for k=1+30:n-30
    Ft_fit(k)=ap*ft*sin(theta(k))*Kt_fit(k);      %Ft的拟合表达式
end

% 结果输出
%-------------------------------------------------------------------------%
figure
set(gcf,'color','w')
plot(theta,Fx,'k')           % Fx
hold on
plot(theta,Fy,'b')           % Fy
grid on
legend('Fx,Fy')
%-------------------------------------------------------------------------%
figure
set(gcf,'color','w')
plot(theta,Ft,'r')          % Ft
hold on
grid on
legend('Ft')
%-------------------------------------------------------------------------%
figure
set(gcf,'color','w')
hold on
plot(theta,Ft_fit,'g')      % Ft_fit
grid on
legend('Ft-fit')
% End
   

5楼: Originally posted by 一段木头 at 2015-01-09 00:44:07
做分母时+eps，输入 help eps可以具体看