图像处理与分析-王伟强-作业题及答案汇总-2016年版.docx

1、【作业1】1、完成课本习题3.2(a)(b), 课本中文版处理第二版的113页。可以通过matlab帮助你分析理解。a: s=Tr=11+(mr)Eb:E控制函数的斜坡，也就是函数的倾斜程度，E越大，函数倾斜程度越大，如下图1，图2所示：图1：E=5图2：E=202、一幅8灰度级图像具有如下所示的直方图，求直方图均衡后的灰度级和对应概率，并画出均衡后的直方图的示意图。（计算中采用向上取整方法，图中的8个不同灰度级对应的归一化直方图为0.17 0.25 0.21 0.16 0.07 0.08 0.04 0.02）【解答】直方图均衡采用公式Sr=Gw=0rPrw-1式中，G为灰度级数，取8，pr(

2、w)为灰度级w的概率，Sr为变换后的灰度，计算过程如下表所示：灰度级r各级概率Pr(r)累积概率w=0rPrw累积概率81向上取整sr00.170.170.36110.250.422.36320.210.634.04530.160.795.32640.070.865.88650.080.946.52760.040.986.84770.02177则新灰度级的概率分别是：Ps(0) = 0Ps(1) = Pr(0) = 0.17Ps(2) = 0Ps(3) = Pr(1) = 0.25Ps(4) = 0Ps(5) = Pr(2) = 0.21Ps(6) = Pr(3) + Pr(4) = 0.23

3、Ps(7) = Pr(5) = Pr(6) = Pr(7) = 0.14编写matlab程序并绘制直方图：s=0:1:7;p=0 0.17 0 0.25 0 0.21 0.23 0.14;bar(s,p);axis(-1 8 0 0.3);可以看出，此图较题目原图更加“均匀”。【作业2】1、完成课本数字图像处理第二版114页，习题3.10。【解答】由图可知prr=-2r+2,0r1pzz=2z,0z1将两图做直方图均衡变换s1=T1r=0rprwdw=0r-2w+2dw=-r2+2rs2=T2z=0zpzwdw=0z2wdw=z2令上面两式相等，则z2=-r2+2r因为灰度级非负，所以z=-r

4、2+2r2、请计算如下两个向量与矩阵的卷积计算结果。（1） 1 2 3 4 5 4 3 2 1 * 2 0 -2 （2）-101-202-101*1320410323041052321431042【解答】（1）设向量a= 1 2 3 4 5 4 3 2 1 ，下标从-4到4，即a(-4)=1，a(-3)=2a(4)=1；设向量b= 2 0 -2 ，下标从-1到1，即b(-1)=2，b(0)=0，b(1)=-2；设向量c=a*b，下标从-5到5。根据卷积公式可知cx=t=-atbx-t=t=-44atbx-t其中，-5x5，则c(-5)=a(-4)b(-1)=1*2=2c(-4)=a(-4)b(

5、0)+a(-3)b(-1)=1*0+2*2=4c(-3)=a(-4)b(1)+a(-3)b(0)+a(-2)b(-1)=1*(-2)+2*0+3*2=4c(-2)=a(-3)b(1)+a(-2)b(0)+a(-1)b(-1)=2*(-2)+3*0+4*2=4c(-1)=a(-2)b(1)+a(-1)b(0)+a(0)b(-1)=3*(-2)+4*0+5*2=4c(0)=a(-1)b(1)+a(0)b(0)+a(1)b(-1)=4*(-2)+5*0+4*2=0c(1)=a(0)b(1)+a(1)b(0)+a(2)b(-1)=5*(-2)+4*0+3*2=-4c(2)=a(1)b(1)+a(2)b

6、(0)+a(3)b(-1)=4*(-2)+3*0+2*2=-4c(3)=a(2)b(1)+a(3)b(0)+a(4)b(-1)=3*(-2)+2*0+1*2=-4c(4)=a(3)b(1)+a(4)b(0)=2*(-2)+1*0=-4c(5)=a(4)b(1)=1*(-2)=-2所以卷积结果为： 2 4 444 0 -4 -4 -4 -4 -2 （2）设矩阵b=-101-202-101下标从(-1,-1)到(1,1)，即b(-1,-1)=-1，b(-1,0)=0b(1,1)=1；设矩阵a=1320410323041052321431042下标从(-2,-2)到(2,2)，即a(-2,-2)=3

