n)
附录matlab 代码:
%%%%%%%%%%%%%%%%%% define the condition
N=200; %%define calculating time
delta=0.4; %%define the step
sn=sqrt(0.05)*randn(1,N);
n=1:N;
N0=sin(1/8*pi*n+pi/10);
N1=sqrt(2)*sin(1/8*pi*n);
input_x=N1;
reference_y=N0+sn;
%%%%%%"误差性能曲面"
x = -2:0.05:4;
y = -4:0.05:2;
[h0,h1] = meshgrid(x,y);
Jn=(0.55 + h0.^2 + h1.^2
+2*h0.*h1.*cos(pi/8)-sqrt(2)*(h0.*cos(pi/10)+h1.*cos(9*pi/40)));
figure(1);
subplot(1,2,1);
surf(h0,h1,Jn);
axis([-2 4 -4 2 0 40]);
xlabel('h0');
ylabel('h1');
zlabel('Jn')
title('误差性能曲面');
subplot(1,2,2);
v=0:0.05:2.5;
contour(x,y,Jn,v);
axis([-2 4 -4 2]);
title('误差等值线');
xlabel('h0');
ylabe
附录matlab 代码:
%%%%%%%%%%%%%%%%%% define the condition
N=200; %%define calculating time
delta=0.4; %%define the step
sn=sqrt(0.05)*randn(1,N);
n=1:N;
N0=sin(1/8*pi*n+pi/10);
N1=sqrt(2)*sin(1/8*pi*n);
input_x=N1;
reference_y=N0+sn;
%%%%%%"误差性能曲面"
x = -2:0.05:4;
y = -4:0.05:2;
[h0,h1] = meshgrid(x,y);
Jn=(0.55 + h0.^2 + h1.^2
+2*h0.*h1.*cos(pi/8)-sqrt(2)*(h0.*cos(pi/10)+h1.*cos(9*pi/40)));
figure(1);
subplot(1,2,1);
surf(h0,h1,Jn);
axis([-2 4 -4 2 0 40]);
xlabel('h0');
ylabel('h1');
zlabel('Jn')
title('误差性能曲面');
subplot(1,2,2);
v=0:0.05:2.5;
contour(x,y,Jn,v);
axis([-2 4 -4 2]);
title('误差等值线');
xlabel('h0');
ylabe