《Wavelet Theory：An elementary Approach With Applications》中关于例6.10的疑惑

我用matlab仿真式(6.79)，式中a0,k直接取书上算出来的系数，即[0.006,-0.056,2.347,3.192,2.009,-0.663]; φ (t)则直接由书上的滤波器系数求得，结果仿真的结果跟书上相差甚远，百思不解。

仿真结果与书上结果比较：

我的结果

书上的结果

现将源码公布如下：

clear clc tic t=-3:0.001:3; a0=[0.006,-0.056,2.347,3.192,2.009,-0.663]; h4=[1+sqrt(3),3+sqrt(3),3-sqrt(3),1-sqrt(3)]/(4*sqrt(2)); jinsi = zeros(1,length(t)); for k=-4:1     jinsi = jinsi + a0(k+5)*phnt(t-k,h4,5);  %phnt是我另写的函数，递归求daubechies尺度函数的5阶逼近 end hold on;PlotAxisAtOrigin(t,jinsi); toc

function [ output ] = phnt( t,h,n ) % t:时间范围 % n:迭代次数 % h:尺度滤波器系数     output = zeros(1,length(t));     if(n==1)         for k=1:length(h)             output = output + h(k)*B0(2*t-k);