7、，a(-2,-1)=2a(2,2)=4；设矩阵c=a*b=b*a，下标从(-3,-3)到(3,3)。根据卷积公式可知cx,y=s=-t=-as,tbx-s,y-t=s=-22t=-22as,tbx-s,y-t其中，-3x3，-3y3，则c(-3,-3)=a(-2,-2)b(-1,-1)=3*(-1)=-3c(0,0)=a(-1,-1)b(1,1)+a(-1,0)b(1,0)+a(-1,1)b(1,-1)+a(0,-1)b(0,1)+a(0,0)b(0,0)+a(0,1)b(0,-1)+a(1,-1)b(-1,1)+a(1,0)b(-1,0)+a(1,1)b(-1,-1)=3*1+4*2+0*1

8、+2*0+1*0+3*0+1*(-1)+0*(-2)+2*(-1)=8c(3,3)=a(2,2)b(1,1)=4*1=4所以卷积结果为： -1 -3 -1 3 -2 0 4 -3 -6 -4 4 -4 2 11 -3 -7 -6 3 -6 4 15 -3 -11 -4 8 -10 3 17 -7 -11 2 5 -10 6 15 -8 -5 6 -4 -6 9 8 -3 -1 3 -3 -2 4 2【作业3】1、高斯型低通滤波器在频域中的传递函数是Hu,v=Ae-u2+v222根据二维傅里叶性质，证明空间域的相应滤波器形式为hx,y=A22e-222x2+y2这些闭合形式只适用于连续变量情况。

9、在证明中假设已经知道如下结论：函数e-x2+y2的傅立叶变换为e-u2+v2【解答】IDFTHu,v=-Ae-u2+v222ej2ux+vydudv=-Ae-u222+j2uxe-v222+j2vydudv=-Ae-122u2-j42uxe-122v2-j42vydudv=-Ae-122u2-j42ux-424x2+424x2e-122v2-j42vy-424y2+424y2dudv=-Ae-u-j22x222e-222x2e-v-j22y222e-222y2dudv令r=u-j22x，s=v-j22y，则du=dr，dv=ds，上式写成：=-Ae-r222e-222x2e-s222e-222

10、y2drds=Ae-222x2+y2-e-r222dr-e-s222ds=Ae-222x2+y22212-e-r222dr12-e-s222ds因为后两项是高斯分布，在-到积分为1，故上式等于：=A22e-222x2+y2=hx,y命题得证。2、第二版课本习题4.6（a）【解答】先来证明结论-1x+y=ejx+y根据欧拉公式展开等式右边ejx+y=cosx+y+jsinx+y因为x，y均为整数，故sinx+y=0，当x+y为奇数时，cosx+y=-1当x+y为偶数时，cosx+y=1故-1x+y=ejx+y再来证明题中等式DFTfx,y-1x+y=DFTfx,yejx+y=1MNx=0M-1y

11、=0N-1fx,yejx+ye-j2uxM+vyN=1MNx=0M-1y=0N-1fx,ye-j2-xM2M-yN2Ne-j2uxM+vyN=1MNx=0M-1y=0N-1fx,ye-j2xu-M2M+yv-N2N=Fu-M2,v-N23、观察如下所示图像。右边的图像这样得到：(a)用左侧图像乘以-1x+y；(b)计算离散傅里叶变换(DFT)；(c)对变换取复共轭；(d)计算离散傅里叶反变换；(e) 结果的实部再乘以-1x+y。用数学方法解释为什么会产生右图的效果。（忽略中间和右侧的黑白条纹，原题没有）【解答】已知IDFTFu,v=1MNu=0M-1v=0N-1Fu,vej2uxM+vyN=f

12、x,y则傅里叶变换的共轭复数进行傅里叶反变换的结果如下：IDFTF*u,v=1MNu=0M-1v=0N-1Fu,ve-j2uxM+vyN=1MNu=0M-1v=0N-1Fu,vej2u-xM+v-yN=f-x,-y设原始图像为fx,y，经过(a)变换后得到-1x+yfx,y经过(b)变换后得到Fu,v=1MNu=0M-1v=0N-1-1x+yfx,ye-j2uxM+vyN经过(c)变换后得到F*u,v=1MNu=0M-1v=0N-1-1x+yfx,yej2uxM+vyN经过(d)变换后得到IDFTF*u,v=1MNu=0M-1v=0N-11MNu=0M-1v=0N-1-1x+yfx,yej2u

