林叶很快就进入了状态。

吃饱了饭，又喝了一杯咖啡提神，此刻精神十分亢奋。

“离散模型中使用代数迭代法（ART）对网格化的离散数据直接建模，用矩阵迭代的方法对像素网格进行处理，实现图像的重构。”

这是林叶的思路。

至于数学模型具体如何，林叶还在思考。

这个离散模型，最重要的还是迭代问题。

FBP（滤波反投影）算法具有的缺陷与限制，只能应用代数迭代才能完美解决。

“这里得先引入投影矩阵的概念...。”

林叶嘀咕说道。

时间一点一点的过，林叶空白的A4纸上现在已经写满了数学公式。

r_ij={1，i号射线通过j号像素中心//0，其他。

记x=[x1，x2，...，xJ]^T为J维图像矢量，x_j表示图像上第j个像素点的吸收强度，

p=[p1，p2...pI]^T为I维投影矢量，p_i表示第i条射线所经过的所有像素的总射线投影;

Rx=p;

...

林叶看着草稿纸上最为核心的迭代模型，

x^(k 1)=x^(k) λ^(k)(q_ik?r_ikx^(k)/||r_ik||^2)

林叶嘀咕道：

“这里还差一点东西，那就是降噪算法。不论何种重建算法，不可避免会带入伪影点。

组委会给的数据含有噪声，成像质量可能并不会达到最好的要求。这个其余学校的建模选手也应该会认识到这个问题才对。”

“这里必须得用算法进行降噪。我得查一查这方面的文献。”

林叶用笔压着笔记本电脑前的草稿纸，然后查阅起了文献。

“这种问题应该在国外医学SCI期刊上比较好寻找合适的算法。”

看了一个多小时，林叶总算找到了一种合适的算法。

“这种算法研究一下能不能改进。”

林叶现在数学水平已经很高了，但是想要在一些顶级SCI期刊上动刀别人的算法模型，

可能性不大。

“看来这里只能引用了，这里用这种算法进行降噪处理已经算是完美答卷了，强行改只会画蛇添足。”

林叶有些丧气，

“早知道这段时间不忙着做第六個任务了。”

就在林叶准备直接照着抄的时候，灵光突然一闪，

到时候辅导员可别迁怒到他们身上。

看着有些精神失常的二人，不少人嘀咕。

“要不小声提醒下他们？”

“算了算了，可能是压力太大了吧，或者想到了思路，我们自己做自己的。”

来自锦大数院的压力，他们好像也十分理解，锦大数院的学科评估毕竟是A，高考收分高得不得了。

机房中别人的低语，林叶等人没有理会。

解释不清楚，等成绩出来再解释。

现在姑且被当做神经病吧。

杨青青已经写完了第一个问题的所有程序，

并且已经运行出了结果。

“小叶子，伱来看看运行的结果符合你的预期不？程序是十分完美的，算法也无可挑剔。

那三种用微元法特殊情况算出来的与我们用最小二乘法拟合的模型结果一模一样，我感觉应该没错。”

林叶已经来到了杨青青身边。

杨青青是用MATLAB写的程序，

林叶也学过一学期的MATLAB，杨青青写的代码林叶也看得懂。

都是基于林叶的数学模型写的。

林叶小声说道：

“程序我这块不如你，但是目前看下来程序十分完美，没有错的地方。”

杨青青十分自恋的说道：

“民大第一程序员，岂是浪得虚名，完美的数学模型与完美的程序，加上最后我们两个一起写论文，高教社杯我觉得稳了。”

林叶说道：

“还有两个问呢，现在先吃饭，吃完了你开始写第二个问题的程序，模型我已经写好了。”

杨青青说道：

“我刚刚看了第三个问，好像跟第二问差不多，只是数据更为复杂了。”

林叶说道：

“我已经有思路了，第三个问需要用到概率论与数理统计的知识，然后反复运用第二个问的数学模型就可以了。

二三问我本来想连着一起做的，不过考虑到论文更完善，更可靠，所以才用概率论与数理统计里面的知识，这可能也是出题人埋下的一个小伏笔。”

杨青青迟疑了一下说道：

“反复运用第二个问模型好像也不是不行，方法与模型我们已经极为完美了，现在缺的就是找到最为精确的十个吸收率。”

林叶说道：

“没错。”