工程热力学上机实验报告

第一篇：工程热力学上机实验报告

工程热力学上机实验报告

姓名：

专业：能源与动力

学号：朗肯循环

对蒸汽动力循环的基本循环——朗肯循环，其工作原理是，从锅炉出来的高温T1，高压p1的过热水蒸气经汽轮机绝热膨胀做工至低压p2的乏汽，在冷凝器中凝结成饱和液体，经水泵升压至p1下的未饱和过冷液体，进入锅炉加热至过热蒸汽，再进入汽轮机绝热膨胀做功，周而复始的将热能转换为机械能，图1为理想朗肯循环工作过程的水蒸气T—S图。循环中：

工质在锅炉中的加热量：q1=h1-h4 在冷凝器中的放热量：q2=h2-h3 在汽轮机中的做功量：w1=h1-h2 在水泵中的耗功量：w2=h4-h3 循环热效率n=(q1-q2)/q1=1-(h2-h3)/(h1-h4)如忽略泵功，h3=h4 则循环效率：n=(h1-h2)/(h1-h3)=1-(h2-h3)/(h1-h3)程序如右图

1、实验初始参数：P1=4Mpa,t1=450摄氏度，p2=6kpa如下图

P1,P2不变改变初温t1，结果如图

t1,P2不变改变初温p1，结果如图

t1,P1不变改变初温p2，结果如图

当汽机绝热内效率nex=0.85，此时实际不可逆循环的热效率nt=0.33226925

二，再热循环

循环热效率：n=((h1-hb)+(ha-h2)-(h4-h3))/((h1-h4)+(ha-hb))如忽略泵功，h3=h4 则循环效率：n=((h1-hb)+(ha-h2))/((h1-h3)+(ha-hb))如果过分提高压力p1，而不响应提高t1,将引起乏汽敢赌x2减小，产生不利后果。为此，将新蒸汽膨胀至某一中间压力pb后撤出汽轮机，导入锅炉中的特设的再热器或其他换热设备中，使之再加热后，再倒入汽轮机继续膨胀至背压p2.即为再热循环。图2为再热循环工作过程的T-S图。程序如右下面3个图

2、初始参数

P1=12Mpa,t1=500摄氏度，终压p2=6kpa, 第一个再热循环的再热压力pa1=3Mpa, 另一个再热循环的压力为pa2=0.6Mpa

顺序1-3

2、初始参数

P1=12Mpa,t1=500摄氏度，终压p2=6kpa, 第一个再热循环的再热压力pa1=3Mpa, 另一个再热循环的压力为pa2=0.6Mpa 结果如右图

经程序运行，比较朗肯循环和在再热循环的效率发现：相同初始参数和相同背压下，再热循环效率较高。

当汽机绝热内效率nex=0.85，此时实际不可逆循环的热效率nt=0.3812862 0.3777485三、一次抽气回热循环

为提高循环平均吸热温度，在汽轮机膨胀做功过程中，取出一部分蒸汽用以回热给水，使循环平均吸热温度有所提高，而平均官方热温度不变，热效率相应也随之提高。

如图为一级抽气回热循环的T-S图。循环热效率：

n=((h1-ha)+(1-a)(ha-h2)-(h6-h5)-(1-a)(h4-h3))/(h1-h6)若忽略泵功，循环效率：

n=((h1-ha)+(1-a)(ha-h2))/(h1-h5)其中a由图回热器的热平衡方程求的: a=(h5-h3)/(ha-h3)=(ha’-h2’)/(ha-h2’)程序如右图

1、实验初始参数： P1=5Mpa,t1=450摄氏度

抽气压力pa=0.5Mpa，p2=6kpa 实验结果如图

当汽机绝热内效率nex=0.85，此时实际不可逆循环的热效率nt=0.72267255

四、心得

经过这两次的工程热力学上机实验，发现可以用电脑来简化计算步骤。巩固了所学热力学知识，掌握了系统热力计算的基本概念和计算方法，得到了所求问题的计算结果，但是这只是一些简单的题目，在更复杂的问题上可更上一层楼。