13、xM+vyNej2uxM+vyN其实部为-1x+yf-x,-y，经过(e)变换后得到-1x+y-1x+yf-x,-y=f-x,-y最终效果是将原图像上下颠倒，左右颠倒，实现了旋转180度的效果。【作业4】1、请用公式列举并描述出你所知道的有关傅里叶变换的性质。【解答】1、时移性2、频移性3、均值4、共轭对称性5、对称性6、周期性7、线性8、微分特性9、卷积定理10、相关定理11、相似性12、几种特殊函数的傅里叶变换2、中文课本173页习题4.21【解答】没有区别。补0延拓的目的是在DFT相邻隐藏周期之间建立一个“缓冲区”。如果把左边的图像无限复制多次，以覆盖整个平面，那么将形成一个棋盘，棋盘中

14、的每个方格都是本图片和黑色的扩展部分。假如将右边的图片做同样的处理，所得结果也是一样的。因此，无论哪种形式的延拓，都能达到相同的分离图像的效果。3、（1）假设我们有一个0,1上的均匀分布随机数发生器U(0,1)，请基于它构造指数分布的随机数发生器，推导出随机数生成方程。（2）若我们有一个标准正态分布的随机数发生器N(0,1)，请推导出对数正态分布的随机数生成方程。【解答】（1）设U(0,1)可生成随机数w，用它来生成具有指数CDF的随机数z，其CDF具有下面的形式Fz=1-e-az0,z0,z0令F(z)=w，当z0时，解方程1-e-az=w得z=-1aln1-w由于w0,1，由上式可知z0，

15、因而不存在z0时，gs,tQ对gs,t有增强作用，由于“胡椒”噪声值较小，对加权平均结果的影响较小，所以滤波后噪声点处(x,y)取值和周围其他值更接近，有利于消除“胡椒”噪声。（b）当Q0时，gs,tQ对gs,t有削弱作用，由于“盐”噪声值较大，取倒数后较小，对加权平均结果的影响较小，所以滤波后噪声点处(x,y)取值和周围其他值更接近，有利于消除“盐”噪声。【作业5】1、请证明带通与带阻的频域关系公式，即课本中的关系公式Hbpu,v=1-Hbru,v【解答】证法一设D0是带宽的径向中心，W是带宽，D是D(u,v)距滤波器中心的距离。理想带通滤波器的频域公式为：Hbpu,v=1,D0-W2DD0

16、+W20,others理想带阻滤波器的频域公式为：Hbru,v=0,D0-W2DD0+W21,others由此可知Hbpu,v=1-Hbru,v图形化表示即为：图中黑色代表0，白色代表1，左图是带阻滤波器，右图是带通滤波器。可以发现两图是“互补”的，若将对应点数值相加，则全为1。证法二一张图像f (x, y)可拆分为带阻部分fr (x, y)和带通部分fp (x, y)，即fx,y=frx,y+fpx,y也可看做由带阻滤波器和带通滤波器分别卷积后叠加所得：fx,y=fx,y*hbrx,y+fx,y*hbpx,y等式两边取傅里叶变换得Fu,v=Fu,vHbru,v+Fu,vHbpu,vFu,v=

17、Fu,vHbru,v+Hbpu,vHbru,v+Hbpu,v=1Hbpu,v=1-Hbru,v2、复习理解课本中最佳陷波滤波器进行图像恢复的过程，请推导出w(x,y)最优解的计算过程，即从公式到的推导过程。【解答】因为2x,y=12a+12b+1s=-aat=-bbgx+s,y+t-wx,yx+s,y+t-gx,y-wx,yx,y2所以对w(x, y)求偏导数，并令其等于02x,ywx,y=wx,y12a+12b+1s=-aat=-bbgx+s,y+t-gx,y+wx,yx,y-x+s,y+t2=0令a=x,y-x+s,y+t，b=gx+s,y+t-gx,y，N=2a+12b+1，则上式可简写

