求分段三次Hermite插值函数
发布网友
共1个回答
热心网友
n=-1:3; %n=-1,0,1,2,3
x=1:5; %x=1,2,....,5
k=0:500; %k=0,1,.....,500
w=(pi/500)*k; %w=pi/500*k，pi是3.1415926
X=x*(exp(-j*pi/500)).^(n'*k); %(.^)中点的意思是元素操作，^是次方， n' 是吧n置换 ，j是虚数
magX=abs(X); %abs(x)就是数学中的|X|
angX=angle(X); %angle(X)是找角度的
subplot(2,1,1) %图片排列想成矩阵，那么就有2行一列，这个图是第一个位置
plot(w./pi,magX); %画图(x轴，y轴）
title('幅度响应'); %标题
grid; %加格子
ylabel('幅度'); %y轴标题
xlabel('以\pi为单位的频率'); %x轴标题
subplot(2,1,2) %第二张图
plot(w/pi,angX)
title('相位响应');
grid;
ylabel('相位/\pi');
xlabel('以\pi为单位的频率');
