林葉很快就進(jìn)入了狀態(tài)。

吃飽了飯，又喝了一杯咖啡提神，此刻精神十分亢奮。

“離散模型中使用代數(shù)迭代法（ART）對網(wǎng)格化的離散數(shù)據(jù)直接建模，用矩陣迭代的方法對像素網(wǎng)格進(jìn)行處理，實(shí)現(xiàn)圖像的重構(gòu)。”

這是林葉的思路。

至于數(shù)學(xué)模型具體如何，林葉還在思考。

這個離散模型，最重要的還是迭代問題。

FBP（濾波反投影）算法具有的缺陷與限制，只能應(yīng)用代數(shù)迭代才能完美解決。

“這里得先引入投影矩陣的概念...。”

林葉嘀咕說道。

時間一點(diǎn)一點(diǎn)的過，林葉空白的A4紙上現(xiàn)在已經(jīng)寫滿了數(shù)學(xué)公式。

r_ij={1，i號射線通過 j號像素中心//0，其他。

記x=[x1，x2，...，xJ ]^T為J維圖像矢量，x_j表示圖像上第j個像素點(diǎn)的吸收強(qiáng)度，

p=[p1， p2 ... pI ]^T為I維投影矢量，p_i表示第i條射線所經(jīng)過的所有像素的總射線投影;

Rx=p;

...

林葉看著草稿紙上最為核心的迭代模型，

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

林葉嘀咕道：

“這里還差一點(diǎn)東西，那就是降噪算法。不論何種重建算法，不可避免會帶入偽影點(diǎn)。

組委會給的數(shù)據(jù)含有噪聲，成像質(zhì)量可能并不會達(dá)到最好的要求。這個其余學(xué)校的建模選手也應(yīng)該會認(rèn)識到這個問題才對。”

“這里必須得用算法進(jìn)行降噪。我得查一查這方面的文獻(xiàn)。”

林葉用筆壓著筆記本電腦前的草稿紙，然后查閱起了文獻(xiàn)。

“這種問題應(yīng)該在國外醫(yī)學(xué)SCI期刊上比較好尋找合適的算法。”

看了一個多小時，林葉總算找到了一種合適的算法。

“這種算法研究一下能不能改進(jìn)。”

林葉現(xiàn)在數(shù)學(xué)水平已經(jīng)很高了，但是想要在一些頂級SCI期刊上動刀別人的算法模型，

可能性不大。

“看來這里只能引用了，這里用這種算法進(jìn)行降噪處理已經(jīng)算是完美答卷了，強(qiáng)行改只會畫蛇添足。”

林葉有些喪氣，

“早知道這段時間不忙著做第六個任務(wù)了。”

就在林葉準(zhǔn)備直接照著抄的時候，靈光突然一閃，

“咦，不對，共形幾何代數(shù)可以在這里應(yīng)用，完美的符合我自己這篇論文的需求。

安明啊安明，我真的太謝謝您辣，這次我真要拿高教社杯了，到時候全國組委會答辯上，我要當(dāng)著所有的媒體狠狠的夸你。”

林葉把改進(jìn)的模型寫了出來，激動得站了起來，

“nice！除非帝都水木的隊(duì)伍也能改進(jìn)這種算法，不然大概率是我拿高教社杯。”

“林葉，你又發(fā)什么瘋？一天發(fā)瘋幾次？”

楊青青忍不住罵道，

“太丟人了，能不能淡定點(diǎn)，你以前沉穩(wěn)的樣子呢，能不能學(xué)學(xué)我，淡定一點(diǎn)。”

林葉興奮的對著楊青青說道：

“你自己過來看我的數(shù)學(xué)模型，完美至極！”

楊青青數(shù)學(xué)不弱，與林葉探討了一個暑假，水平是突飛猛進(jìn)，

一看林葉的數(shù)學(xué)模型，不由得愣住，思維飛速的運(yùn)轉(zhuǎn)。

特別是看到林葉改進(jìn)的算法，楊青青手也不由得顫抖了起來，

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);