18、成2x,ywx,y=wx,y1Ns=-aat=-bbawx,y+b2=1Ns=-aat=-bbawx,y+b2wx,y=0已知awx,y+b2wx,y=2awx,y+ba=0的解为wx,y=-ba所以wx,y=1Ns=-aat=-bbgx+s,y+t-gx,yx+s,y+t-x,y=1Ns=-aat=-bbgx+s,y+t-gx,yx+s,y+t+x,yx+s,y+t-x,yx+s,y+t+x,y令g=gx+s,y+t，=x+s,y+t，则上式等于=1Ns=-aat=-bbg-gx,y+gx,y-gx,yx,y2x+s,y+t-2x,y=gx,yx,y-gx,yx,y+gx,yx,y-gx,y

19、x,y2x,y-2x,y=gx,yx,y-gx,yx,y2x,y-2x,y推导成立。3、考虑在x方向均匀加速导致的图像模糊问题。如果图像在t = 0静止，并用均匀加速x0(t) = at2/2加速，对于时间T, 找出模糊函数H(u, v)，可以假设快门开关时间忽略不计。【解答】由定义可知Hu,v=0Te-j2ux0t+vy0tdt将x0t=at22，y0t=0代入得Hu,v=0Te-j2uat22dt=0Te-juat2dt上述积分难以化简，故不再推导。4、已知一个退化系统的退化函数H(u, v)，以及噪声的均值与方差，请描述如何利用约束最小二乘方算法计算出原图像的估计。【解答】频域中原图像的

20、估计由下式给出Fu,v=H*u,vHu,v2+Pu,v2Gu,v其中P(u, v)是拉普拉斯算子的傅里叶变换。定义“残差”向量r=g-Hf，由于Fu,v是的函数，则f和r都是的函数。令=rTr=r2，则它是的单调递增函数。再调整使r2=2a，a是一个精确度因子。已知噪声的均值为mu，方差2，和Hu,v、Pu,v，（1）设定一个的初始值，（2）计算r2，（3）若满足r2=2a则执行第4步，若不满足，则调整大小，然后返回第2步。（4）使用最新的，计算Fu,v=H*u,vHu,v2+Pu,v2Gu,v（5）再通过傅里叶反变换即可得到估计图像。【作业6】1、r, g, b是RGB彩色空间沿R,G,B轴

21、的单位向量，定义向量u=Rxr+Gxg+Bxbv=Ryr+Gyg+Byb将gxx，gyy，gxy定义为这两个向量的点乘：gxx=uu=uTu=Rx2+Gx2+Bx2gyy=vv=vTv=Ry2+Gy2+By2gxy=uv=uTv=RxRy+GxGy+BxBy推导出最大变换率方向和点(x, y)在方向上变化率的值F()。【解答】要求最大变换率方向，即求使ucos+vsin2取最大值的。ucos+vsin2=u2cos2+2uvsincos+v2sin2=12gxx1+cos2+gxysin2+12gyy1-cos2=12gxx+gyy+12gxx-gyycos2+gxysin2将上式对求偏导得1

22、2gxx+gyy+12gxx-gyycos2+gxysin2=-gxx-gyysin2+2gxycos2令上式等于0，解方程求出的值-gxx-gyysin2+2gxycos2=0gxx-gyysin2=2gxycos2tan2=2gxygxx-gyy=12tan-12gxygxx-gyy此即是最大变换率方向，对应的变化率值为F=12gxx+gyy+12gxx-gyycos2+gxysin2【作业7】1、请根据课本中Z变换的定义，证明如下结论。（1）若xn的Z变换为Xz，则-1nxn的Z变换为X-z。（2）若xn的Z变换为Xz，则x-n的Z变换为X(1/z)。【解答】（1）由Z变换定义可知，xn

23、的Z变换是：Xz=-xnz-n则-1nxn的Z变换为：-1nxnz-n=-1-nxnz-n=-xn-z-n=X-z（2）x-n的Z变换为：-x-nz-n令n=-n，则上式可写成：-xnzn=-xn1z-n=X1z2、若G1z=-z-2K+1G0-z-1成立，请证明g1n=-1n+1g02K-1-n。【解答】已知g1n的Z变换是G1z，g0n的Z变换是G0z，根据Z变换的时间翻转性x-nXz-1可得g0-nG0z-1根据Z变换的平移性xn-kz-kXz可得g0-n-2K+1z-2K+1G0z-1根据Z变换的Z域翻转性-1nxnX-z可得-12K-1-ng02K-1-n-z-2K+1G0-z-1-

