调整幅度调制/波形整形

Question

我有一些工作代码可以进行幅度调制并绘制它。 但是我试图改变调制方式（y变量） 所以它看起来像一个蛋形。我找到了一个看起来不错的等式/网站 http://www16.ocn.ne.jp/~akiko-y/Egg/index_egg_E.html 但我不确定如何将其转换为matlab / octave代码来更改y变量

%test_amplitude modultaion
fs=1000;
t=linspace(0,2*pi,fs);

mt=1*sin(100*t); %signal you want to use
y=mt.*(1+cos(1*t+pi));%modulation equation, use pi to shift over 90 deg to start at 0

y=y';
y_norm=(y(:,1)/max(abs(y(:,1)))*.8); %normalize signal

plot(y_norm)

PS：这是matlab / octave代码

Answer 1

使用您链接的页面上给出的公式：

a = 2*pi;
b = a; % change this depending on the shape of the egg you want

mt=1*sin(100*t); %signal you want to use
y = mt.*sqrt((a-b)-2*t + sqrt(4*b*t + (a-b)^2)).*sqrt(t)/sqrt(2); % modulation

你的其余代码是A-OK，虽然最后我可能会使用plot(t,y_norm)。