一直就想好好学习下Mablab的，不过总给各种各样的琐事打扰，现在好了，有整整一个暑假可以研习这个经典的软件了，图书馆借了本书，翻开就有点恐怖的架式，一堆堆似懂非懂的公式和框图，不管了，先看了再说。依我的能力看懂应该是没有问题的哈，嘻嘻。
下面就开始咯，先为自己的这门自修课来道开菜：上采样和内插
当然咯，算法么还是没有研究过得啦，不过是拿采样来练练手的。
对标准sin函数的采样，x=sin(2*pi*t/10)
采样主程序：
1
% File: c3_upsampex.m
2M=6;
3h=c3_lininterp(M);
4t=0:10;
5tu=0:60;
6x=sin(2*pi*t/10);
7xu=c3_upsamp(x,M);
8subplot(3,1,1)
9stem(t,x,'k.')
10ylabel('X')
11subplot(3,1,2)
12stem(tu,xu,'k.')
13ylabel('Xu')
14xi=conv(h,xu);
15subplot(3,1,3)
16stem(xi,'k.')
17ylabel('Xi')
% File: c3_upsampex.m2
M=6;3
h=c3_lininterp(M);4
t=0:10;5
tu=0:60;6
x=sin(2*pi*t/10);7
xu=c3_upsamp(x,M);8
subplot(3,1,1)9
stem(t,x,'k.')10
ylabel('X')11
subplot(3,1,2)12
stem(tu,xu,'k.')13
ylabel('Xu')14
xi=conv(h,xu);15
subplot(3,1,3)16
stem(xi,'k.')17
ylabel('Xi')上采样函数：
1% File: c3_upsample.m
2function out=c3_upsamp(in,M)
3L=length(in);
4out=zeros(1,(L-1)*M+1);
5for j=1:L
6 out(M*(j-1)+1)=in(j);
7end
8% End of function file.
% File: c3_upsample.m2
function out=c3_upsamp(in,M)3
L=length(in);4
out=zeros(1,(L-1)*M+1);5
for j=1:L6
out(M*(j-1)+1)=in(j);7
end8
% End of function file.线性内插函数：
1% File: c3_lininterp.m
2function h=c3_lininterp(M)
3h1=zeros(1,(M-1));
4for j=1:(M-1)
5 h1(j)=j/M;
6end
7h=[0,h1,1,fliplr(h1),0]
8% End of function file.
% File: c3_lininterp.m2
function h=c3_lininterp(M)3
h1=zeros(1,(M-1));4
for j=1:(M-1)5
h1(j)=j/M;6
end7
h=[0,h1,1,fliplr(h1),0]8
% End of function file.生成的图形：
不错，很能说明问题了，上采样的开销较大，线性内插的就好多了的说，嘻嘻，到这里结束^_^
