matlab傅里叶变换急求！

Question

clear %清除内存
clc %清除命令窗口显示内容
x=0:1e-3:20; %以1e-3为单位间隔从0到20生成向量x的元素，离散取值
h=cos(2*x).^2; %将题设赋值给变量h
w=[0:1e-2:2];
y=w;
a=w;
j=sqrt(-1);
for i=1:length(w) %取得w的长度for循环从1开始到length(w)为止
h=(1/sqrt(2*pi))*trapz(x,h.*exp(-j*w(i)*x));%
y(i)=abs(h);
a(i)=angle(h);
end
subplot(3,1,1), %在显示图像窗口中，将窗口划分为上下三块区域（3行1列），在第1个区域显示图像
plot(x,h) %以x为横坐标，以h为纵坐标画出的图形
subplot(3,1,2), %在显示图像窗口中，将窗口划分为上下三块区域（3行1列），在第2个区域显示图像
plot(w,y)
subplot(3,1,3), %在显示图像窗口中，将窗口划分为上下三块区域（3行1列），在第3个区域显示图像
plot(w,a)

求大神看一下我这个程序有什么问题没有。。并且把我没有用%注释的地方注释一下！！感激不尽！！！！！！！！！！！！！！！！！！！！看我真挚的眼神★，★
在线等~~~完成后经验直接50加，不给这辈子没人品！！

Answer 1

clear %清除内存
clc %清除命令窗口显示内容
x=0:1e-3:20; %以1e-3为单位间隔从0到20生成向量x的元素，离散取值
h=cos(2*x).^2; %将题设赋值给变量h
w=0:1e-2:2; %区间[0,2]
y=w; %y赋初值
a=w; %a赋初值
j=sqrt(-1); %虚数单位
for i=1:length(w) %取得w的长度for循环从1开始到length(w)为止
g=(1/sqrt(2*pi))*trapz(x,h.*exp(-j*w(i)*x));%求梯形数值积分
y(i)=abs(g); %求矢径
a(i)=angle(g); %求幅角
end
subplot(3,1,1), %在显示图像窗口中，将窗口划分为上下三块区域（3行1列），在第1个区域显示图像
plot(x,h) %以x为横坐标，以h为纵坐标画出的图形
title('h=(cos(2x))^2') %图名
subplot(3,1,2), %在显示图像窗口中，将窗口划分为上下三块区域（3行1列），在第2个区域显示图像
plot(w,y) %以w为横坐标，以g为纵坐标画出的图形
title('y=abs(g)') %图名
subplot(3,1,3), %在显示图像窗口中，将窗口划分为上下三块区域（3行1列），在第3个区域显示图像
plot(w,a) %以w为横坐标，以a为纵坐标画出的图形
title('a=angle(g)') %图名

追问

虚数单位，矢径，和幅角是什么意思啊？
这样子就是完成了一个函数的傅里叶变换么？是属于哪一种傅里叶变换？

追答

其实，你这个程序没有进行傅里叶变换，只是进行梯形积分，傅里叶变换是用命令fft()

追问

那可以麻烦你写一个fft()的吗？
要变形的函数就是这个h=cos(2*x).^2

可以么~~