发现这段代码的问题并改正

根据你的指示，我来找出并修改代码中的问题。

问题1：在计算 x 的时候，x1 的部分被重复计算了。

修改1：删除重复计算的部分，即删除以下代码：

x1=2.3*sin(2*pi*(0:L-1)'*2*f0/fs) + 1*sin(2*pi*(0:L-1)'*3*f0/fs);
x1=x1-mean(x1);

问题2：plot 函数调用时传入的向量长度不相同。

修改2：将 plot 函数的参数改为 c 的长度，即使用 `plot(c)`。

修改后的代码如下：

fs = 2e3;               % 采样频率 (单位：Hz)
L = 2e4;                % 信号长度 (样本数)
f0 = 25;                % 周期频率 (单位：Hz)

x1 = randn(L,1);        % 随机信号

x = x1.*(3 + 2.3*sin(2*pi*(0:L-1)'*f0/fs)+ 2.3*sin(2*pi*(0:L-1)'*2*f0/fs) + 1*sin(2*pi*(0:L-1)'*3*f0/fs));
x = x - mean(x);

c = xcorr(x,x1);
f = (0:length(c)-1) * fs/L;

plot(f,c)

请注意，我还根据修改后的代码重新计算了频率向量 f 的长度，以保证与相关向量 c 的长度一致。这样就可以正确地绘制频谱图了。