調整幅度調製/波形整形

Question

我有一些可以進行幅度調製並繪製它的工作代碼。 但是我試圖改變調製的方式看起來像（y變量） 所以它看起來像一個雞蛋形狀。我發現了一個公式/網站看起來不錯 http://www16.ocn.ne.jp/~akiko-y/Egg/index_egg_E.html 但我不知道如何將其轉換成MATLAB /倍頻代碼改變變量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 /倍頻代碼

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)。