【OSQP】系列(1)——OSQP的安裝與使用

一旺坠、OSQP與osqp-eigen

1.1 OSQP

OSQP是一種二次規(guī)劃求解器乔遮,可用求解線性組合或二次規(guī)劃問(wèn)題,在同類問(wèn)題中求解效率極高取刃!

OSQP官網(wǎng):https://osqp.org/蹋肮,其C語(yǔ)言版本由oxford大學(xué)開(kāi)發(fā)出刷。

其特點(diǎn)如下:

  • Efficient:采用了ADMM-based first-order算法.
  • Robust:只需要問(wèn)題本身為凸即可, 對(duì)問(wèn)題數(shù)據(jù)無(wú)要求.
  • Free:Apache 2.0 license
  • Embaddable:生成嵌入式 C 代碼, 不需要管理內(nèi)存.
  • Interface:多語(yǔ)言,跨平臺(tái),支持 Matlab, Python, R, Julia
官方osqp接口
  • Library-Free:無(wú)需安裝依賴.
  • 同類性能 bechmark 如下圖


    bechmark

1.2 C++接口

OSQP是用純C實(shí)現(xiàn)的,但是純C中對(duì)于矩陣和向量的操作不是特別友好坯辩,因此為了更方便使用OSQP馁龟,在構(gòu)造C++接口時(shí),引入了Eigen庫(kù)漆魔,方便矩陣和向量的操作坷檩。

注意官方支持python和matlab,但是無(wú) C++ 接口改抡,但給出兩個(gè)推薦的第三方社區(qū)維護(hù)的接口

第三方接口

osqp-cpposqp-eigen類似矢炼, 區(qū)別如下:

  • license:osqp-cpp 為 MIT(最寬松),osqp-eigen 為 LGPL(傳染式開(kāi)源)
  • dependency:兩者除都依賴 Eigen 以外, osqp-cpp 依賴 google 自家的 abseil-cpp 庫(kù)阿纤,osqp-eigen 沒(méi)有其他依賴句灌。

推薦使用osqp-eigen作為osqp的C++接口

1.3 osqp-eigen介紹

(1) osqp-eigen軟件組織形式:

軟件組織形式

(2) 相關(guān)類
osqp-eigen主要是構(gòu)造了三個(gè)類,Data類欠拾、Setting類胰锌、Solver類:

  • Data類:對(duì)C語(yǔ)言O(shè)SQP中的OSQPData(struct)進(jìn)行包裝
  • Setting類:對(duì)C語(yǔ)言O(shè)SQP中的OSQPSetting(struct)進(jìn)行包裝
  • Solver類:對(duì)C語(yǔ)言O(shè)SQP中的OSQPWorkspace(struct)進(jìn)行包裝

(3) 求解過(guò)程

  • 1.創(chuàng)建求解器Solver的實(shí)例
  • 2.設(shè)置Setting參數(shù)
  • 3.設(shè)置Data參數(shù)
  • 4.初始化求解器Solver
  • 5.求解優(yōu)化問(wèn)題
  • 6.提取最優(yōu)解

osqp-eigen求解過(guò)程對(duì)應(yīng)代碼如下:

osqp-eigen求解過(guò)程對(duì)應(yīng)代碼

二、osqp-eigen的安裝與使用

本文主要介紹OSQP的C++接口藐窄,即OSQP和osqp-eigen安裝與使用方法匕荸。

2.1 OSQP安裝

OSQP安裝很簡(jiǎn)單,直接把github的庫(kù)拉下來(lái)編譯枷邪,安裝即可。

git clone --recursive https://github.com/osqp/osqp
cd osqp
mkdir build
cd build
cmake .. -DBUILD_SHARED_LIBS=ON
make
sudo make install

注意:如果不將OSQPsudo make install到系統(tǒng)目錄下诺凡,則需要修改osqp-eigen的CMakeLists.txt文件东揣,手動(dòng)找到OSQP庫(kù)的位置,建議執(zhí)行安裝命令腹泌。

2.2 OSQP測(cè)試

OSQP內(nèi)置demo示例嘶卧,可以進(jìn)入build/下進(jìn)行測(cè)試

demo位于~/osqp/build/out/路徑下

直接在終端中運(yùn)行
./osqp_demo

osqp_demo終端輸出結(jié)果

2.3 osqp-eigen安裝

osqp-eigen依賴eigen3,ubuntu系統(tǒng)里內(nèi)置安裝了eigen3包凉袱。

查看是否安裝了eigen3
locate eigen3

osqp-eigen源碼編譯安裝

git clone https://github.com/robotology/osqp-eigen.git
cd osqp-eigen
mkdir build 
cd build
cmake ..    #默認(rèn)安裝在/usr/local/include中,非apt安裝的包芥吟,不在/usr/include/下,而是/usr/local/include下
make
sudo make install

