时变信号频域极值三维绘图

问题遇到的现象和发生背景

本人matlab入门，学校下方任务中有一道练习题，关于时变信号频域极值三维绘图，但是本人所学内容尚未涉及题目方面内容

问题相关代码，请勿粘贴截图

clc
clear
close all
fs=10000;
T=1/fs;
N=500;
t=(0:N-1)T;
t1=0:0.001:0.5;
y1=cos(250pi.t1)+(1/5).cos(2250pi.t1)+(1/7).cos(2350pi.t1)+randn(1)/5;
t2=0.5:0.001:1;
y2=cos(250pi.t2)+(1/7).cos(2350pi.t2)+(1/11).cos(2550pi.t2)+randn(1)/5;
figure
plot(y1);
xlabel("y1加白噪声后原始信号");
figure
plot(y2);
xlabel("y2加白噪声后原始信号");
y1=fft(y1);
y2=fft(y2);
P1_2 = abs(y1/N);
P1_1 = P1_2(1:N/2+1);
P1_1(2:end-1) = 2P1_1(2:end-1)；
P2_2 = abs(y2/N);
P2_1 = P2_2(1:N/2+1);
P2_1(2:end-1) = 2P2_1(2:end-1);
f = fs(0:(N/2))/N;
figure
plot(f,P1_1);
xlabel("y1白噪声后fft");
figure
plot(P2_1);
xlabel("y2白噪声后fft");

运行结果及报错内容

运行结果是二维图，不是三维

我的解答思路和尝试过的方法

无解答思路，对不起，由于本人水平实在有限，面对题目实在不知道从何入手，下附题目图片

我想要达到的结果

完成题目，提供一些入手思路

您好，请用“插入代码块”重新提交下代码，乘号*都被当做加粗省略了

时域转频域不是要傅里叶变换吗