有没有大佬真的我这个程序运行总是出错，该怎么改啊？

朱威斯

% 发射5MHz的超声波信号
fs = 50e6; % 采样率为50MHz
t = 0:1/fs:10e-6; % 发射时长为10微秒
f0 = 5e6; % 中心频率为5MHz
s = sin(2pif0*t); % 生成超声波信号

% 加入噪声
snr = 30; % 信噪比为30dB
noise = randn(size(s));
noise = noise/norm(noise)*norm(s)/10^(snr/20);
x = s + noise;

% 接收信号并进行放大分析
y = amplification(x); % 放大分析函数
plot(y); % 绘制放大后的信号图像

% 快速傅里叶变换
N = 1024; % 数据点数
Y = fft(y,N); % 进行快速傅里叶变换

% 对探头中心频率的幅值取对数
f = fs*(0:(N/2))/N; % 构造频率坐标
Ylog = 20*log10(abs(Y(1:N/2+1))/max(abs(Y)));
% 最小二乘拟合