matlab常见的错误提示及其解决方法 (2011-07-23 09:08:09) 分类：matlab学习 标签： matlab 常见错误 原因 解决方法 it 今天在百度知道上收到求助，说想知道matlab常见错误。个人虽然也学了点matlab皮毛，并且碰过不少问题，但从没有认真总结过常见的错误。于是上网搜索一番，整理出常见的错误，当前以整理出的内容如下，不足之处待以后补充，也期待看到此篇博文的朋友帮忙指出不足。 1.Undefined function or variable "U" 中文解释：函数或变量 U 没有定义 出错原因：函数或变量 U 没有定义 解决办法：仔细检查错误所指向的位置，寻找没有定义的函数或者变量，加以定义。 2.Unexpected MATLAB expression. 中文解释：不规范的matlab表达式 出错原因：少了运算符号，或者括号等情况 解决办法：仔细检查错误所指向的位置，检查是否少了乘号，或者逗号、分号、括号等。 3.Subscript indices must either be real positive integers or logicals
中文解释：下标索引必须是正整数类型或者逻辑类型 出错原因：在访问矩阵（包括向量、二维矩阵、多维数组，下同）的过程中，下标索引要么从 0 开始，要么出 现了负数。注：matlab 的语法规定矩阵的索引从 1 开始，这与 C 等编程语言的习惯不一样。 解决办法：自己调试一下程序，把下标为 0 或者负数的地方修正。 4.Matrix dimensions must agree 中文解释：矩阵的维数必须一致 出错原因：这是由于运算符(= + - / * 等)两边的运算对象维数不匹配造成的，典型的出错原因是错用了矩阵运 算符。matlab 通过“.”来区分矩阵运算和元素运算。 解决办法：自己调试一下程序，保证运算符两边的运算对象维数一致。 5.Function definitions are not permitted at the prompt or in scripts 中文解释：不能在命令窗口或者脚本文件中定义函数 出错原因：一旦在命令窗口写 function c = myPlus(a,b)，此错误就会出现，因为函数只能定义在 m 文件中。 关于脚本文件和 m 文件的区别可查阅 matlab 基础书。 简言之： 1) 如果你写成 function 的形式，那么必须写在 m 文件中，且以function 开头（即 function 语句前不能包含其他语句，所有语句必须放在function 中，当然，function 的定义可以有多个，各 function 之间是并列关系，不能嵌套）； 2) 如果你写成脚本的形式，则既可以写在命令窗口中，也可以写在 m 文件中，但两者均不能包含 function 语句（即不能进行函数的定义） 解决办法：新建一个 m 文件，然后再进行函数的定义
滨江学院 《信息论与编码》课程论文题目香农编码及其应用改善 院系电子工程系 专业班级通信班 学生姓名 学号 教师杨玲 成绩 二Ｏ一四年十二月二十二日
香农编码及其应用改善 摘要：香农编码作为变长信源编码的重要方法之一，具有重要的理论指导意义，但其在实际应用中存在效率较低的缺点。本文对香农编码方法进行阐述，及运用MATLAB实现香农编码操作，并找出香农编码的不足，针对其缺陷，通过判断码字之间是否互为前缀来确定码字的方法对其编码算法进行了优化，给出了优化算法的实现步骤。最后，通过具体实例分析得出本文提出的改善算法能有效地提高编码效率。 关键词：香农码方法；MATLAB；编码效率；优化编码；
引言：1948年，美国工程师香农在贝尔实验室杂志上发表了长文《通讯的数学原理》他用概率测度和数理统计的方法系统地讨论了通信的基本问题，得出了几个重要而带有普遍意义的结论，并由此奠定了现代信息论的基础。香农编码理论揭示了在通信系统中，采用适当的编码后能够实现高效率和高可靠地传输信息的规律，并给出了相应的信源编码定理和信道编码定理。从数学观点看，这些定理是最优编码的存在定理。它们给出了编码的性能极限，在理论上阐明了通信系统中各种因素的相互关系，为寻找最佳通信系统提供了重要的理论依据。 在多媒体数据的传输和存储过程中，为了确保通信的顺利进行，必须要通过信源编码技术必须对多媒体信息进行压缩处理。香农编码技术作为变长信源编码的重要方法之一，具有重要的理论指导意义。但由于在香农编码的过程中先限定每个码字的码长，以至于在码字的选取中是以每个码字的码长作为先决条件而不考虑各个码字之间的相关性，因此编出的码字往往存在较大的冗余，影响了整个通信系统的传输效率。就这一缺陷，本文提出了通过剔除先限定每个码字的码长这一过程，通过判断码字之间是否互为前缀来确定码字的方法对其编码算法进行了改善。 1.香农编码的方法 在写香农编码之前先简单介绍下信源编码： 编码分为信源编码和信道编码，其中信源编码又分为无失真和限失真。由于这些定理都要求符号数很大，以便其值接近所规定的值，因而这些定理被称为极限定理。一般称无失真信源编码定理为第一极限定理；信道编码定理（包括离散和连续信道）称为第二极限定理；限失真信源编码定理称为第三极限定理。完善这些定理是香农信息论的主要内容。 信源编码的基础是信息论中的两个编码定理：无失真编码定理和限失真编码地宫里，前者是可逆编码的基础。可逆是指当信源符号转换成代码后，可从代码无失真的恢复原信源符号。当已知信源符号的概率特性时，可计算它的符号熵，这表示每个信源符号所载有的信息量。编码定理不但证明了必定存在一种编码方法，可使代码的平均长度可任意接近但不低于符号熵，而且还阐明达到这目标的途径，就是使概率与码长匹配。无失真编码或可逆编码只适用于离散信源。本节讨论离散信源编码。首先从无失真编码定理出发，重点讨论以香农码为代表的最
matlab常用操作备忘（1）2007-11-30 22:01:06 分类： 北京理工大学 20981 陈罡 帮助朋友做几个数据的卷积的仿真，一用才知道，呵呵，发现对不住偶的导师了。。。好多matlab的关键字和指令都忘记了。特意收集回顾一下： （1）管理命令和函数 addpath :添加目录到MATLAB搜索路径 doc :在Web浏览器上现实HTML文档 help :显示Matlab命令和M文件的在线帮助 helpwin helpdesk :help 兄弟几个 lookfor :在基于Matlab搜索路径的所有M文件中搜索关键字 partialpath:部分路径名 8*) path :所有关于路径名的处理 pathtool :一个不错的窗口路径处理界面 rmpath :删除搜索路径中指定目录 type :显示指定文件的内容 ver :版本信息 version :版本号 web :打开web页 what :列出当前目录吓所有的M文件 Mat文件和 Mex文件 whatsnew :显示readme文件 which :显示文件位置 （2）管理变量和工作区 clear :从内存中删除所有变量，clear x y z是删除某个变量 disp :显示文本或数组内容 length :数组长度(最长维数) load :重新载入变量(从磁盘上) mlock :锁定文件，防止文件被错误删除 munlock :解锁文件 openvar :在数组编辑器中打开变量 pack :整理内存空间 save :保存变量到文件 8*) size :数组维数
who whos :列出内存变量 workspace :显示工作空间窗口 （3）管理命令控制窗口(command窗口) clc :清空命令窗口 echo :禁止或允许显示执行过程 format :设置输出显示格式 home :光标移动到命令窗口左上角 more :设置命令窗口页输出格式 （4）文件和工作环境 cd :改变工作目录 copyfile :复制文件 delete :删除文件和图形对象 diary :把命令窗口的人机交互保存到文件 dir :显示目录 edit :编辑文本文件 fileparts :返回文件的各个部分 fullfile :使用指定部分建立文件全名 inmem :返回内存(伪代码区)的matlab函数名 ls :在unix系统中列出目录(win中亦可) matlabroot :根目录 mkdir :新建目录 open :打开文件 pwd :显示当前目录 tempdir :返回系统临时目录的名字 tempname :随机给出一个临时字符串(可用作文件名) ! :直接调用操作系统command命令 （5）启动和推出matlab matlabrc :Matlab的启动M文件 exit quit :退出Matlab startup :运行matlab启动文件 （6）程序设计 builtin :从可重载方法中调用内置函数 eval :执行包含可执行表达式的字符串
Matlab期中作业 标题：基于matlab的信息隐藏代码实现姓名：******** 班级：******** 学号：********
目录 摘要----------------------------------------------------------------------------------------------------------------------------- - 2 -前言----------------------------------------------------------------------------------------------------------------------------- - 2 - 一、图像隐藏技术的研究现状 ------------------------------------------------------------------------------------ - 2 - 二、实现图像隐藏的主要算法 ------------------------------------------------------------------------------------ - 4 - 三、基于行列变换的图像置乱与反置乱算法 ---------------------------------------------------------------- - 4 - 3.1 图像置乱算法 ----------------------------------------------------------------------------------------------------- - 4 - 3.2 图像反置乱算法-------------------------------------------------------------------------------------------------- - 5 - 四、基于m序列的图像加密与解密算法----------------------------------------------------------------------- - 5 - 五、相关程序代码及运行结果 ------------------------------------------------------------------------------------ - 5 - 5.1 源代码 --------------------------------------------------------------------------------------------------------------- - 5 - 5.2 运行结果 ------------------------------------------------------------------------------------------------------------ - 7 - 六、参考文献------------------------------------------------------------------------------------------------------------ - 8 - 七、总结 ------------------------------------------------------------------------------------------------------------------ - 8 -
《信息论与编码》教学大纲 一课程简介 课程编号：04254002 课程名称：信息论与编码Informatics & Coding 课程类型：基础课必修课 学时：32 学分：2 开课学期：第六学期 开课对象：通信、电子专业 先修课程：概率论与数理统计、信号与系统、随机信号原理。 参考教材：信息论与编码，陈运，周亮，陈新，电子工业出版社，2002年8月 二课程性质、目的与任务 信息论在理论上指出了建立最佳编码、最佳调制和最佳接收方法的最佳系统的理论原则，它对通信体制和通信系统的研究具有指导意义。提高信息传输的可靠性和有效性始终是通信工作所追求的目标。因此，信息论与编码是从事通信、电子系统工程的有关工程技术人员都必须掌握的基本理论知识。 内容提要：本课程包括狭义相对论和提高通信可靠性的差错控制编码理论。信息论所研究的主要问题是在通信系统设计中如何实现有效性和可靠性。 三教学基本内容与基本要求 本课程总学时为32。其中理论教学为28，实验学时为4。 主要的理论教学内容包括：离散信源和连续信源的熵、条件熵、联合熵和平均互信息量的概念及性质；峰值功率受限和平均功率受限下的最大熵定理和连续信源熵的变换；变长码的霍夫曼编码方法，熟悉编码效率和平均码长的计算；最大后验概率准则和最大似然译码准则等。 实验内容主要包括：离散无记忆信道容量的迭代算法，循环码的编译码。 四教学内容与学时分配 第3章离散信源无失真编码
第6章网络信息论 （教学要求：A—熟练掌握；B—掌握；C—了解） 五实习、实验项目及学时分配 1．离散无记忆信道容量的迭代算法2学时 要求用Matlab编写计算离散信道容量的实用程序并调试成功，加深对信道容量的理解。 2．循环码的编译码2学时 要求用Matlab编写程序，用软件完成循环码的编译码算法。 六教学方法与手段 常规教学与多媒体教学相结合。
matlab常见错误命令 >> a=【0 1 2 3 4 5 6 7 8 9】 ??? a=【0 1 2 3 4 5 6 7 8 9】 Error: The input character is not valid in MATLAB statements or expressions. error ['er?]n. 错误,过失input ['input]n. 输入character ['k?rikt?]n. 个性,字符,人物valid ['v?lid]a。有确实根据的,有效的,正当的statement ['steitm?nt]n. 声明,陈述expression [iks'pre??n]n. 表达,表示,表现,表情,措辞,词句expressions n.表示式公式>> x=-5:5;y=-5:5;z=x^2+y^2;polt(x,y,z) Error using ==> mpower Matrix must be square. matrix ['meitriks]n. 母体,子宫，细胞，脉石,字母n.矩阵 square [skw??]n. 正方形,街区,平方a. 正方形的,正直的,公正的v. 一致,符合,使...成方形>> x=-5:5;y=-5:5;z=x.^2+y.^2;polt(x,y,z) Undefined function or method 'polt' for input arguments of type 'double'. undefined ['?ndi'faind]未定义的function ['f??k??n]n.功能,函数,重大聚会vi. 运行,起作用method ['meθ?d]n. 方法,办法argument ['ɑ:gjum?nt]n. 辩论,争论,论据,理由
文章编号:1006-2475(2005)11-0028-04 收稿日期:2004-12-20作者简介:叶天语(1982-),男,浙江温州人,汕头大学电子工程系硕士研究生,研究方向:数字图像处理与信息安全。 基于图像LSB 的秘密信息隐藏技术 叶天语 (汕头大学电子工程系,广东汕头 515063) 摘要:基于图像的信息隐藏技术有着广泛的应用,而最低有效位算法是其中有效的方法之一。本文首先对载体图像进行位平面分析,发现LSB 位具有最好的隐藏效果;然后分别实现了三个以二值图像和文本文件作为秘密信息的信息隐藏和提取实验。实验结果表明:基于图像LSB 的秘密信息隐藏不仅具有算法简单、透明性良好的特点,而且还具有隐藏信息容量大、无失真提取、对文本文件具有双重隐藏功能等优点。关键词:图像;LSB;信息隐藏 中图分类号:TP391.41 文献标识码:A Secret Information Hiding Technology Based on LSB of Image YE Tian -yu (Department of Electronic Engineering,Shantou University,Shantou 515063,China) Abstract:Secret i nformation hiding technology based on i mage has broad applications.LSB algorithm is one of effective algori th ms.This paper first introduces the method of bit plane analysis,finding LSB has the bes t hiding effect,then respectively realizes three experi ments on information embedding and extracting:the first one takes an image with two different values as secret information and another two take a text file as secret information.The experi mental resul ts show that secret information hiding technology based on LSB of image not only has the advantages of si mplicity,good imperceptibility,bu t also has great hiding capaci ty,extraction wi th no distortion,dual hiding effect for text file,and so on. Key words:image;LSB;information hiding 0 引 言 信息隐藏是在网络环境下把机密信息隐藏在其它无关紧要的信息中形成隐秘信道,除通信双方以外的任何第三方并不知道秘密通信这个事实的存在。它与密码加密方法进行的保密通信相比,信息加密从/看不懂0变为/看不见0,转移了攻击者的目标。这种技术具有良好的透明性,即隐藏了秘密的信息,即使你看见了、听见了还是不能有所发现,加入了秘密信息的媒体在传输时与未加入秘密信息前的媒体没有什么不同。 信息隐藏技术发展到今天,已出现了各式各样的隐藏方法,最低有效位算法LSB(least significant bit )是一种简单而通用的信息隐藏方法,大部分的多媒体文件(如图像、音频和视频文件等)都可作为LSB 算法中 的载体。本文讨论以图像为载体的LSB 隐藏算法。 1 数字图像的最不重要位嵌入方法 (1)位平面分析。 根据图像采用的颜色数,可以将图像分为2位、8 位、24位图。图像中每个像素点的颜色值的某一位共同构成的一个新的二值图像就称为该图像的一个位平面图像。可以定义图像的第0个位平面到第7个位平面依次为最不重要位平面(LSB)到最重要位平面。一个8位灰度图graylena.tif 的8个位平面分解如图1所示。注意到较高位(特别是前四位),包含了大多数在视觉上很重要的数据,其他位平面对图像中更多的微笑细节有作用。把数字图像分解成位平面对于分析每一位在图像中的相对重要性是非常有用的。可以看出,图像的位平面存在一定的规律 计算机与现代化 2005年第11期 JISUANJI YU XIANDAIHUA 总第123期
信源熵实验程序： clc; close all; clear; linwidd=1 fontt=20 p0=0; pd=1; N=20 p=linspace(p0,pd,N); I=-log2(p); plot(p,I,'k'); title('I=-log2(p)函数图'); xlabel('p');ylabel('I'); clc; close all; clear; linwidd=1 fontt=20 p0=0; pd=1; N=20 p=linspace(p0,pd,N); H=-p.*log2(p)-(1-p).*log2(1-p); plot(p,H,'k'); title('H=-p.*log2(p)-(1-p).*log2(1-p)函数图'); xlabel('p');ylabel('H'); 信道容量实验程序： clc; close all; clear; linwidd=1 fontt=20 p0=0; pd=1; N=20 p=linspace(p0,pd,N); r=4 c=log2(r)+(1-p).*log2(1-p)+p.*log2(p/(r-1)); plot(p,c,'k'); title('强对称信道容量数值模拟图');
有噪信道编码--费诺不等式程序：结果图clc;close all;clear; r=3;p0=0.00001;pd=0.99999;N=2000; p=linspace(p0,pd,N); q=1-p; H=-p.*log2(p)-q.*log2(q); hold on HH=H+p.*log2(r-1) title('费诺不等式示意图');box on xlabel('PE'); ylabel('H(X/Y)'); plot(p,HH,'k:') hold on hold on fill([p,1],[HH,0],[0.6,0.6,0.6]) stem((r-1)/r,1.59,'--.r') text(0.66,1.6,'最大值') 香农编码程序： clc;clear all;close all; p=[0.2 0.19 0.18 0.17 0.15 0.1 0.01]; if sum(p)<1||sum(p)>1 error('输入概率不符合概率分布') end [p index]=sort(p,'descend'); n=length(p); pa=zeros(n,1); for ii=2:n pa(ii)=pa(ii-1)+p(ii-1); end k=ceil(-log2(p));%码字长度计算 c=cell(1,n);%生成元胞数组，用来存不同长度的码字 for ii=1:n c{ii}=''; tmp=pa(ii); for jj=1:k(ii) tmp=tmp*2; if tmp>=1 tmp=tmp-1; %c{ii}{jj}='1'; c{ii}=[char(c{ii}),'1']; else %c{ii}{jj}='0'; c{ii}=[char(c{ii}),'0']; end end end c(index)=c;%换回原来的顺序 codelength=zeros(1,n);%码长初始化 for ii=1:n fprintf(['第',num2str(ii),'个消息对应为']); disp(c{ii});%显示码字 codelength(ii)=length(c{ii});% end n_average=sum(codelength.*p) %平均码长 fprintf('平均码长为');disp(n_average); H=-sum(p.*log2(p)); fprintf('信源熵');disp(H); x=H/(n_average.*log2(2)) fprintf('编码效率');disp(x); figure h=stem(1:n,codelength);% axis([0 n+1 0 n+1]); set(h,'MarkerFaceColor','blue','linewidth',2) 实验结果 第1 个消息对应为000 第2个消息对应为001 第3个消息对应为011 第4个消息对应为100 第5个消息对应为101 第6个消息对应为1110 第7个消息对应为1111110 n_average = 3.1400 平均码长为 3.1400 信源熵 2.6087 x =0.8308 编码效率 0.8308
1.Subscript indices must either be real positive integers or logicals 中文解释：下标索引必须是正整数类型或者逻辑类型 出错原因：在访问矩阵（包括向量、二维矩阵、多维数组，下同）的过程中，下标索引要么从 0 开始，要么出现了负数。注：matlab 的语法规定矩阵的索引从1 开始，这与 C 等编程语言的习惯不一样。 解决办法：自己调试一下程序，把下标为 0 或者负数的地方修正。 2.Undefined function or variable "U" 中文解释：函数或变量 U 没有定义 .出错原因及解决办法：可能变量名输入错误，仔细检查 3.Matrix dimensions must agree 中文解释：矩阵的维数必须一致 出错原因：这是由于运算符(= + - / * 等)两边的运算对象维数不匹配造成的，典型的出错原因是错用了矩阵运算符。matlab 通过“.”来区分矩阵运算和元素运算。 解决办法：自己调试一下程序，保证运算符两边的运算对象维数一致。 4.Function definitions are not permitted at the prompt or in scripts 中文解释：不能在命令窗口或者脚本文件中定义函数 出错原因：一旦在命令窗口写 function c = myPlus(a,b)，此错误就会出现，因为函数只能定义在 m 文件中。关于脚本文件和 m 文件的区别请查阅 matlab 基础书。简言之： 1) 如果你写成 function 的形式，那么必须写在 m 文件中，且以 function 开头（即 function 语句前不能包含其他语句，所有语句必须放在 function 中，当然，function 的定义可以有多个，各 function 之间是并列关系，不能嵌套）； 2) 如果你写成脚本的形式，则既可以写在命令窗口中，也可以写在 m 文件中，但两者均不能包含 function 语句（即不能进行函数的定义） 解决办法：新建一个 m 文件，然后再进行函数的定义 5.One or more output arguments not assigned during call to '...' 中文解释：在调用...函数过程中，一个或多个输出变量没有被赋值 出错原因：函数如果带有输出变量，则每个输出在返回的时候都必须被赋值。容易出现这个错误的两个地方是：在部分条件判断语句（如 if ）中没有考虑到输出变量的返回值。 在循环迭代过程中部分变量的维数发生了变化。 解决办法：调试程序，仔细查看函数返回时各输出变量的值。更好的方法是：在条件判断或者执行循环之前对所使用的变量赋初值。 6.??? Index exceeds matrix dimensions.
《信息隐藏技术》实验 实验三：基于LSB的信息隐藏算法 学生姓名:学号： 学院: 计算机学院 专业: 信息安全 班级: 指导教师： 2015年 12 月 16 日
目录 1 实验目的 (2) 2 实验环境及内容 (2) 3 实验原理 (2) 4 实验步骤 (3) 5 实验思考和总结 (12)
1实验目的 该实验为验证性实验。目的是通过实验使学生掌握经典隐藏算法，在MATLAB环境下，编写基于图像的LSB信息隐藏算法程序。用MATLAB 函数实现LSB信息隐藏，并进行分析。 2实验环境及内容 ●实验环境 安装MATLAB7.0的计算机 ●实验内容 首先学会提取图片的八个位平面，然后再将信息隐藏在最低位平面中。 3实验原理 LSB(least significant bit)算法是最早提出的一种典型的空间域信息隐藏算法。它使用特定的密钥通过伪随机序列发生器产生随机信号，然后按一定的规则排列成二维水印信号，并逐一插到原始图像相应像素值的最低几位。由于水印信号隐藏在最低位，相当于叠加了一个能量微弱的信号，因此在视觉和听觉上很难察觉。作为大数据量的信息隐藏方法，LSB在隐藏通信中仍占据相当重要的地位。 隐秘算法核心是将我们选取的像素点的最不重要位依次替换成秘密信息,以达到信息隐秘的目的。嵌入过程包括选择一个图像载体像
素点的子集{j1,…,jl(m)},然后在子集上执行替换操作像素cji←→mi,即把cji的LSB与秘密信息mi进行交换(mi可以是1或0)。一个替换系统也可以修改载体图像像素点的多个比特,例如,在一个载体元素的两个最低比特位隐藏两比特、三比特信息,可以使得信息嵌入量大大增加但同时将破坏载体图像的质量。在提取过程中,找出被选择载体图像的像素序列,将LSB(最不重要位)排列起来重构秘密信息,算法描述如下: 嵌入过程:for(i=1;i<=像素序列个数;i++) si←ci for(i=1;i<=秘密消息长度;i++) //将选取的像素点的最不重要位依次替换成秘密信息 sji←cji←→mi 4实验步骤 提取位平面源代码 a = imread('1.jpg'); b = rgb2gray(a); %将彩色图像转换为灰度图像 c = bitand(b, 1); figure(1); imshow(c); d = bitand(b, 2);
第一章插值方法 1.1Lagrange插值 1.2逐步插值 1.3分段三次Hermite插值 1.4分段三次样条插值 第二章数值积分 2.1 Simpson公式 2.2 变步长梯形法 2.3 Romberg加速算法 2.4 三点Gauss公式 第三章常微分方程德差分方法 3.1 改进的Euler方法 3.2 四阶Runge-Kutta方法 3.3 二阶Adams预报校正系统 3.4 改进的四阶Adams预报校正系统 第四章方程求根 4.1 二分法 4.2 开方法 4.3 Newton下山法 4.4 快速弦截法 第五章线性方程组的迭代法 5.1 Jacobi迭代 5.2 Gauss-Seidel迭代 5.3 超松弛迭代 5.4 对称超松弛迭代 第六章线性方程组的直接法 6.1 追赶法 6.2 Cholesky方法 6.3 矩阵分解方法 6.4 Gauss列主元消去法
第一章插值方法 1.1Lagrange插值 计算Lagrange插值多项式在x=x0处的值. MATLAB文件：（文件名：Lagrange_eval.m）function [y0,N]= Lagrange_eval(X,Y,x0) %X,Y是已知插值点坐标 %x0是插值点 %y0是Lagrange插值多项式在x0处的值 %N是Lagrange插值函数的权系数 m=length(X); N=zeros(m,1); y0=0; for i=1:m N(i)=1; for j=1:m if j~=i; N(i)=N(i)*(x0-X(j))/(X(i)-X(j)); end end y0=y0+Y(i)*N(i); end 用法》X=[…];Y=[…]; 》x0= ; 》[y0,N]= Lagrange_eval(X,Y,x0) 1.2逐步插值 计算逐步插值多项式在x=x0处的值. MATLAB文件：（文件名：Neville_eval.m）function y0=Neville_eval(X,Y,x0) %X,Y是已知插值点坐标 %x0是插值点 %y0是Neville逐步插值多项式在x0处的值 m=length(X); P=zeros(m,1); P1=zeros(m,1); P=Y; for i=1:m P1=P; k=1; for j=i+1:m k=k+1;
实验一 MATLAB 的基本操作 一、实验目的 1、掌握Matlab 软件使用的基本方法； 2、熟悉Matlab 的数据表示、基本运算方法； 3、熟悉Matlab 绘图命令及基本绘图控制。 二、实验仪器与软件 1、PC 机 1台 2、MATLAB7.0环境 三、实验原理 MATLAB 环境是一种为数值计算、数据分析和图形显示服务的交互式的环境。MATLAB 有3种窗口，即：命令窗口（The Command Window ）、m-文件编辑窗口（The Edit Window ）和图形窗口（The Figure Window ），而Simulink 另外又有Simulink 模型编辑窗口。 1、命令窗口（The Command Window ） 当MATLAB 启动后，出现的最大的窗口就是命令窗口。用户可以在提示符“>>”后面输入交互的命令，这些命令就立即被执行。 在MATLAB 中，一连串命令可以放置在一个文件中，不必把它们直接在命令窗口内输入。在命令窗口中输入该文件名，这一连串命令就被执行了。因为这样的文件都是以“.m ”为后缀，所以称为m-文件。 2、m-文件编辑窗口（The Edit Window ） 我们可以用m-文件编辑窗口来产生新的m-文件，或者编辑已经存在的m-文件。在MATLAB 主界面上选择菜单“File/New/M-file ”就打开了一个新的m-文件编辑窗口；选择菜单“File/Open ”就可以打开一个已经存在的m-文件，并且可以在这个窗口中编辑这个m-文件。如，计算3 [(12)34)]2+?-÷，只需在提示符“>>”后输入“((1+2)*3-4)/2^3”，然后按Enter 键。该命令行涉及加、减、乘、除及幂运算符，MATLAB 运算的执行次序遵循的优先规则为：从左到右执行；幂运算具有最高的优先级，乘法和除法具有相同的次优先级，加法和减法有相同的最低优先级；使用括号可以改变前述优先次序，并由最内层括号向外执行。 3、常用命令 MATLAB 可以把多条命令放在同一行，各命令之间用逗号“，”或分号“；”隔开，逗号告诉MATLAB 系统显示本命令的结果，分号告诉系统取消结果的显示（只是不显示，该命令行仍正常执行）。 clear 命令 若想清除MATLAB 当前工作空间中所有定义过的变量，使用clear 命令；若只想清除其中某几个变量，只需在命令clear 后写入变量的名称即可，如想清除变量x 、y ，在“>>”后键入“clear x y ”即可。 help 命令 在MA TLAB 命令窗口使用help 命令寻求帮助。例如， >> help ceil ? lookfor 命令 lookfor 命令来查询根据用户提供的关键字搜索到的相关函数，然后再使用help 命令与其配合使用，可了解查询到的确切函数的具体用法。 常用的数学函数如表1所示。 表1 常用数学函数表 函数 功能 函数 功能 sin(x) 正弦函数 asin(x) 反正弦函数
Matlab常见语法错误 Subscript indices must either be real positive integers or logicals 中文解释：下标索引必须是正整数类型或者逻辑类型 出错原因：在访问矩阵（包括向量、二维矩阵、多维数组，下同）的过程中，下标索引要么从0 开始，要么出现了负数。注：matlab的语法规定矩阵的索引从1 开始，这与C 等编程语言的习惯不一样。 解决办法：自己调试一下程序，把下标为0 或者负数的地方修正。关于矩阵的各种访问方式。 举例： 【错误代码】： 1.for s = 0:0.2:2 2. a(s) = 4 * s - 1; 3.end 复制代码 【正确代码】： 1. 2.for s = 1:10 3. a(s) = 4 * s - 1; 4.end 复制代码 当然，这样的赋值使用循环是低效的，因此不提倡
Undefined function or variable "a" 中文解释：函数或变量 a 没有定义 出错原因及解决办法： 1. 如果 a 是变量，则参见[原创]写给学习 matlab 的新手们一贴中技巧 6 后半部分 2. 如果 a 是函数，则这是因为matlab在所有已添加的路径中都无法找到该函数对应的m 文件而导致的，对此，把该m文件移动到当前路径下再运行便可。（注：何谓“所有已添加的路径”——在命令窗口中输入 pathtool 可以查看；何谓“当前路径”——在命令窗口中输入 pwd 可以查看） 举例： myPlus.m： 1.function d = myPlus(a,b) 2. d = a + b; 复制代码 【错误调用】： 命令窗口： 1.>> z = myPlus(a,b) 复制代码 【正确调用】： 命令窗口： 1.>> x = 2; 2.>> y = 5; 3.>> z = myPlus(x,y) 复制代码
实验报告 课程名称：信息论与编码姓名： 系：专 业：年 级：学 号：指导教 师：职 称：
年月日 目录 实验一信源熵值的计算 (1) 实验二Huffman 信源编码. (5) 实验三Shannon 编码 (9) 实验四信道容量的迭代算法 (12) 实验五率失真函数 (15) 实验六差错控制方法 (20) 实验七汉明编码 (22)
实验一信源熵值的计算 、实验目的 1 进一步熟悉信源熵值的计算 2 熟悉Matlab 编程 、实验原理 熵(平均自信息)的计算公式 q q 1 H(x) p i log2 p i log2 p i i 1 p i i 1 MATLAB实现：HX sum( x.* log2( x))；或者h h x(i)* log 2 (x(i )) 流程：第一步：打开一个名为“ nan311”的TXT文档，读入一篇英文文章存入一个数组temp，为了程序准确性将所读内容转存到另一个数组S，计算该数组中每个字母与空格的出现次数( 遇到小写字母都将其转化为大写字母进行计数) ，每出现一次该字符的计数器+1；第二步：计算信源总大小计算出每个字母和空格出现的概率；最后，通过统计数据和信息熵公式计算出所求信源熵值(本程序中单位为奈特nat )。 程序流程图： 三、实验内容 1、写出计算自信息量的Matlab 程序 2、已知：信源符号为英文字母(不区分大小写)和空格输入：一篇英文的信源文档。输出：给出该信源文档的中各个字母与空格的概率分布，以及该信源的熵。 四、实验环境 Microsoft Windows 7
五、编码程序 #include"stdio.h" #include
MATLAB常用指令 1．常用命令-->管理命令和函数 addpath 添加目录到MATLAB搜索路径 doc 在Web浏览器上现实HTML文档 help 显示Matlab命令和M文件的在线帮助 helpwin helpdesk help lookfor 在基于Matlab搜索路径的所有M文件中搜索关键字partialpath 部分路径名 path 所有关于路径名的处理 pathtool 一个不错的窗口路径处理界面 rmpath 删除搜索路径中指定目录 type 显示指定文件的内容 ver 版本信息 version 版本号 web 打开web页 what 列出当前目录吓所有的M文件Mat文件和Mex文件whatsnew 显示readme文件 which 显示文件位置 （返回） 2．常用命令-->管理变量和工作区(输入输出、内存管理等) clear 从内存中删除 disp 显示文本或数组内容 length 数组长度(最长维数) load 重新载入变量(从磁盘上) mlock 锁定文件，防止文件被错误删除 munlock 解锁文件 openvar 在数组编辑器中打开变量 pack 整理内存空间 save 保存变量到文件8*) size 数组维数 who 列出内存变量 whos 列出内存变量,同时显示变量维数 workspace 显示工作空间窗口 （返回） 3．常用命令-->管理命令控制窗口(command窗口) clc 清空命令窗口 echo 禁止或允许显示执行过程 format 设置输出显示格式 home 光标移动到命令窗口左上角 more 设置命令窗口页输出格式 （返回）