生物光学原理与成像最后作业.docx

华科 生物光学原理与成像 专业选修课基本都是我自己做的，比往年的版本应该细致很多，学弟学妹需要就下载，觉得可以给点好评。1.7 本来想做一个动态小球行走的过程，后来因为总是之前的轨迹会消失，退而求其次就做静态轨迹显示，技术含量就没有那么高了。下面是循环三次的循环主体：结果：1.9 我采用的是从excel导入Hb、HbO2、波长数据结果完全和书上一样，美化了好久1.10 water用的是1994年的那个数据画图之前先要用如下公式换算成mua：μa (lambda) = (2.303) e (lambda) (150 g/liter)/ (64,500 g Hb/mole) = 0.0054e (lambda)黑圈标出来了那个红外底2.6 （1）计算代码和结果如下（2）调用mie,大概用了6min首先把mie改成function,然后再调用主程序部分代码截图，自认为这个编的不错！下面是结果：和书上非常相似，比某些往年版本好看多了，不足在于g的后半部分的褶皱没有很好显示，实际数据中是有这个抖动的。(3)利用的公式是：σ这个我用的是双对数坐标，画出来是个线性的3.10(1)首先对TEMPLATE.MCI里面的程序改变参数，层数，光子数等，与表格3.2 规定的参数相同下图显示了我主要设置的参数，由于要run 10次求平均值，所以每次的参数都是一样的，注意改输出文件的名称就好了。运行，得到十次结果，然后得出每次的Total Reflectance= Specular reflectance + Diffuse reflectanceTemp1: Rd1=0.04+0.222321=0.262321Temp2:Rd2=0.265571Temp3:Rd3=0.266965Temp4:Rd4=0.261868同理后面的：Rd5=0.253162, Rd6=0.258984, Rd7=0.258818, Rd8=0.25673, Rd9=0.258067, Rd10=0.246791计算出来均值为：Rd=0.258928 R_err=0.005916 和书中TABLE3.2相比误差较小（2）按题目要求改了参数Run后结果：(3)利用公式：Fn_rel=1时的计算结果（部分截图）：(用了671秒，1000000个光子，一百万个光子，这个时间最久)n_rel=1.37时的计算结果（部分截图）：(用了1048秒，1000000个光子，一百万个光子)MATLAB拿数据作图：作图代码：图：标准的和书上一样，哈哈，比往年版本好3.11把原来算法的这个地方改成单独的mua或者mus.把这个地方的dwa改成w*(mua/mus) 或者w*[1-exp(-mua*s)]5.9画公式5.36，根据题意条件，十分简单的在极坐标系中就是画L5.10公式为：P（PS, 5.12和5.15做的比较匆忙，出来的数据居然有抖动，后来没有太明白抖动如何出来的，但是趋势应该是对的）5.12 两次n_rel均取的1.（去掉界面反射带来的区别）会有两组结果，然后画图画图得到：5.15：和5.12区别不大