“小葉子，我們...我們...我們可能要發(fā)財了，哈哈哈哈哈哈。高教社杯一人獎勵一萬...我程序也寫得十分完美。”

楊青青笑得波瀾起伏，兩眼都是冒的金光，縱然是穿的寬大T恤，美好的曲線也展現(xiàn)得淋漓盡致。

林葉：...。

“死財迷，剛剛你不是還嫌丟人嗎。”

“是本天才膚淺了。”

楊青青收斂起了笑容，一臉和善的說道，比起之前的態(tài)度，完全是兩個人。

“學(xué)長~，學(xué)姐姐~，到底發(fā)生了什么，人家也想?yún)⑴c進(jìn)來嘛，想要與你們一起高興高興，委屈屈...。”

林思憶一臉委屈的說道。

“小學(xué)妹躺好，看你青青學(xué)姐帶你飛，另外能不能麻煩你去拿一下午飯，好餓...，程序還在跑，我得看著。”

林思憶的活比較輕松，只是寫第一個問題數(shù)學(xué)模型的代碼，其余時間都是在看著二人忙碌，倒倒水，遞一遞東西。

機(jī)房其余人已經(jīng)有些麻了，林葉精神這兩天不對，你楊青青怎么也跟著不對啊。

你們別嚇人啊，趙院長、輔導(dǎo)員整個數(shù)院可是把你們兩個寶貝得不得了。

到時候輔導(dǎo)員可別遷怒到他們身上。

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

“要不小聲提醒下他們？”

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

來自錦大數(shù)院的壓力，他們好像也十分理解，錦大數(shù)院的學(xué)科評估畢竟是A，高考收分高得不得了。

機(jī)房中別人的低語，林葉等人沒有理會。

解釋不清楚，等成績出來再解釋。

現(xiàn)在姑且被當(dāng)做神經(jīng)病吧。

楊青青已經(jīng)寫完了第一個問題的所有程序，

并且已經(jīng)運(yùn)行出了結(jié)果。

“小葉子，你來看看運(yùn)行的結(jié)果符合你的預(yù)期不？程序是十分完美的，算法也無可挑剔。

那三種用微元法特殊情況算出來的與我們用最小二乘法擬合的模型結(jié)果一模一樣，我感覺應(yīng)該沒錯。”

林葉已經(jīng)來到了楊青青身邊。

楊青青是用MATLAB寫的程序，

林葉也學(xué)過一學(xué)期的MATLAB，楊青青寫的代碼林葉也看得懂。

都是基于林葉的數(shù)學(xué)模型寫的。

林葉小聲說道：

“程序我這塊不如你，但是目前看下來程序十分完美，沒有錯的地方。”

楊青青十分自戀的說道：

“民大第一程序員，豈是浪得虛名，完美的數(shù)學(xué)模型與完美的程序，加上最后我們兩個一起寫論文，高教社杯我覺得穩(wěn)了。”

林葉說道：

“還有兩個問呢，現(xiàn)在先吃飯，吃完了你開始寫第二個問題的程序，模型我已經(jīng)寫好了。”

楊青青說道：

“我剛剛看了第三個問，好像跟第二問差不多，只是數(shù)據(jù)更為復(fù)雜了。”

林葉說道：

“我已經(jīng)有思路了，第三個問需要用到概率論與數(shù)理統(tǒng)計的知識，然后反復(fù)運(yùn)用第二個問的數(shù)學(xué)模型就可以了。

二三問我本來想連著一起做的，不過考慮到論文更完善，更可靠，所以才用概率論與數(shù)理統(tǒng)計里面的知識，這可能也是出題人埋下的一個小伏筆。”

楊青青遲疑了一下說道：

“反復(fù)運(yùn)用第二個問模型好像也不是不行，方法與模型我們已經(jīng)極為完美了，現(xiàn)在缺的就是找到最為精確的十個吸收率。”

林葉說道：

“沒錯。”