第二篇：上机实验报告格式

请各位同学按照下列格式写：否则按照0分记

上机实验报告格式

Matlab学习第？次上机实验报告（第？次）

姓名：？？？班级:???上机时间：？？？？？？

1.上机内容

2.建模及算法分析

3.程序

4.小结

5.参考文献

将文档发到我的邮箱里面

learn_matlab@163.com

第三篇：上机实验报告

一． 题目1． 建立一个学生档案，内容包括学号，姓名，年龄，性别，数学，物理和英语3门功课成绩。要求实现以下功能：1)数据输入；2)查询某个学生的成绩；3)按平均排列输出；4)统计某门课各分数段人数；5)删除某个学生记录；2． 编程实现对二位数进行加，减，乘运算，每运行一次程序做10道题，完成后给出成绩（每题10分）。二． 设计思想和模块划分1.1．先定义所需要的条件，例如，姓名，学号，性别，三门功课的成绩，平均分等。2．输入姓名，学号，性别，年龄，三门功课的成绩后，就输出相应的姓名，学号，性别，年龄，三门功课的成绩。3．打印表头；4．在求出平均分，并打印出平均分。5．在查找学生，并打印出学生的相关资料；并按照平均分的高低排列；6．在统计学生各科成绩是否》=60分，如果是就是通过，否就是没通过，最后打印出最后的结果。7．删除学生，选中删除的学生，如果确定要，那么该学生的一切记录就变为0；2．1．随机调用函数，产生两个其值为10到99的随机整数分别放在c,d中用作运算时的运算数。2．随机调用函数，产生一个值为1到3 的随机整数放入b中，用来选择不同的运算。根据b的值选择不同的运算。当b=1时，去进行c+d的运算练习；当b＝2 时，去进行c-d的运算练习；当b=3时，去进行c×d的运算练习。每完成一个算题就给出是否正确的信息。3．根据答案正确与否统计。正确时，显示’right’;错误时，显示’wrong’。4．步骤1.2.3.4.重复10次。5．输出学生成绩。三． 运行结果1.Input the 1 student :num:03name:wfage:18sex:fmath:98phy:75eng:65Input the 2 student :num:06name:scage:17sex:fmath:78phy:45eng:65Input the 3 student :num:09name:ytage:17sex:mmath:69phy:75eng:501.readin 2.finds 3.del 4.tj 5.exit1-----------------------------table----num name age sex math phy eng ave------3 wf 18 f 98.0 75.0 65.0 79.3------9 yt 17 m 69.0 75.0 50.0 64.7-------6 sc 17 f 78.0 45.0 65.0 62.7-------1.readin 2.finds 3.del 4.tj 5.exit2Into number:3---num name age sex math phy eng ave---3 wf 18 f 98.0 75.0 65.0 79.3----1.readin 2.finds 3.del 4.tj 5.exit3Into number:6---num name age sex math phy eng ave---6 sc 17 f 78.0 45.0 65.0 62.6---true? Y/Ny1----------------table------------------num name age sex math phy eng ave----3 wf 18 f 98.0 75.0 65.0 79.3-----9 yt 17 m 69.0 75.0 50.0 64.7-----0 0 0 f 0.0 0.0 0.0 0.0-----1.readin 2.finds 3.del 4.tj 5.exit4 1.input the math 2.input the phy 3.input the eng 1math: pass num is :3 no pass num is :02.32+82=114right!97*15=1455right!99-58=41right!76-66=10right!58+13=71right!37+49=86right!99+52=151right!49*69=5462wrong!73-38=31wrong!86+34=120right!mark is 80四． 主要错误改正方法1.在本题的编写过程中常出现begin 和end 不配对，只要始它们的个数相同就可以了，就可以了。在程序中常出现的’unknown identifier’没有定义的错误，只要在程序前加上定义，就可以了。还有就是type mismatch’的类型不匹配，只要根据前后内容，重新定义，就可以了。2.在本题编写中，题目要求是两个两位数运算，如果写成c:=random(100)在运行过程中，就会出现一位数运算，只要把它写成c:=random(90)+10，就可以了。五． 实习小结1．第一题的程序编写比较复杂，需要较多的定义，因而在后面的程序部分，就常出现标识符重复，而导致的类型不匹配。还有就是每个部分编好后，上下不能连接，而无法调用。2．第二题较第一题简单，和书上的例题类似，所以，比较简单，但从中也出现一些不问题。以上两大题中出现的问题，要在同学的帮助下，才能完成。总的来说两题程序比较困难六． 程序清单1.program twins;const m=3;n=3;typestudent=recordnum:integer;name:string[3];age:integer;sex:char;s:array[1..n] of real;math,eng,phy:real;ave:real;end;sarr=array[1..m] of student;ta=array [1..m] of student;f=string[6];varstu:sarr;ct:ta;k,d:integer;procedure readin(var stu:sarr);const wrong='0<=data<=100 ,again!';var j,i:integer;a:student;beginwriteln;writeln('Input ',m,' num name sex math phy eng',' of student.');for j:=1 to m dobeginwriteln;writeln('Input the ',j,' student :');with a dobeginwrite('num:');readln(num);write('name:');readln(name);write('age:');readln(age);write('sex:');readln(sex);while(sex<>'m')and(sex<>'f')dobeginwrite('sex:');readln(sex)end;for i:=1 to n do begincase i of1:begin write('math:');readln(s[i]);while(s[i]<0)or(s[i]>100)dobeginwriteln(wrong);write('math:');readln(s[i]);end;end;2:begin write('phy:');readln(s[i]);while(s[i]<0)or(s[i]>100)dobeginwriteln(wrong);write('phy:');readln(s[i]);end;end;3:begin write('eng:');readln(s[i]);while(s[i]<0)or(s[i]>100)dobeginwriteln(wrong);write('eng:');readln(s[i]);end;end;end;end;stu[j]:=a;end;writeln;end;end;procedure ave(var stu:sarr);var j,i:integer;k:real;beginfor j:=1 to m dobegink:=0;with stu[j] dobeginfor i:=1 to n do k:=k+s[i];ave:=k/nend;end;end;procedure px(var stu:sarr);var tm:ta;j,i,p:integer;beginfor j:=1 to m dobeginp:=j;for i:=j+1 to m doif stu[i].ave>stu[p].ave then p:=i;tm[1]:=stu[p];stu[p]:=stu[j];stu[j]:=tm[1];end;end;procedure head(ct:ta);var i:integer;beginwrite(' ');for i:=1 to 72 do write('-');writeln;write('',' num':4,' ','name':9,' ','age':8,' ','sex':5,' ','math':6,' ','phy':9,' ','eng':8,' ');writeln('ave':6,' ');write(' ');for i:=1 to 72 do write('-');writeln;with ct[1] dobeginwrite('',num:7,' ',name:8,' ',age:8,' ',sex:5,' ');for i:=1 to n do write(s[i]:7,' ');write(ave:10,' ');writeln;end;write(' ');for i:=1 to 72 do write('-');writeln;end;procedure find(var stu:sarr);var n,i:integer;u:boolean;beginu:=true;while u dobeginwrite('Into number:');readln(n);for i:=1 to m dobeginif stu[i].num=n then beginct[1]:=stu[i];u:=false;d:=iend;end;end;end;procedure del(var stu:sarr);var i:integer;b:char;beginhead(ct);write('true? Y/N');readln(b);if b='y' then begin write('1');with stu[d] dobeginnum:=0;name:='0';age:=0;sex:='f';for i:=1 to n do s[i]:=0;ave:=0end;endelseend;procedure tj(var stu:sarr);var w1,q1,q2,num,j:integer;h1,h2,h3:boolean;beginh1:=false;h2:=false;h3:=false;q1:=0;q2:=0;writeln(' ':30,'1.input the math 2.input the phy 3.input the eng ');readln(w1);case w1 of1:h1:=true;2:h2:=true;3:h3:=true;end;if h1 thenbeginbeginfor j:=1 to m dowith stu[j] doif stu[j].math>=60 thenq1:=q1+1elseq2:=q2+1;writeln('math:');end;if h2 thenbeginbeginfor j:=1 to m dowith stu[j] doif stu[j].phy>=60 thenq1:=q1+1elseq2:=q2+1;writeln('phy:');end;if h3 thenbeginbeginfor j:=1 to m dowith stu[j] doif stu[j].eng>=60 thenq1:=q1+1elseq2:=q2+1;writeln('eng:');end;end;end;end;writeln(' ':5,'pass num is :',q1);writeln(' ':5,'no pass num is :',q2);end;procedure head1(var stu:sarr);var j,i:integer;beginwriteln('----------------table--------------':60);writeln;write(' ');for i:=1 to 72 do write('-');writeln;write('',' num':4,' ','name':9,' ','age':8,' ','sex':5,' ','math':6,' ','phy':9,' ','eng':8,' ');writeln('ave':6,' ');write(' ');for i:=1 to 72 do write('-');writeln;for j:=1 to m do beginwith stu[j] dobeginwrite('',num:7,' ',name:8,' ',age:8,' ',sex:5,' ');for i:=1 to n do write(s[i]:7:1,' ');write(ave:10:1,' ');writeln;end;write(' ');for i:=1 to 72 do write('-');writeln;end;end;beginreadin(stu);ave(stu);k:=0;while k<>5 dobeginwriteln('1.readin 2.finds 3.del 4.tj 5.exit');read(k);readln;while(k<>1)and(k<>2)and(k<>3)and(k<>4)and(k<>5)dobeginwriteln(' mistake,please again');read(k)end;case k of1: begin px(stu);head1(stu)end;2: begin find(stu);head(ct)end;3: begin find(stu);del(stu);px(stu);head1(stu)end;4: begin tj(stu);end;5: writeln('exit!');end;end;end..2.program shadow;var a,b,c,d,e:integer;right:boolean;procedure s1;beginwriteln('right!');right:=trueend;procedure s2;beginwriteln('wrong!');right:=falseend;procedure s3(x,y:integer);var h:integer;beginwrite(x,'+',y,'=');readln(h);if h=x+y then s1else s2 end;procedure s4(x,y:integer);var h,w:integer;beginif x

