自己编写卷积程序 matlab

function y=myconv（hx） 
n=length（x）;  
k=length（h）;      %定义序列长度 
for m=1:1:（n-1）       
    A（m）=0;        
end
for m=n:1:n+k-1      
    A（m）=h（m-n+1）; %给h序列添（n-1）个0 
end
B=hankel（A）;
    for i=1:1:n      
        C（:i）=B（:i）;%抽取B矩阵中从第1列到第n列的矩阵
    end
    D=flipl