摘要：多項(xiàng)式回歸（Polynomial Regression）是一種回歸分析方法，通過(guò)擬合一個(gè)多項(xiàng)式方程來(lái)模擬自變量與因變量之間的非線性關(guān)系。多項(xiàng)式回歸的目標(biāo)是找到一組多項(xiàng)式系數(shù)，使得擬合曲線盡可能地接近數(shù)據(jù)點(diǎn)。這種方法可以用于預(yù)測(cè)、設(shè)計(jì)和推理、其他數(shù)據(jù)分析任務(wù)。多項(xiàng)式回歸是一種統(tǒng)計(jì)學(xué)和機(jī)器學(xué)習(xí)領(lǐng)域常用的回歸分析方法，屬于線性回歸的一種形式，將自變量x和因變量y之間的關(guān)系建模為n次多項(xiàng)式。

目前求解多項(xiàng)式回歸問(wèn)題，主要是把多項(xiàng)式回歸轉(zhuǎn)化成線性回歸并用最小二乘法得到損失函數(shù)，令偏導(dǎo)等于0來(lái)求解。另外也可以直接使用梯度下降等方法來(lái)最小化損失函數(shù)。這些優(yōu)化技術(shù)的時(shí)間復(fù)雜度是多項(xiàng)式級(jí)別的，隨著問(wèn)題規(guī)模的增長(zhǎng)，很難在短時(shí)間找到滿意的解。

鑒于摩爾定律即將走到盡頭，就要探索非常規(guī)計(jì)算范式（如量子計(jì)算）解決大規(guī)模優(yōu)化問(wèn)題的適用性。量子計(jì)算機(jī)擅長(zhǎng)解決困難的優(yōu)化問(wèn)題，并為加速求解問(wèn)題提供了更快、更精準(zhǔn)有效的替代方案。

在場(chǎng)景應(yīng)用上，多項(xiàng)式回歸分析可用來(lái)估計(jì)功率失配問(wèn)題。如孤島配電系統(tǒng)減載的策略用例上，精準(zhǔn)制定負(fù)載優(yōu)先級(jí)以準(zhǔn)確估計(jì)頻率失衡問(wèn)題，穩(wěn)定系統(tǒng)頻率，高效克服配電系統(tǒng)因分布式發(fā)電（DG）集成而帶來(lái)的系統(tǒng)穩(wěn)定性不足造成的停電等問(wèn)題。此外，在煉鋼場(chǎng)景下，如將生鐵和回收廢料的混合物還原為低碳鋼，對(duì)這一過(guò)程的預(yù)測(cè)及演變能力操作員可以通過(guò)過(guò)程目標(biāo)預(yù)測(cè)系統(tǒng)來(lái)增強(qiáng)。而多項(xiàng)式分析方法就可以從過(guò)程參數(shù)的較小子集來(lái)預(yù)測(cè)過(guò)程，對(duì)復(fù)雜預(yù)測(cè)模型實(shí)現(xiàn)高精確推理。

同樣，在預(yù)測(cè)農(nóng)作物產(chǎn)量和降雨量的關(guān)系、每日用水量預(yù)測(cè)、光伏發(fā)電預(yù)測(cè)等場(chǎng)景下，多項(xiàng)式回歸分析都能高度適配。

5月16日，北京玻色量子科技有限公司（以下簡(jiǎn)稱“玻色量子”）在新品發(fā)布會(huì)上推出的100量子比特相干光量子計(jì)算機(jī)真機(jī)——“天工量子大腦”，旨在快速、高效地求解NP-hard的Ising問(wèn)題。而多項(xiàng)式回歸系數(shù)的擬合就可以轉(zhuǎn)化為一個(gè)Ising/QUBO模型，由“天工量子大腦”快速求解?！疤旃ち孔哟竽X”可以極大簡(jiǎn)化求解步驟并在毫秒級(jí)的時(shí)間里給出較好的擬合系數(shù)取值。同時(shí)，多項(xiàng)式回歸是機(jī)器學(xué)習(xí)領(lǐng)域的基礎(chǔ)問(wèn)題，量子計(jì)算可以在統(tǒng)計(jì)學(xué)和機(jī)器學(xué)習(xí)領(lǐng)域的其他問(wèn)題發(fā)揮更好的效果。

建模思路

首先給出符號(hào)的定義：

下面是多項(xiàng)式回歸的原問(wèn)題，其中d為多項(xiàng)式的次數(shù)

為了將多項(xiàng)式回歸問(wèn)題轉(zhuǎn)換為QUBO形式，我們先將（1）式進(jìn)行改寫(xiě)

接著，引入一個(gè)精度向量P=[p1，p2，...，pk]T，k∈N。p中的每個(gè)項(xiàng)可以是2的整數(shù)次冪，冪指數(shù)可以為正也可以為負(fù)。同時(shí)，引入k個(gè)二元變量。

則每個(gè)系數(shù)w可以表示為：

式（3）可以用矩陣表示如下

式（4）帶入式（2）替換后可得

該問(wèn)題本身沒(méi)有額外約束，我們已經(jīng)將問(wèn)題的變量進(jìn)行二值化處理，式（5）即為多項(xiàng)式回歸的QUBO模型。

案例求解

我們以如下三次曲線為例進(jìn)行多項(xiàng)式回歸分析，

通過(guò)在目標(biāo)值y中加入符合正態(tài)分布的隨機(jī)噪聲, 生成的隨機(jī)樣本圖像如圖1所示。

圖1三次曲線隨機(jī)樣本

在本案例中，回歸系數(shù)w=[w0，w1，w2]，即

我們可以設(shè)置P=[21,20,2-1,...,2-4]，每個(gè)系數(shù)設(shè)置6個(gè)比特，共設(shè)置18個(gè)比特，求解后可以得到

所以，各個(gè)系數(shù)取值結(jié)果如下：

擬合結(jié)果如圖2所示。

圖2擬合結(jié)果

問(wèn)題拓展

多項(xiàng)式回歸是機(jī)器學(xué)習(xí)等領(lǐng)域的基礎(chǔ)的問(wèn)題，量子計(jì)算還可以在機(jī)器學(xué)習(xí)的其他問(wèn)題中發(fā)揮優(yōu)勢(shì)。

未來(lái)，玻色量子將依托100計(jì)算量子比特相干光量子計(jì)算機(jī)真機(jī)——“天工量子大腦”，聚焦“實(shí)用化量子計(jì)算”，不斷深入研究包括NP-Complete問(wèn)題在內(nèi)的眾多問(wèn)題，拓展更多可實(shí)用化量子計(jì)算的真實(shí)應(yīng)用場(chǎng)景。

玻色量子還將啟動(dòng)“燎原計(jì)劃”開(kāi)發(fā)者平臺(tái)，并持續(xù)對(duì)外開(kāi)放“天工量子大腦”的真機(jī)測(cè)試，熱忱歡迎更多不同領(lǐng)域的研究伙伴前來(lái)了解相干量子計(jì)算的原理和能力，在此基礎(chǔ)上展開(kāi)共同研發(fā)，用量子計(jì)算去解決更多真實(shí)場(chǎng)景中的問(wèn)題，讓量子計(jì)算的超強(qiáng)算力能真正服務(wù)于各行各業(yè)，滿足未來(lái)時(shí)代對(duì)于計(jì)算的需求。