《吉林大学22春《过程控制与自动化仪表》综合作业一答案参考47》由会员分享，可在线阅读，更多相关《吉林大学22春《过程控制与自动化仪表》综合作业一答案参考47（17页珍藏版）》请在金锄头文库上搜索。
1、吉林大学22春过程控制与自动化仪表综合作业一答案参考1. 对于给定的任意符号序列，给出其算术编码和解码的程序，并给出编解码示例。对于给定的任意符号序列，给出其算术编码和解码的程序，并给出编解码示例。正确答案：算术编码步骤如下：rn (1)统计信息流中每个符号的概率。rn (2)设置“当前区间”为01。rn (3)对输出字符串或信息流的每个符号编码器按步骤和进行处理：rn 将“当前区间”分成子区间该子区间的长度正比于符号的概率。rn 选择子区间对应于下一个信源符号并使它成为新的“当前区间”。rn (4)将整个输入信息流处理后输出的“当前区间”中任意取一个数就是该给输入信息流的算术编码。rn 以信
2、息流33214为例首先统计各符号的概率：P1=02P2=02rnP3=04P4=02。图62是该信息流的算术编码过程图解。rnrnMatlab源程序如下：rn(1)调用编码、解码函数进行算术编码和解码的示例程序(函数)。rn 调用编码、解码函数进行算术编码和解码的示例程序(函数)rnfunctionoutinv_x=test_arth(in)rnclc；rnclear all；rnclose all；rnformat longrn in：输入信息流rn out：输出小数rn inv_X解码输出信息流rnin=3 3 2 1 4； 举例rntempP_spacesymbolL=en_code(i
3、n)； 编码rninv_x=de_code(tempP_spacesymbolL) 解码rn(2)算术编码程序(函数)。rn 编码函数rnfunctiontempp_spacesymbolL=en_code(in)rnsymbol=unique(in)；rnL=length(in)；rnn=length(symb01)；rninrn 符号的概率统计rnfor i=1：nrn p(i)=length(find(in=symbol(i)L；rnendrn 建立符号的概率间隔空间rntemp=0；rnfor i=1：nrn p_space(i1：2)=temptemp+p(i)；rn temp=te
4、mp+p(i)；rnendrn 编码rnout=P_space(sym_ind(symbolin(1)1：2)；rnfor i=2：Lrn out=out(1)+(out(2)一out(1)*pl-space(sym_ind(symbolin(i)1：2)；rnendrn 输出间隔区间rn outrn 寻找最后区间内的一个二进制小数temprn temp1=short_ab(out(1)out(2)；rn temprn 编码二进小数temp需要的比特数rn 1rn function i=sym_ind(symbola)rn 符号到序号的索引rn i=find(symbol=a)；rn 寻找区间
5、ab中二进制表示最少的数rn functionc1=short_ab(ab)rn C=0；rn i=1；rn while(1)rn a=2*a：rn b=2*b：rn if floor(a)=floor(b)rn break；rn endrn C=C+floor(a)*2一i：rn i=i+1：rn a=afloor(a)；rn b=bfloor(b)；rn endrn C=C+2一(i)；rn l=i；rn (3)算术解码程序(函数)。rn 解码函数rn functioninv_x=de_code(tempP_spacesymbolL)rn 解码rn for i=1：Lrnk=spac_in
6、d(p_spacetemp)；rn inv_x(i)=symbol(k)；rn range=P_space(k2)一P_space(k1)；rn temp=(tempP_space(k1)range；rn endrn 解码输出rn function b=spac_ind(spacela)rn 间隔区间到序号的索引rn M=size(spacel1)；rn for i=1：Mrn if a=spacel(i1)&a算术编码步骤如下：(1)统计信息流中每个符号的概率。(2)设置“当前区间”为0,1。(3)对输出字符串或信息流的每个符号,编码器按步骤和进行处理：将“当前区间”分成子区间,该子区间的长
7、度正比于符号的概率。选择子区间对应于下一个信源符号,并使它成为新的“当前区间”。(4)将整个输入信息流处理后,输出的“当前区间”中任意取一个数就是该给输入信息流的算术编码。以信息流3,3,2,1,4为例,首先统计各符号的概率：P1=02,P2=02,P3=04,P4=02。图62是该信息流的算术编码过程图解。Matlab源程序如下：(1)调用编码、解码函数进行算术编码和解码的示例程序(函数)。调用编码、解码函数进行算术编码和解码的示例程序(函数)functionout,inv_x=test_arth(in)clc；clearall；closeall；formatlongin：输入信息流out：
8、输出小数inv_X解码输出信息流in=33214；举例temp,P_space,symbol,L=en_code(in)；编码inv_x=de_code(temp,P_space,symbol,L)解码(2)算术编码程序(函数)。编码函数functiontemp,p_space,symbol,L=en_code(in)symbol=unique(in)；L=length(in)；n=length(symb01)；in符号的概率统计fori=1：np(i)=length(find(in=symbol(i)L；end建立符号的概率间隔空间temp=0；fori=1：np_space(i,1：2)=
9、temp,temp+p(i)；temp=temp+p(i)；end编码out=P_space(sym_ind(symbol,in(1),1：2)；fori=2：Lout=out(1)+(out(2)一out(1)*pl-space(sym_ind(symbol,in(i),1：2)；end输出间隔区间out寻找最后区间内的一个二进制小数temptemp,1=short_ab(out(1),out(2)；temp编码二进小数temp需要的比特数1functioni=sym_ind(symbol,a)符号到序号的索引i=find(symbol=a)；寻找区间a,b中二进制表示最少的数functio
10、nc,1=short_ab(a,b)C=0；i=1；while(1)a=2*a：b=2*b：iffloor(a)=floor(b)break；endC=C+floor(a)*2一i：i=i+1：a=afloor(a)；b=bfloor(b)；endC=C+2一(i)；l=i；(3)算术解码程序(函数)。解码函数functioninv_x=de_code(temp,P_space,symbol,L)解码fori=1：Lk=spac_ind(p_space,temp)；inv_x(i)=symbol(k)；range=P_space(k,2)一P_space(k,1)；temp=(tempP_sp
11、ace(k,1)range；end解码输出functionb=spac_ind(spacel,a)间隔区间到序号的索引M=size(spacel,1)；fori=1：Mifa=spacel(i,1)&a2. 自动化系统中，进行参数检测，传送到控制室，进行显示、记录的是：( )A.自动检测系统B.自动控制系统C.自动操作系统D.自动联锁保护系统参考答案：A3. 架空光缆接续时，接头盒内光缆的加强芯应( )A.加强芯在接头盒内要余留B.加强芯要分别固定在两架空光缆接续时，接头盒内光缆的加强芯应( )A.加强芯在接头盒内要余留B.加强芯要分别固定在两边的紧固件上，固定要牢靠C.对加强芯无要求D.两端
12、的加强芯要确保联通答案：B4. 金属材料在载荷作用下，断裂前发生塑性变形(永久变形)而不被破坏的能力称为塑性。( )金属材料在载荷作用下，断裂前发生塑性变形(永久变形)而不被破坏的能力称为塑性。( )答案：对5. 编写一个RGB与HSI相互转换的计算程序，并对一实际真彩图像进行变换。编写一个RGB与HSI相互转换的计算程序，并对一实际真彩图像进行变换。正确答案：RGB与HSI的相互转换的Matlab程序段如下其中程序段(1)将文件IMGRGBDAT中的R、G、B分量转换为H、S、I分量；程序段(2)将文件IMGHSLDAT中的H、S、I分量转换为R、G、B分量。rn (1)RGB转换为HSI如
13、下：rn function hsi=rgb2hsi(image)rn image=imread(IMGRGBDAT)；rn image=im2double(image)；rn r=image(：1)；rn g=xmagel(：z)；rn b=image(：3)；rn num=05*(rg)+(rb)；rn den=sqrt(rg)2+(rb)*(gb)；rn theta=acos(num(den+eps)；rn H=theta；rn H(bg)=2*piH(bg)；rn H=H(2*pi)；rn num=min(min(rg)b)；rn den=r+g+b：rn den(den= =0)=eps；rn S=13*numden：rn H(S=0)=0；rn I=(r+g+b)3；rn hsi=cat(3HSI)；rn figure(1)；rn imshow(hsi)；rn (2)HSI转换为RGB如下：rn function rgb=hsi2rgb(image)rn image=imread(IMGHSIDAT)；rn image=im2double(image)；rnH=image(：1)*2*pi；rn S=image(：2)；rn I=image(：3)；rn R=zeros(size(image1)size(image2)；rn G=zero
