原创 matlab PID程序

%設一被控對象G(s)=50/(0.125s^2+7s) %用增量式PID控制算法編寫仿真程序 %輸入分別爲單位階躍、正弦信號,採樣時間爲1ms clear; close all; ts=0.001; %

原创 從文件名中刪除下劃線

在網上下載的英文文獻,經常在文件名中有很多下劃線,每個單詞之間都是用下劃線隔開的,看着不舒服。 想把下劃線全部替換爲空格。用下面的腳本實現: #%% delete_.py重命名某個文件夾下的文件 # 從文件名中刪除下劃線 import

原创 動態條形圖bar chart race

學習一下網上很火的動態條形圖。 學習自:https://pypi.org/project/bar-chart-race/ https://github.com/dexplo/bar_chart_race https://www.cnblog

原创 運籌學與最優化手寫筆記

第一章運籌學思想與運籌學建模   第二章基本概念和基本理論 第三章線性規劃 第四章搜索算法 第五章無約束最優化 第六章約束最優化

原创 matlab既約梯度法習題

既約梯度法的約束條件爲線性,目標函數可爲非線性。 驗證: oldFmt = format("rat"); f=@(x) x(1)^2+x(2)^2-x(1)*x(2)-2*x(1)-3*x(2); A = [1,1 1,5];

原创 matlab求單峯函數最小值

測試函數是-x*sin(x),求其在0~4.5上的最小值。 首先繪圖觀察一下: fplot(@(x) -x*sin(x),[0,4.5]); hold on; fplot(@(x) -sin(x)-x*cos(x),[0,4.5]); h

原创 matlab近似泛函優化

泛函優化可以近似轉化爲各點座標的優化。 如以下題目 用理論方式(歐拉方程)求解比較繁瑣,門檻高,適應性差。當邊界條件或目標函數發生變化時,需要新的理論。 該題目本質是求從(1,1)到(2,2)的最優路徑曲線,使目標函數最小。因此可以轉換爲

原创 matlab簡單線性規劃&單純形法

1. f=[-3,-1,0,0]; A=[2,-1,0,0]; b=[12]; Aeq=[3,3,1,0 4,-4,0,1]; beq=[30,16]; lb=[0,0,0,0]; ub=[]; [x,y] = linprog(

原创 如何求矩陣的逆矩陣

求逆矩陣最有效的方法是初等變換法(雖然還有別的方法)。如果要求方陣 \(A\) 的逆矩陣,標準的做法是: 將矩陣 \(A\) 與單位矩陣 \(I\) 排成一個新的矩陣 \((A \quad I)\) 將此新矩陣 \(( A \qua

原创 Maxwell轉矩對角度直觀理解

之前在這篇https://www.cnblogs.com/dingdangsunny/p/14327019.html#_label0中提到了轉矩隨定轉子相對偏移角度的變化關係,但是其實這該模型中,轉子只有鐵心,即定子對轉子只有吸引力,沒有排

原创 過程選擇模型

前段時間看到趙玉平老師講的關於相親選擇的問題,感覺比較有趣,希望通過概率模擬驗證一下該方法的有效性。 原鏈接如下,感興趣可先了解原講解:管理學博士是怎麼硬核相親的,過程太真實了,最後居然選到這麼好的“對象”!強烈建議大家學學…… %相親  

原创 matlab kriging模型

kriging工具箱:https://orbit.dtu.dk/en/publications/dace-a-matlab-kriging-toolbox x=rand(1,100)*5; y=rand(1,100)*5; z=x./

原创 matlab等高線圖

利用三列數據,創建矩陣,以繪製等高線圖。只適用於對齊的數據,否則需要插值。 函數: function colplot(data,nums,p,leg) %data:待繪圖的數據,每列一個變量 %nums:繪圖所用數據所在列1×3,依次爲顏

原创 matlab多變量繪圖函數(類似ggplot2)

1.多變量繪圖 由於在matlab繪圖時遇到多個變量,需要逐一繪製(還沒找到別的好方法),在此過程中使用hold on保持圖形,比較繁瑣,所以將其封裝成爲以下函數,簡化該過程。目的是以更接近統計學的思維去繪圖,當然還十分粗糙,和ggplot

原创 關於maxwell中電機功率流的梳理

軟件爲Ansys Maxwell,所用模型爲示例庫中的ipm_1,測試方法爲: 打開示例模型; 運行RMxprt仿真; 生成Maxwell2D模型; 運行2D模型(運行前可修改仿真時長和步長); 繪出希望查看的變量並分析。 1.功率梳理