第四篇：材料热力学实验报告

材料热力学实验报告

姓名：XXX 学号：201401049 时间：2014.12.25

实验题目：Al-Si合金DSC-TG及热膨胀分析

一、实验目的

1.掌握DSC（差示扫描量热）（功率补偿型）、TGA（热重分析）和TMA实验原理及测量方法；

2.学会使用上述三种实验仪器，并能够绘制样品的DSC、TGA曲线、计算实验材料线膨胀系数；

3.掌握实验结果的处理方法，并能够对谱图进行分析。

二、实验仪器及实验样品

同步热分析仪 型号：NETZSCH STA 409PC/PG

制造厂家：德国耐弛公司

热机械分析仪 型号：TMA/SDTA840

制造厂家：瑞士梅特勒-托利多公司 样品：Al-Si合金

三、实验原理 1.DSC原理

差示扫描量热仪DSC是在程序控温下，测量物质和参比物之间的能量差随温度变化关系的一种技术。根据测量方法的不同，又分为功率补偿型DSC和热流型DSC两种类型。常用的功率补偿DSC是在程序控温下，使试样和参比物的温度相等，测量每单位时间输给两者的热能功率差与温度的关系的一种方法。DSC是在控制温度变化情况下，以温度(或时间)为横坐标，以样品与参比物间温差为零所需供给的热量为纵坐标所得的扫描曲线。当试样在加热过程中由于热效应与参比物之间出现温差ΔT时，通过差热放大电路和差动热量补偿放大器，使流入补偿电热丝的电流发生变化，当试样吸热时，补偿放大器使试样一边的电流立即增大；反之，当试样放热时则使参比物一边的电流增大，直到两边热量平衡，温差ΔT消失为止。换句话说，试样在热反应时发生的热量变化，由于及时输入电功率而得到补偿，所以实际记录的是试样和参比物下面两只电热补偿的热功率之差随时间t的变化的关系。2.TGA原理