24、12K-1-1-ng02K-1-n-1-2K+1z-2K+1G0-1-1z-1因为2K-1为奇数，所以-12K-1=-1，-1-n=-1n，-2K+1为奇数，故-1-2K+1=-1，上式继续推导得-1n+1g02K-1-n-z-2K+1G0-z-1=G1z所以G1z的Z反变换是-1n+1g02K-1-n又因为g1n也是G1z的Z反变换，所以g1n=-1n+1g02K-1-n3、假设课本中给出完美重建滤波器的正交族对应的三个滤波器间的关系式是正确的，请以此为基础，推导h0，h1的关系。【解答】已知g0n=-1n+1h1ng1n=-1nh0ng1n=-1n+1g02K-1-n将第1、2个式子代入第

25、3个式子得到-1nh0n=-1n+1-12K-1-n+1h12K-1-nh0n=-12K-n+1h12K-1-nh0n=-1n+1h12K-1-n若g0n=-1nh1ng1n=-1n+1h0ng1n=-1n+1g02K-1-n则将第1、2个式子代入第3个式子得到-1n+1h0n=-1n+1-12K-1-nh12K-1-nh0n=-12K-1-nh12K-1-nh0n=-1n+1h12K-1-n因此，两种情况的解一致，都是：h0n=-1n+1h12K-1-n4、请证明完美重建滤波器组的双正交性质，即课本282页的公式7.1.20。【解答】已知书上公式7.1.9，7.1.10，7.1.18H0zG

26、0z+H1zG1z=2H0-zG0z+H1-zG1z=0H0zG0z+H0-zG0-z=2调制矩阵Hmz=H0zH0-zH1zH1-z用调制矩阵的行列式表示G0z和G1zG0z=2detHmzH1-zG1z=-2detHmzH0-z（1）先来证第1个公式：g1k,h12n-k=n令Pz=G1zH1z，将式代入得Pz=G1zH1z=-2detHmzH0-zH1z由于detHm-z=-detHmz，则将式两端同时乘以H0z得G0zH0z=2detHmzH0zH1-z=P-z因此G0zH0z=G1-zH1-z，代入式得G1zH1z+G1-zH1-z=2根据Z变换的Z域翻转性-1nxnX-z，对上式反

27、Z变换得到kg1kh1n-k+-1nkg1kh1n-k=2n由于冲激函数n在n=0时等于1，其他情况等于0，且n为奇数时，上式奇次方项可以相互抵消，因此只取n为偶数的情况，用2n代替上式的n可得2kg1kh12n-k=22n=2nkg1kh12n-k=g1k,h12n-k=n（2）再来证第2个公式：g0k,h12n-k=0由式可得H0-z=-H1-zG1zG0z得H0-zG0-z-H1zG1z=0将上面两式合并得-H1-zG1zG0zG0-z-H1zG1z=0H1-zG1zG0-z+H1zG1zG0z=0由于G1z是滤波器，不为0，故可将上式两端同时除以G1z得H1-zG0-z+H1zG0z=

28、0根据Z变换的Z域翻转性-1nxnX-z，对上式反Z变换得到kg0kh1n-k+-1nkg0kh1n-k=0由于n为奇数时，上式奇次方项可以相互抵消，因此只取n为偶数的情况，用2n代替上式的n可得2kg0kh12n-k=0kg0kh12n-k=g0k,h12n-k=0（3）最后证第3个公式：g1k,h02n-k=0由式可得H1-z=-H0-zG0zG1z得H1-zG1-z-H0zG0z=0将上面两式合并得-H0-zG0zG1zG1-z-H0zG0z=0H0-zG0zG1-z+H0zG0zG1z=0由于G0z是滤波器，不为0，故可将上式两端同时除以G0z得H0-zG1-z+H0zG1z=0根据Z

