matlab正弦波合成方波

可以使用傅里叶级数进行正弦波合成方波。

步骤如下：

1. 定义方波周期 T 和采样点数 N。
2. 定义傅里叶级数中的基频率 f0，即第一个谐波的频率。
3. 计算每个谐波的频率 fn = nf0，其中 n 为谐波次数。
4. 计算每个谐波的振幅 An，公式为 An = 2/(n*pi) * (1 - cos(n*pi)).
5. 对于每个采样点，计算出该点的函数值。对于第 i 个采样点，函数值为：

f(i) = A1*sin(2*pi*f1*t(i)) + A2*sin(2*pi*f2*t(i)) + ... + AN*sin(2*pi*fn*t(i))

其中 t(i) 为第 i 个采样点的时间，0 <= t(i) <= T。

6. 绘制合成后的正弦波方波。

下面是 MATLAB 代码示例：

% 定义方波周期 T 和采样点数 N
T = 1;
N = 1000;

% 定义傅里叶级数中的基频率 f0
f0 = 1/T;

% 计算每个谐波的频率和振幅
n = 1:2:25;
fn = n*f0;
An = 2./(n*pi).*(1-cos(n*pi));

% 计算每个采样点的函数值
t = linspace(0, T, N);
f = zeros(1, N);
for i = 1:length(n)
    f = f + An(i)*sin(2*pi*fn(i)*t);
end

% 绘制合成后的正弦波方波
plot(t, f, 'LineWidth', 2);
axis([0 T -1.5 1.5]);
title('Synthesized Square Wave');
xlabel('Time (s)');
ylabel('Amplitude');
grid on;

运行该代码可以得到合成后的正弦波方波图像。