注意:如果不將osqp-eigen sudo make install到系統(tǒng)目錄下专甩,則需要修改自己的CMakeLists.txt文件钟鸵,手動(dòng)找到osqp-eigen庫(kù)的位置,建議執(zhí)行安裝命令涤躲。

2.4 osqp-eigen測(cè)試

osqp-eigen內(nèi)置example示例棺耍,但該示例未與osqp-eigen同時(shí)編譯,需要單獨(dú)編譯种樱,該example位于~/osqp-eigen/example路徑

編譯

cd 
osqp-eigen/example
mkdir build 
cd build
cmake ..
make

編譯完后蒙袍,在終端中運(yùn)行
./MPCExample

MPCExample終端輸出結(jié)果

2.5 CMake工程中調(diào)用osqp-eigen

與調(diào)用其它庫(kù)類似俊卤,需要三個(gè)部分,分別是(1)在頭文件中包含osqp-eigen的頭文件害幅,(2)在cpp文件中實(shí)例化osqp-eigen中對(duì)應(yīng)的類消恍,(3)在CMakeLists.txt文件中添加osqp-eigen庫(kù)路徑。

(1) 頭文件

// eigen庫(kù)
#include <Eigen/Dense>  

// osqp-eigen庫(kù)
#include "OsqpEigen/OsqpEigen.h"

(2) cpp文件


實(shí)例化與求解

(3) CMakeLists.txt文件

cmake_minimum_required(VERSION 3.0)
project(myproject)
find_package(OsqpEigen REQUIRED)
add_executable(example example.cpp)
target_link_libraries(example OsqpEigen::OsqpEigen)
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末以现,一起剝皮案震驚了整個(gè)濱河市狠怨,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌叼风,老刑警劉巖取董,帶你破解...
    沈念sama閱讀 211,265評(píng)論 6 490
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異无宿,居然都是意外死亡茵汰,警方通過(guò)查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 90,078評(píng)論 2 385
  • 文/潘曉璐 我一進(jìn)店門孽鸡,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)蹂午,“玉大人,你說(shuō)我怎么就攤上這事彬碱《剐兀” “怎么了?”我有些...
    開(kāi)封第一講書(shū)人閱讀 156,852評(píng)論 0 347
  • 文/不壞的土叔 我叫張陵巷疼,是天一觀的道長(zhǎng)晚胡。 經(jīng)常有香客問(wèn)我,道長(zhǎng)嚼沿,這世上最難降的妖魔是什么估盘? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 56,408評(píng)論 1 283
  • 正文 為了忘掉前任,我火速辦了婚禮骡尽,結(jié)果婚禮上遣妥,老公的妹妹穿的比我還像新娘。我一直安慰自己攀细,他們只是感情好箫踩,可當(dāng)我...
    茶點(diǎn)故事閱讀 65,445評(píng)論 5 384
  • 文/花漫 我一把揭開(kāi)白布。 她就那樣靜靜地躺著谭贪,像睡著了一般境钟。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上俭识,一...
    開(kāi)封第一講書(shū)人閱讀 49,772評(píng)論 1 290
  • 那天吱韭,我揣著相機(jī)與錄音,去河邊找鬼。 笑死理盆,一個(gè)胖子當(dāng)著我的面吹牛痘煤,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播猿规,決...
    沈念sama閱讀 38,921評(píng)論 3 406
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼衷快,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼!你這毒婦竟也來(lái)了姨俩?” 一聲冷哼從身側(cè)響起蘸拔,我...
    開(kāi)封第一講書(shū)人閱讀 37,688評(píng)論 0 266
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤,失蹤者是張志新(化名)和其女友劉穎环葵,沒(méi)想到半個(gè)月后调窍,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 44,130評(píng)論 1 303
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡张遭,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,467評(píng)論 2 325
  • 正文 我和宋清朗相戀三年邓萨,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片菊卷。...
    茶點(diǎn)故事閱讀 38,617評(píng)論 1 340
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡缔恳,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出洁闰,到底是詐尸還是另有隱情歉甚,我是刑警寧澤,帶...
    沈念sama閱讀 34,276評(píng)論 4 329
  • 正文 年R本政府宣布扑眉,位于F島的核電站纸泄,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏腰素。R本人自食惡果不足惜聘裁,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,882評(píng)論 3 312
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望耸弄。 院中可真熱鬧,春花似錦卓缰、人聲如沸计呈。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 30,740評(píng)論 0 21
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)捌显。三九已至,卻和暖如春总寒,著一層夾襖步出監(jiān)牢的瞬間扶歪,已是汗流浹背。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 31,967評(píng)論 1 265
  • 我被黑心中介騙來(lái)泰國(guó)打工, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留善镰,地道東北人妹萨。 一個(gè)月前我還...
    沈念sama閱讀 46,315評(píng)論 2 360
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像炫欺,于是被迫代替她去往敵國(guó)和親乎完。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 43,486評(píng)論 2 348

推薦閱讀更多精彩內(nèi)容