倍频程：绘制离散信号，奈奎斯特采样问题

如何解决倍频程：绘制离散信号，奈奎斯特采样问题

我在Octave上创建了一个程序，目的是在其中创建一个带有噪声的连续时间信号，然后应该对该噪声进行离散采样（显示为正下方的橙色线适合正弦波）。这应该是完成代码后的预期输出。

但是，尽管我确实设法获得了相同的精确波，但我没有得到的是连续时间信号的离散采样，因此得到了它。

由于我对它们应用了相同的奈奎斯特定理，因此我不确定这些线为何不适合图形，下面的代码将更好地解释我的意思。

%Time Base
t = 0:0.001:1.8;

%Nyquist Frequencies
Fn1 = 1;
Fn2 = 6;

%Nyquist Rates 
Fnr1 = 2*(Fn1);
Fnr2 = 2*(Fn2);

%Sampling Period
Sp1 = 5*(Fnr1);
Sp2 = 5*(Fnr2);

Ts1 = 1/(Sp1);
Ts2 = 1/(Sp2);
T1 = 1/(Fn1);
T2 = 1/(Fn2);

%Number of Samples
N1 = (T1/Ts1); 
n1 = 0:1:N1-1; 
N2 = (T2/Ts2); 
n2 = 0:1:N2-1; 

nTs1 = n1 * Ts1; 
nTs2 = n2 * Ts2; 

x_c = sin(2*pi*Fn1*nTs1);
x_c2 = sin(2*pi*Fn2*nTs2);

x_c1 = sin(2*pi*Fn1*t);
x_2 = sin(2*pi*Fn2*t);

signal = x_c1 + x_2;
ct = nTs1 + nTs2;
nqsignal = x_c + x_c2;

%Second Part
h = stem(ct,nqsignal,'linewidth',2);
hold
plot(t,signal,2)
lgd = legend('Discrete Data','Continuous Data');
set (lgd,"fontsize",12)
set(gca,'XTick',[0:0.2:1.8]);
set(gca,'YTick',[-2:0.5:2]);
title('Time vs Magnitude','fontweight','bold','fontsize',16);
xlabel('Time(s)',14)
ylabel('Magnitude',14)
grid

我只想知道哪里出了问题以及如何解决这种错误。