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

华科 生物光学原理与成像 专业选修课 基本都是我自己做的，比往年的版本应该细致很多，学弟学妹需要就下载，觉得可以给点好评。 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 reflectance Temp1: Rd1=0.04+0.222321=0.262321 Temp2: Rd2=0.265571 Temp3: Rd3=0.266965 Temp4: 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) 利用公式：F n_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，根据题意条件，十分简单的在极坐标系中就是画L 5.10 公式为：P （PS, 5.12和5.15做的比较匆忙，出来的数据居然有抖动，后来没有太明白抖动如何出来的，但是趋势应该是对的） 5.12 两次n_rel均取的1.（去掉界面反射带来的区别） 会有两组结果，然后画图 画图得到： 5.15：和5.12区别不大