林叶很快就进入了状态。
吃饱了饭,又喝了一杯咖啡提神,此刻精神十分亢奋。
“离散模型中使用代数迭代法(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)r_ik;
林叶嘀咕道:
“这里还差一点东西,那就是降噪算法。不论何种重建算法,不可避免会带入伪影点。
组委会给的数据含有噪声,成像质量可能并不会达到最好的要求。这个其余学校的建模选手也应该会认识到这个问题才对。”
“这里必须得用算法进行降噪。我得查一查这方面的文献。”
林叶用笔压着笔记本电脑前的草稿纸,然后查阅起了文献。
“这种问题应该在国外医学SCI期刊上比较好寻找合适的算法。”
看了一个多小时,林叶总算找到了一种合适的算法。
“这种算法研究一下能不能改进。”
林叶现在数学水平已经很高了,但是想要在一些顶级SCI期刊上动刀别人的算法模型,
可能性不大。
“看来这里只能引用了,这里用这种算法进行降噪处理已经算是完美答卷了,强行改只会画蛇添足。”
林叶有些丧气,
“早知道这段时间不忙着做第六個任务了。”
就在林叶准备直接照着抄的时候,灵光突然一闪,
“咦,不对,共形几何代数可以在这里应用,完美的符合我自己这篇论文的需求。
安明啊安明,我真的太谢谢您辣,这次我真要拿高教社杯了,到时候全国组委会答辩上,我要当着所有的媒体狠狠的夸你。”
林叶把改进的模型写了出来,激动得站了起来,
“nice!除非帝都水木的队伍也能改进这种算法,不然大概率是我拿高教社杯。”
“林叶,你又发什么疯?一天发疯几次?”
杨青青忍不住骂道,
“太丢人了,能不能淡定点,你以前沉稳的样子呢,能不能学学我,淡定一点。”
林叶兴奋的对着杨青青说道:
“你自己过来看我的数学模型,完美至极!”
杨青青数学不弱,与林叶探讨了一个暑假,水平是突飞猛进,
一看林叶的数学模型,不由得愣住,思维飞速的运转。
特别是看到林叶改进的算法,杨青青手也不由得颤抖了起来,
NLM(x_i)=∑_(j∈N_i)w(i,j)x_j;
w(i,j)=e^(?||xV_i?xV_j||^2/h^2)/∑_(j∈N_i)e^(?||xV_i?xV_j||^2/h^2);
“小叶子,我们...我们...我们可能要发财了,哈哈哈哈哈哈。高教社杯一人奖励一万...我程序也写得十分完美。”
杨青青笑得波澜起伏,两眼都是冒的金光,纵然是穿的宽大T恤,美好的曲线也展现得淋漓尽致。
林叶:...。
“死财迷,刚刚你不是还嫌丢人吗。”
“是本天才肤浅了。”
杨青青收敛起了笑容,一脸和善的说道,比起之前的态度,完全是两个人。
“学长~,学姐姐~,到底发生了什么,人家也想参与进来嘛,想要与你们一起高兴高兴,委屈屈...。”
林思忆一脸委屈的说道。
“小学妹躺好,看你青青学姐带你飞,另外能不能麻烦你去拿一下午饭,好饿...,程序还在跑,我得看着。”
林思忆的活比较轻松,只是写第一个问题数学模型的代码,其余时间都是在看着二人忙碌,倒倒水,递一递东西。
机房其余人已经有些麻了,林叶精神这两天不对,你杨青青怎么也跟着不对啊。
你们别吓人啊,赵院长、辅导员整个数院可是把你们两个宝贝得不得了。
到时候辅导员可别迁怒到他们身上。
看着有些精神失常的二人,不少人嘀咕。
“要不小声提醒下他们?”
“算了算了,可能是压力太大了吧,或者想到了思路,我们自己做自己的。”
来自锦大数院的压力,他们好像也十分理解,锦大数院的学科评估毕竟是A,高考收分高得不得了。
机房中别人的低语,林叶等人没有理会。
解释不清楚,等成绩出来再解释。
现在姑且被当做神经病吧。
杨青青已经写完了第一个问题的所有程序,
并且已经运行出了结果。
“小叶子,伱来看看运行的结果符合你的预期不?程序是十分完美的,算法也无可挑剔。
那三种用微元法特殊情况算出来的与我们用最小二乘法拟合的模型结果一模一样,我感觉应该没错。”
林叶已经来到了杨青青身边。
杨青青是用MATLAB写的程序,
林叶也学过一学期的MATLAB,杨青青写的代码林叶也看得懂。
都是基于林叶的数学模型写的。
林叶小声说道:
“程序我这块不如你,但是目前看下来程序十分完美,没有错的地方。”
杨青青十分自恋的说道:
“民大第一程序员,岂是浪得虚名,完美的数学模型与完美的程序,加上最后我们两个一起写论文,高教社杯我觉得稳了。”
林叶说道:
“还有两个问呢,现在先吃饭,吃完了你开始写第二个问题的程序,模型我已经写好了。”
杨青青说道:
“我刚刚看了第三个问,好像跟第二问差不多,只是数据更为复杂了。”
林叶说道:
“我已经有思路了,第三个问需要用到概率论与数理统计的知识,然后反复运用第二个问的数学模型就可以了。
二三问我本来想连着一起做的,不过考虑到论文更完善,更可靠,所以才用概率论与数理统计里面的知识,这可能也是出题人埋下的一个小伏笔。”
杨青青迟疑了一下说道:
“反复运用第二个问模型好像也不是不行,方法与模型我们已经极为完美了,现在缺的就是找到最为精确的十个吸收率。”
林叶说道:
“没错。”