%%修改后的带注释
clc;
clear;
close all;
warning off;
% addpath(genpath(pwd)); % 如果你不需要改变路径,可以注释掉这一行
% 三维图SBS温度变化
z1 = 0:30:3000; % 第一个z轴区间
z2 = 3000:5:3500; % 第二个z轴区间
z3 = 3500:5:4000; % 第三个z轴区间
f = 10000:1:12000; % 频率范围
fB1 = 11200; % 第一个布里渊频移
fB2 = 11260; % 第二个布里渊频移
fB3 = 11200; % 第三个布里渊频移
fBi = 35; % 布里渊线宽
W = 150e-9; % 波束宽度
P