[求助] 拟合式中有sinθ，图形怪异求帮助。。。。泪。。。。

Fx=load('Fxx.txt');
Fy=load('Fyy.txt');

n=length(Fx);                %Fx,Fy数据数量
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;
delta=70;
for k=1+delta:n-delta
%   Kt(k)=Ft(k)/(ft*ap);     %
    Kt(k)=Ft(k)/(ap*ft*sin(theta(k)));  
end
E=10;                                             %设定多项式拟合次数
a=polyfit(theta,Kt,E);                            %拟合多项式
Kt_fit=polyval(a,theta);
for k=1+delta:n-delta
%   Ft_fit(k)=ap*ft*Kt_fit(k);      %
    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
plot(theta,Ft_fit,'g')      % Ft_fit
hold on
grid on
legend('Ft','Ft-fit')

2楼: Originally posted by taokaijoshua at 2014-12-02 11:11:36
顶顶

Fx=load('Fxx.txt');
Fy=load('Fyy.txt');

n     =length(Fx);                                     % Sample quantity of Fx,Fy
theta =zeros(n,1);                                     % Rotation angle
Ft    =zeros(n,1);                                     % Tangential force initialization

for k=1:n
    theta(k)=pi/n*k;                                   % Discretize the rotation angle
end
for k=1:n
    Ft(k)=-cos(theta(k))*Fx(k)+sin(theta(k))*Fy(k);    % Evaluate Ft
end

% Calculate Kt,Kr
ft          =0.1;                                      % Feed per tooth
ap          =3;                                        % Axial cutting depth
delta_theta =10;                                       % Set the angle-cut length at both ends (Unit: degree)
delta       =round(delta_theta/180*n);                 % Quantity of samples deserted   
theta_Mid   =zeros(n-2*delta,1);                       % Angle data for fitting initialization
Kt          =zeros(n-2*delta,1);                       % Coefficient of tangential force
Ft_fit      =zeros(n-2*delta,1);                       % Fitting of Tangential force initialization


for k=delta+1:n-delta
%   Kt(k)=Ft(k)/(ft*ap);      
    Kt(k)=Ft(k)/(ap*ft*sin(theta(k)));                 % Evaluate Kt within the given interval
end

for k=1:n-2*delta
   theta_Mid(k)=theta(k+delta);                        % Cut angle data near 0 and pi
end

Degree=5;                                              % Set the polynomial degree
p_Ft=polyfit(theta_Mid,Kt,Degree);                     % Evaluate the coefficients
Kt_fit=polyval(p_Ft,theta_Mid);

for k=delta+1:n-delta
%   Ft_fit(k)=ap*ft*Kt_fit(k);      %
    Ft_fit(k)=ap*ft*sin(theta(k))*Kt_fit(k);           % Fitting expression of Ft
end

2楼: Originally posted by taokaijoshua at 2014-12-02 11:11:36
顶顶