% Parameters
omega_a = 0.1*pi;
omega_b = 0.05*pi;
omega_I = 0.5*pi;
omega_Q = 0.5*pi;
fs = 1000; % Sampling frequency
T = 1/fs; % Sampling period
t = 0:200; % Time vector
% Generate input signals
x1 = cos(omega_a*t);
x2 = sin(omega_b*t);
% Modulation
modulated_signal = x1 .* cos(omega_I*t) - x2 .* sin(omega_Q*t);
% Demodulation
demodulated_x1 = modulated_signal .* cos(omega_I*t);
demodulated_x2 = -modulated_signal .* sin(omega_Q*t);
% Butterworth lowpass filter design
cutoff_freq = 0.15*pi; % Adjust as needed
order = 4;
[b, a] = butter(order, cutoff_freq);
% Apply filter to regenerate input message signals
filtered_x1 = filter(b, a, demodulated_x1);
filtered_x2 = filter(b, a, demodulated_x2);
% Plot the results
figure;
subplot(1,1,1);
plot(t, x1,'r',t, x2,'g');
title('Original x1[n] and x2[n]');
xlabel('Time');
ylabel('Amplitude');
legend('x1[n]', 'x2[n]');
figure;
subplot(1,1,1);
plot(t, modulated_signal);
title('Modulated Signal');
xlabel('Time');
ylabel('Amplitude');
figure;
subplot(1,1,1);
plot(t, filtered_x1, 'r', t, filtered_x2, 'g');
title('Filtered Demodulated Signals');
xlabel('Time');
ylabel('Amplitude');
legend('Filtered x1[n]', 'Filtered x2[n]');
figure;
subplot(2,1,1);
plot(t, demodulated_x1);
title('demodulated x1');
xlabel('Time');
ylabel('Amplitude');
2023-07-28
最后編輯于 ：
?著作權歸作者所有,轉載或內容合作請聯系作者
- 文/潘曉璐 我一進店門，熙熙樓的掌柜王于貴愁眉苦臉地迎上來有勾，“玉大人疹启，你說我怎么就攤上這事“ǎ” “怎么了喊崖？”我有些...
- 正文 為了忘掉前任节仿，我火速辦了婚禮，結果婚禮上掉蔬，老公的妹妹穿的比我還像新娘廊宪。我一直安慰自己，他們只是感情好女轿，可當我...
- 文/花漫 我一把揭開白布挤忙。 她就那樣靜靜地躺著，像睡著了一般谈喳。 火紅的嫁衣襯著肌膚如雪册烈。 梳的紋絲不亂的頭發(fā)上，一...
- 文/蒼蘭香墨 我猛地睜開眼驾中，長吁一口氣：“原來是場噩夢啊……” “哼唉堪！你這毒婦竟也來了？” 一聲冷哼從身側響起肩民，我...
- 正文 年R本政府宣布，位于F島的核電站攻人，受9級特大地震影響取试，放射性物質發(fā)生泄漏。R本人自食惡果不足惜怀吻，卻給世界環(huán)境...
- 文/蒙蒙 一瞬浓、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧蓬坡，春花似錦猿棉、人聲如沸。這莊子的主人今日做“春日...
- 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至兆龙，卻和暖如春杖爽，著一層夾襖步出監(jiān)牢的瞬間，已是汗流浹背紫皇。 一陣腳步聲響...
推薦閱讀更多精彩內容
- 好好生活赶盔，讓自己知足常樂企锌，如此，才能享受生活的幸福于未。我們都知道撕攒，人生沒有重來時，歲月也不會為誰而停留烘浦，人生沒有如果...
