, (2.4)
s(t) . - ( 2.4).
function s=cosinobn1(t,Um,T,Uo) % s=cosinobn1(t,Um,T,Uo) % t - % Um - % T - % Uo - | T2=T/2; % dt=abs(mod(t+T2,T)-T2); % % s=cosinob1(dt,Um,T,Uo); |
1.4 -
2.3.3
. , (1.5) :
, (1.5)
(1.6)
(1.6) ( ) MATLAB.
function B=BergN(n,Um,Uo,O)
%
% B=BergN(n,O,Um,Uo)
% n
% O
% Um
% Uo
% B0(O)=(sin(O)-h*cos(O))/pi
% B1(h)=(O-sin(O)*cos(O))/pi
% Bn(h)=2*(sin(n*O)*cos(O)-n*sin(O)*cos(n*O))/(pi*n*(n*n-1))
if nargin == 3
O = acos(Uo/Um);
end
k=length(n);
B=zeros(1,k);
for i=1:k
switch n(i)
case 0, B(i)=Um*(sin(O)-O*cos(O))/pi;
case 1, B(i)=Um*(O-sin(O)*cos(O))/pi;
otherwise
m=n(i);
B(i)=2*Um*(sin(m*O)*cos(O)-...
m*sin(O)*cos(m*O))/(pi*m*(m*m-1));
end
end
MATLAB
n = 0:10; a = BergN(n,2,1);
, n = 010:
0.2180 0.3910 0.2757 0.1378 0.0276 -0.0276 -0.0315 -0.0098 0.0098 0.0138 0.0050.
( 1.5), stem(n,a*T). . .
1.5 -