29、变换的Z域翻转性-1nxnX-z，对上式反Z变换得到kg1kh0n-k+-1nkg1kh0n-k=0由于n为奇数时，上式奇次方项可以相互抵消，因此只取n为偶数的情况，用2n代替上式的n可得2kg1kh02n-k=0kg1kh02n-k=g1k,h02n-k=05、请围绕本周课堂讲授的内容编写至少一道习题，并给出自己的分析解答。题目形式可以是填空题、选择题、判断对错题、计算题、证明题。发挥你的创造力吧。【题目】设计一个系统，用下图的编码器对预测残差金字塔进行编码，并画出框图。【解答】对应的解码器如下图所示。【作业8】1、哈尔变换可以用矩阵的形式表示为：T=HFHT其中，F是一个NN的图像矩阵，H

30、是NN变换矩阵，T是NN变换结果。对于哈尔变换，变换矩阵H包含基函数hkz，它们定义在连续闭区间z0,1，k=0,1,2N-1，其中N=2n。为了生成H矩阵，定义整数k，即k=2p+q-1（这里0pn-1，当p=0时q=0或1；当p0时，1q2p）。可得哈尔基函数为：h0z=h00z=1N,z0,1hkz=hpqz=1N2p2,q-12pzq-0.52p-2p2,q-0.52pzq2p0,其它，z0,1NN哈尔变换矩阵的第i行包含了元素hi(z)，其中z=0N,1N,N-1N。计算当N=16时的H16矩阵。【解答】由N=16可知，k=015，根据公式k=2p+q-1计算p、q如下表所示：k01

31、23456789101112131415p0011222233333333q0112123412345678计算hkz=hpqz，得到H矩阵：2、课本322页习题7.10的（a）与（b）小题。【解答】（a）由于展开函数0和1构成正交基：0,1=1212-1212=00,0=1212+1212=11,1=1212+-12-12=1所以展开系数0=0T,f=121232=5221=1T,f=12-1232=22可以保证00+11=5221212+2212-12=32=f（b）由于展开函数0和1双正交：0,1=10+01=01,0=11+1-1=00,0=11+0-1=11,1=10+11=1所以展

32、开系数0=0T,f=1-132=11=1T,f=0132=2可以保证00+11=10+211=32=f3、课本323页习题7.11【解答】令0,0x=x，则1,0x=22x，1,1x=22x-1，三个函数的图像如下所示由上图可知，0,0x无法用1,0x和1,1x的线性加权和表示出来，因此本题给定的尺度函数x不满足多分辨率分析的第2个要求。4、课本323页习题7.16【解答】（a）因为本题是单尺度变换，开始尺度j0=1，所以j只能是1，相应的k=0或1，根据书上公式（7.3.5）和（7.3.6）计算M=4的一维DWT系数。W1,0=1Mn=03fn1,0n=1212+42-30+00=522W1

33、,1=1Mn=03fn1,1n=1210+40-32+02=-322W1,0=1Mn=03fn1,0n=1212-42-30+00=-322W1,1=1Mn=03fn1,1n=1210+40-32-02=-322所以DWT系数为522,-322,-322,-322，函数fn的展开形式为fn=2451,0n-31,1n-31,0n-31,1n（b）根据上式结果f1=2452-30-3-2-30=45、请围绕本周课堂讲授的内容编写至少一道习题，并给出自己的分析解答。题目形式可以是填空题、选择题、判断对错题、计算题、证明题。发挥你的创造力吧。【题目】计算图像F=3-162的哈尔变换。【解答】根据公式

34、T=HFHT，取变换矩阵H为H2=12111-1H2T=12111-1T=H2FH2T=12111-13-162111-1=1291-3-3111-1=12108-60=54-30【作业9】1、课本323页习题7.21【解答】（a）根据书上公式7.3.5和7.3.6Wj0,k=1Mxfxj0,kxWj,k=1Mxfxj,kx可得当尺度J=3，j0=0，M=8，f(n)=1(n=0,1,.7)时W0,0=18nfn0,0n=122118=22W0,0=18nfn0,0n=122114+1-14=0W1,0=18nfn1,0n=122122+1-22+104=0W1,1=18nfn1,1n=122

