请根据你说的步骤写出MATLAB代码

根据你提供的SAR成像参数，以下是一个简单的MATLAB代码示例，用于实现SAR图像重建：

% 读取原始回波数据
data = fread(fid, [2, N], 'uint8'); % 假设原始数据为[N,2]的矩阵，其中第一列为实部，第二列为虚部
data = data(1,:) + 1j*data(2,:); % 将实部和虚部合成复数值

% 距离向压缩
win = kaiser(N, beta); % 数据窗函数
data = data .* win'; % 加窗
data = fft(data); % FFT变换

% 多普勒向压缩
fd = linspace(-0.5, 0.5, N) * prf; % 多普勒频率
f0 = -doppler_center; % 多普勒中心频率
data = data .* exp(-1j*2*pi*fd*(f0/prf)); % 相移
data = fftshift(data); % 频域移位
data = sum(reshape(data, [], num_lines), 2); % 多普勒向累加

% 逆FFT变换
data = ifft(data); % 逆FFT变换

% SAR图像重建
range = linspace(0, c/2/B, N); % 距离向采样点
azimuth = linspace(-pi/2, pi/2, num_lines); % 方位向采样点
[X, Y] = meshgrid(azimuth, range); % 采样网格
data = reshape(data, num_lines, []); % 转换为[num_lines, N/num_lines]的矩阵
figure, imagesc(azimuth, range, abs(data)); % 显示SAR图像

注：以上代码仅为演示用途，实际应用中可能需要更多的错误处理和优化。另外，其中的变量和参数需要根据具体情况进行调整。