matlab程序运行中特别慢,程序运行很慢，怎么优化？

clear all;

close all;

clc;

format long;

tic

%%

fs=10e9;%采样率设置

fr_int=1e3;%频率分辨率，fft时的频率分辨率，N=fs/fr_int，得到总点数，同时得到采样的总时间，N/fs

t0=40e-6;%激励持续时间

N0=ceil(fs*t0)-1;%激励点数

N=fs/fr_int;%总点数

t_total=N/fs;%采样持续时间

f_sig=433.92e6;%16.493422089129357 + 3.193299237956759i,16.8

d_t=1/fs;

k=1:N;

y=sin(f_sig*2*pi*k*d_t);%总共发射信号时间

y=[y(1:N0),zeros(1,N-N0)];%激励信号加窗，t0长度

y_f=fft(y);%对激励信号fft，获得sinc函数

length_y_f=length(y_f)/2;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%    等效电路模拟导纳曲线

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

R=16;L=98.501e-6;C=1.3658e-15;C0=1.5e-12;

ffr=2*pi*sqrt(L*C);

ffr=1/ffr;%谐振器谐振频率

Qr=(L/C)^(1/2)/R;%谐振器