35、104+122+1-22=0W2,0=18nfn2,0n=1222-2+06=0W2,1=18nfn2,1n=12202+2-2+04=0W2,2=18nfn2,2n=12204+2-2+02=0W2,3=18nfn2,3n=12206+2-2=0所以变换系数为：W0,0,W0,0,W1,0,W1,1,W2,0,W2,1,W2,2,W2,3=22,0,0,0,0,0,0,0（b）当输入变为f(n)=1,1,1,1,-1,-1,-1,-1,(n=0,1,.7)时，上面8个公式可以算得：W0,0=18nfn0,0n=122114+1-14=0W0,0=18nfn0,0n=122114+-1-14=

36、22W1,0=18nfn1,0n=122122+1-22-104=0W1,1=18nfn1,1n=122104-122-1-22=0W2,0=18nfn2,0n=1222-2+06=0W2,1=18nfn2,1n=12202+2-2+04=0W2,2=18nfn2,2n=12204-2+2+02=0W2,3=18nfn2,3n=12206-2+2=0所以变换系数为：W0,0,W0,0,W1,0,W1,1,W2,0,W2,1,W2,2,W2,3=0,22,0,0,0,0,0,0（c）因为W2,2=18nfn2,2n=12204+2f4+-2f5+02=B所以设输入序列f(n)=0,0,0,0,x

37、,-x,0,0,(n=0,1,.7)，则12204+2x+-2-x+02=2x=Bx=22B所以输入序列fn=0,0,0,0,22B,-22B,0,02、课本325页习题7.26【解答】（a）（b）3、请围绕本周课堂讲授的内容编写至少一道习题，并给出自己的分析解答。题目形式可以是填空题、选择题、判断对错题、计算题、证明题。发挥你的创造力吧。【题目】对于函数f(n)=1,4,-3,0，当起始尺度j0=1时，在区间0,3上的一维DWT系数如下：W1,0,W1,1,W1,0,W1,1=522,-322,-322,-322画出上述变换所需的FWT滤波器组，标记所有的输入和输出。【解答】【作业10】1、

38、信息论的相关概念也是进行图像分析的一个重要的工具，请形式化简要叙述信息论中的如下五个重要概念：（1）信息量；（2）信息熵；（3）条件熵；（4）互信息；（5）信道的容量。【解答】（1）信息量IE=log1PE=-logPE（2）信息熵Hz=-j=1JPajlogPaj（3）条件熵Hz|bk=-j=1JPaj|bklogPaj|bk（4）互信息Iz,v=Hz-Hzv=j=1Jk=1KP(aj,bk)logPaj,bkPajPbk（5）信道的容量C=maxzIz,v2、对于一个二值的信息源，可以用一个参量Pbs来表示二值信源字幕的概率分布，假设二值信息在一个二值对称信道上传送，该信道由单一错误率Pe

39、来刻画，请证明信道传递的信息，即互信息可由如下公式来计算：Iz,v=Hbspbspe+pbspe-Hbspe请参考课本341页。【解答】为书写方便，我们约定下面的证明过程用a代替pbs，用b代替pe，用a代替pbs，用b代替pe。根据书上公式8.3.12可知Iz,v=j=1Jk=1KP(aj)qkjlogqkji=1JPaiqki本题中Pa1=a,Pa2=a,q11=b,q12=b,q21=b,q22=b，代入上式得Iz,v=ablogbab+ab+ablogbab+ab+ablogbab+ab+ablogbab+ab=ablogb-ablogab+ab+ab logb-ab logab+ab

40、+ab logb-ab logab+ab+ablogb-ablogab+ab=ab+ablogb+ab+ablogb-ab+ablogab+ab-ab+ablogab+ab=blogb+blogb-ab+ablogab+ab-ab+ablogab+ab根据公式Hbst=-tlog2t-tlog2t可知要证等式的右端等于Hbsab+ab-Hbsb=-ab+ablogab+ab-1-ab-ablog1-ab-ab+blogb+blogb其中1-ab-ab=1-ab-1-a1-b=-ab+a+b-ab=a1-b+1-ab=ab+ab所以代入前式可得Hbsab+ab-Hbsb=-ab+ablogab+ab-ab+ablogab+ab+blogb+blogb=Iz,v等式Iz,v=Hbspbspe+pbspe-Hbspe成立。29