是指在程序控制温度下测量待测样品的质量与温度变化关系的一种热分析技术，用来研究材料的热稳定性和组份。TGA在研发和质量控制方面都是比较常用的检测手段。热重分析在实际的材料分析中经常与其他分析方法连用，进行综合热分析，全面准确分析材料。

3.同步热分析（TG-DSC）将热重分析与示差热分析结合测量。4.热膨胀原理

物体的体积或长度随温度的变化而发生相对变化的现象称为热膨胀；本质是晶体点阵结构间的平均距离随温度变化而变化。

一般的普通材料，通常所说膨胀系数是指线膨胀系数，其意义是温度升高1℃时单位长度上所增加的长度，单位为厘米/厘米·度。假设物体原来的长度

为L0，温度升高后长度的增加量为L，它们之间存在如下关系：

L╱L0＝1t

（1）

式中，1称为线膨胀系数，也就是温度每升高1℃时，物体的相对伸长。当物体的温度从T1上升到T2时，其体积也从V1变化为V2，则该物体在T1至T2的温度范围内，温度每上升一个单位，单位体积物体的平均增长量为

=(V1－V2）╱V1（T1－T2）

（2）

式中，为平均体膨胀系数。从测试技术来说，测体膨胀系数较为复杂。因此，在讨论材料的热膨胀系数时，常常采用线膨胀系数：

＝（L1－L2）╱L1（T1－T2）

(3)式中：―试样的平均线膨胀系数；L1―在温度为T1时试样的长度； L2―在温度为T2时试样的长度；

必须指出，由于膨胀系数实际上并不是一个恒定的值，而是随温度变化的，所以上述膨胀系数都是具有在一定温度范围t内的平均值的概念，因此使用时要注意它适用的温度范围，如：(0～300)= 5.7×10-7 / k。本实验主要测量材料的平均线膨胀系数。

四、实验参数

DSC-TGA

参数: 温度范围：40-740℃

加热速率：20K/min

吹扫气体：氩气 20mL/min

坩埚：DSC/TG pan Al2O

3样品量：15.434mg 热膨胀参数： 铝硅合金材料厚度为2.6342mm

温度范围：30-300℃

升温速率：10℃/min

五、实验步骤

1、DSC-TGA实验

a.开机 打开恒温水浴、STA449C主机、TASC414/4控制器与计算机电源。一般在水浴与热天平打开4小时后，可以开始测试。如果配有低温系统，打开冷却控制器CC200L电源。打开Proteus测试软件。b.进行抽真空 进行3次以上抽真空操作，确保氧气排除干净，然后在通入惰性气体保护。

c.基线测试（若已有原先做好的基线文件，可跳过此步骤）d.样品测试 先将空坩埚放在天平上称重，随后将样品加入坩埚中，称取样品重量，本实验中样品重量为15.434mg。将装有样品的坩埚放入炉体内，关闭炉体，点击“编辑”菜单下的“测量向导”，在弹出的“STA449C测量向导”对话框中选择测量模式为“样品+修正”，输入样品名称、编号与样品质量，设定完成后点击“继续”，打开温度校正文件，打开灵敏度校正文件，编辑设定温度程序，本实验温度范围为40~740℃。设定测量文件名，点击“初始化工作条件”，软件将根据实验设置自动打开各路气体。转动流量计上的旋钮，调节各路气体流量。本实验所用吹扫气体为氩气，流速为20ml/min。随后点击“诊断”菜单下的“炉体温度”与“查看信号”，调出相应的显示框，点击“清零”，对天平进行清零。随后观察仪器状态满足如下条件：

1）炉体温度、样品温度、TG信号及DSC信号处于稳定状态； 2）样品温度与设定起始温度相吻合。即可点击“开始”开始测量。

2、热膨胀实验

a.开机启动TMA系统，进入主程序；

b.设置起始温度、终止温度、升温速率，输入新建文件名开始运行；

c.达到预定温度后，打开炉子，归零，抬PROBE，加入样品，放PROBE，再次归零，得到长度，本次实验材料的长度为2.6342mm；ADJUST XU旋转炉子使POS位于0%，对中，按OK确认； d.回到编辑程序处：开始实验；终止实验。

六、实验分析及结果 1.DSC-TGA实验

从上图中观察到由548.7℃开始至562.8℃时的吸热峰，峰面积为19.05J/g，初步分析是由于铝硅合金的亚共晶转变形成的相变潜热。因为在铝硅相图中，共晶转变是发生在5771℃，而此时的562.8℃还未达到共晶转变温度，但实验用材料的成分不确定，合金中可能存在其他影响共晶转变温度变化的成分，所以推测在此温度下是发生了亚共晶转变的吸热过程。

从上图中观察到由618.3℃开始的至642.0℃时出现了很强的吸热峰，峰面积为187.6J/g，初步分析是由于铝硅合金的熔融产生的。理论中铝的熔点为660.4℃，硅的熔点更高，而任何合金的熔点都是小于合金中任意成分的熔点，所以由此推测是随着温度逐渐升高，铝硅合金达到熔点，熔融时需要吸收大量的热，从而引起明显的吸热峰。由上面两个图中还能观察到除上面分析的两个峰外，DSC曲线整体还有明显的下降趋势，初步分析是由于仪器本身的原因导致的。2.热重实验分析

由上图中可观察到TGA曲线有小幅度的上升，质量增重为1.09％，初步分析是由于在测量过程中铝硅合金被氧气氧化引起的；虽然在实验开始前已经进行三次抽真空，但是不排除有氧气混入。3.线膨胀实验分析

热（线）膨胀系数：α=斜率/材料厚度（单位：μm）由拟合图像得 斜率B=0.05971 μm/℃

长度=2.6342mm=2634.2μm



0.05971m℃2.27105/℃2634.2m

计算得到Al-Si合金的热膨胀系数为2.27×10-5/℃。

七、误差分析

1、影响DSC结果的因素

a、样品量 样品的分辨率高，但灵敏度下降，一般根据样品热效应大小调节样品量，一般3～5mg。另一方面样品量多少对所测转变温度也有影响，随样品量的增加，峰起始温度基本不变，但峰顶温度增加，峰结束温度也提高，因此如同类样品要相互比较差异，最好量相同。b、升温速率 通常升温速率范围在5～20℃/min。升温速率越快，灵敏度提高，分辩率下降。灵敏度和分辨率是一对矛盾，一般选择较慢的升温速率以保持好的分辨率，而适当增加样品量来提高灵敏度。随着升温速率的增加。融化峰起始温度变化不大，而峰顶和峰结束温度提高。峰型变宽。

c、气氛 一般使用惰性气体，如氮气、氩气、氦气等，就不会产生氧化反应峰，同时又减少试样挥发对监测器的腐蚀。气流流速必须恒定，否则会引起基线波动。气体性质对测定有显著影响，要根据不同气体的热导率和所做的实验合理选择。有时可通过比较氮气和氧气中DSC曲线，来解释一些氧化反应。

2、影响TG的因素 仪器影响

a、浮力的影响

热天平在热区中，其部件在升温过程中排开空气的重量不断减小，即浮力在减小，会产生试样表观增重。

热天平试样周围气氛受热变轻会上升，形成向上热气流，作用在热天平上相当于减重。b、坩埚影响

热分析作用的坩埚材质，要求对试样、中间产物、最终产物和气氛都是惰性的，其几何形状和重量对热分析也有影响。c、挥发物再冷却的影响

实样热分析过程溢出的挥发物有可能在热天平其他部分再冷凝，这会污染仪器、使测的的失重量偏低，待温度进一步上升，这些冷凝物可能再次挥发产生假失重，使曲线变形，使测定不准，也不能重复。为了解决这个问题可适当向热天平通适量气体。操作条件的影响 a、升温速率的影响

升温速率越大温度滞后越严重，开始分解温度Ti及终止分解温度Tf都高，温度区间宽。

一般热重法测定不要采用太高的升温速率，对传热差的高分子试样一般用5～10K/min，对传热好的无机物、金属试样可用10～20K/min，对作用动力学分析还要低一些。b、气氛的影响

热天平周围气氛的改变对TG曲线的影响也非常显著。在流动气氛中进行TG测定时，流速大小、气氛纯度、进气温度等是否稳定，对TG曲线都有影响。一般，气流速度大小，对传热和溢出气体扩散都有利。使热分解温度降低。对于真空和高压热天平，气氛压力对TG也有很大影响。c、试样用量、粒度和装填情况的影响

试样用量多时，要过较长时间内部才能达到分解温度。试样粒度与用量对TG影响相似，粒度越小、反映面积越大，反映越容易进行，使TG曲线的Ti和Tf都低，反应区间也窄。试样装填情况首先要求颗粒均匀，必要时要过筛。

第五篇：信息论上机实验报告

华中科技大学文华学院

信息论上机实验报告

专业年级：10级通信工程** 学号：1001******** 姓名：** 指导教师：***

2012年 04月30日

目录 等长编码.........................2 1.1 求信源的熵...................2 1.2 求等长码码长.................2 1.3 求等长长等编码...............3 1.4 等长编码全部源代码...........3 1.5MATLAB仿真及其截图..........5 2霍夫曼编码.......................6 2.1 霍夫曼编码源代码.............6 2.2MATLAB仿真及其截图..........9 3游程编码.........................9 3.1 游程编码源代码...............9 3.2 MATLAB仿真及其截图..........11 4 学习心得........................12 5 参考书籍........................14

for i=1:1: 4;s=i;t=lc;for j=1:1:lc;if s>power(2,t)/2;A(i,lc-t+1)=1;s=s-power(2,t)/2;else A(i,lc-t+1)=0;end t=t-1;end end 1.4全部源代码

s=[0.25,0.25,0.25,0.25];len=size(s,2);if sum(s)~=1 error('s is not a complete set');end h=0;for i=1:1:len h=h-s(i)*log(s(i))/log(2);

结果：h=

2、lc=

2、A=00 01 10 11 1.5 MATLAB仿真及其截图

2、霍夫曼编码

2.1霍夫曼编码源代码

T(t-1)=r;T(t)=0;T=fliplr(sort(T));t=t-1;end B;%输出编码表

END1=sym('[0,1]');%给最后一列的元素编码 END=END1;t=3;d=1;for j=n-2:-1:1%从倒数第二列开始依次对各列元素编码 for i=1:t-2 if i>1 & B(i,j)==B(i-1,j)d=d+1;else d=1;end B(B(n,j+1),j+1)=-1;temp=B(:,j+1);

x=find(temp==B(i,j));END(i)=END1(x(d));

2.2MATLAB仿真及其截图

3、游程编码

3.1、游程编码源代码 I=checkerboard(10,2);[m,n]=size(I);J=[];for i=1:m

压缩图像大小：

Name Size Bytes Class J 1x400 3200 double array 图像压缩比：4 3.2 MATLAB仿真及其截图

就直接向他们求教，虽然有些时候他们也不一定能解答出来，但最少也能给我一些提示。最后有时候也靠自己，慢慢的想，慢慢的翻资料一遍一遍的输入，校正，最后就做出来了。虽然有些题目我也不知道自己的对不对，但是通过自己动手把它做出结果了我就觉得有收获了。

Matlab的课比较少，课堂上学到的东西也不是特别的多，但是通过实验，我自学了许多东西，也完成了学习任务，同时也锻炼了自己的个人能力

2012/5/3

314-