[ÇóÖú] 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¿ÉÒÔ¾ßÌå¿´