[gy022]

function [y]=W(x)
% 本函数为一窗口函数
% T0为单齿作用产生的转矩
% x=0-2*pi 摩擦齿轮的旋转角度，w是设计变量
%c2=w3; %摩擦齿轮基圆半径
%c3=w2; %齿顶高度
w=0:0.01:2*pi;
c4=0.371; %轮胎外径
c10=2; %公式M中的常系数
e1=2; %公式M中的常指数
w0=pi/2-acos(((x(3)+c4)^e1+(x(3)+x(2))^e1-c4^e1)/(c10*(x(3)+c4)*(w(3)+w(2))));%单齿旋转初始角
w1=pi-w0;%单齿旋转结束角度
%a1=w1;%轮齿的齿侧面和齿顶面的夹角，也是摩擦轮胎的切入角
z=zeros(size(w));
N=length(w);
for k=1:N
if w(k)<w0;
z(k)=0;
elseif w(k)>=w0&&w(k)<=w1;
z(k)=1;
else w(k)>w1;
z(k)=0;
end
end
y=z;