1. 從“搭積木”到“造積木”為什么我們需要自定義S函數模塊如果你用過Simulink肯定熟悉那些琳瑯滿目的庫模塊增益、積分、傳遞函數、各種信號源和示波器。用這些現成的模塊搭建系統就像玩高級樂高效率很高。但當你需要實現一個標準庫里沒有的算法或者需要與外部硬件、特定數據格式、甚至是C/C代碼深度交互時你就會發現這些“標準積木”不夠用了。比如你想在Simulink里仿真一個自定義的通信調制解調算法或者實現一個復雜的非線性控制器標準模塊要么組合起來非常繁瑣要么根本無法實現。這時你就需要自己“造積木”。在Simulink中自定義模塊的終極武器就是S函數System-Function。它不是一個具體的模塊而是一個編程接口規范。你可以用MATLAB語言、C、C甚至Fortran來編寫S函數描述一個動態系統的輸入、輸出、狀態以及它們隨時間演化的關系。然后Simulink就能像調用內置模塊一樣在仿真的每個時間步調用你的代碼。這相當于你為Simulink擴展了一套專屬的、功能強大的新“積木”。本次我們要“造”的就是一個實現QPSK正交相移鍵控調制與解調的通信鏈路“積木”。QPSK是數字通信中非常經典的一種調制方式它通過改變載波的相位來傳遞信息一個符號可以攜帶2個比特的信息頻譜效率是BPSK的兩倍。在Simulink中雖然有Communications Toolbox提供專業的通信模塊但通過S函數親手實現一遍其價值遠超“仿真出一個波形”。你能徹底理解數據如何在仿真中流動、狀態如何更新、離散系統如何與連續仿真環境交互這是深入掌握Simulink仿真引擎和數字通信原理的絕佳路徑。2. 解剖S函數理解Simulink仿真引擎的“回調機制”在深入代碼之前我們必須先理解S函數是如何與Simulink仿真引擎協同工作的。Simulink的仿真可以看作一個按時間步推進的循環。在每個時間步引擎需要知道系統當前的狀態是什么輸入信號的值是多少基于當前狀態和輸入輸出應該是什么下一個狀態如果有應該更新為什么仿真初始化、結束時要做什么S函數通過一系列預定義好的“回調函數”Callback Methods來回答這些問題。你可以把S函數想象成一個你編寫的“演員”Simulink引擎是“導演”。在仿真的不同階段“導演”會喊出不同的指令調用對應的回調函數“演員”則執行相應的動作。主要的回調函數包括mdlInitializeSizes: 這是“演員”的第一次亮相。它告訴“導演”我這個系統的基本信息有幾個輸入端口幾個輸出端口是連續系統、離散系統還是混合系統有沒有內部狀態離散狀態或連續狀態采樣時間是多少這個函數在仿真一開始只調用一次用于定義系統的“骨架”。mdlInitializeSampleTimes: 緊接著它告訴“導演”我的“節奏”。是繼承上一級的采樣時間是固定的離散采樣時間比如0.001秒還是連續采樣對于我們的QPSK系統通常處理的是離散的數字符號所以會定義一個固定的離散采樣時間。mdlOutputs: 這是仿真循環中最重要的回調之一。“導演”問“給定當前時間、當前輸入和當前狀態你的輸出是什么”這個函數負責計算并輸出當前時間步的端口值。對于無狀態的靜態系統輸出只依賴于當前輸入對于有狀態的動態系統輸出依賴于當前狀態和/或輸入。mdlUpdate: 對于離散狀態系統這是另一個核心回調。“導演”問“在計算完輸出之后基于當前輸入你的下一個狀態在下一個采樣時刻的狀態應該更新為什么”這個函數負責更新離散狀態。注意mdlOutputs使用的是t時刻的狀態而mdlUpdate計算的是tTs時刻下一個采樣點的狀態。這個順序保證了仿真的因果性。mdlDerivatives: 如果是連續系統這個函數用來計算連續狀態的導數dx/dt供Simulink的數值積分器如ode45使用。我們的QPSK是離散系統通常用不到這個函數。mdlTerminate: 仿真結束時調用用于執行一些清理工作如關閉文件、釋放內存等。對于用MATLAB語言編寫的S函數即M語言S函數這些回調函數被寫在一個巨大的switch語句里由一個標志flag來區分。這是M語言S函數的經典結構模板。理解了這個“導演-演員”的交互模型再看代碼就不會覺得是一團亂麻了。注意S函數有Level-1和Level-2兩種版本。Level-1較為古老使用flag切換Level-2面向對象更強大靈活支持多實例、更豐富的端口和數據類型。對于復雜模塊推薦使用Level-2。但為了清晰理解基本原理我們從經典的Level-1 M語言S函數入手。3. QPSK鏈路核心原理與S函數設計映射在動手寫代碼前我們需要把通信原理“翻譯”成S函數能理解的“系統描述”。一個最簡單的QPSK調制解調鏈路可以分解為以下幾個步驟我們將其映射到S函數的設計要素上發射端調制二進制數據源產生隨機的0/1比特流。這可以作為S函數的輸入也可以直接在S函數內部生成。我們設計為外部輸入更通用。串并轉換將連續的比特流兩兩一組分成同相I路和正交Q路兩支路。例如比特對[0 0]映射到(I, Q) (-1, -1)。這屬于數據處理邏輯將在mdlOutputs中實現。脈沖成型可選但重要為了限制信號帶寬需要將離散的I/Q符號通過一個成型濾波器如升余弦滾降濾波器。這會引入記憶性即當前輸出不僅取決于當前符號還取決于之前的符號。這需要S函數內部有**狀態State**來存儲濾波器狀態或之前的符號。狀態更新在mdlUpdate中完成。載波調制將基帶I/Q信號分別與載波cos(2πfct)和-sin(2πfct)相乘后相加得到帶通信號。在仿真中我們通常直接處理復基帶信號即I j*Q載波調制隱含在后續處理中。因此S函數的輸出可以是復基帶信號。信道簡化 5.加性高斯白噪聲AWGN在信號中加入噪聲。這可以在另一個S函數或使用Simulink內置的AWGN Channel模塊實現。為了簡化我們可以將噪聲添加集成到調制S函數中或者作為獨立的鏈路環節。接收端解調6.匹配濾波/采樣接收端通過匹配濾波器通常與成型濾波器相同最大化信噪比并在最佳時刻采樣。這同樣需要狀態濾波器狀態。 7.判決根據采樣得到的I/Q值判斷它離哪個星座點(±1, ±1)最近并映射回比特對。這屬于無狀態的決策邏輯在mdlOutputs中實現。我們的S函數設計決策為了聚焦S函數本身我們將實現一個包含成型濾波的QPSK調制器和一個包含匹配濾波與判決的QPSK解調器作為兩個獨立的S函數模塊。它們通過一個簡單的加法器模擬噪聲連接。這樣我們可以清晰地看到調制器S函數輸入為二進制比特流輸出為復基帶信號I/Q兩路。內部需要狀態來存儲成型濾波器的延遲線。解調器S函數輸入為含噪的復基帶信號輸出為解調后的二進制比特流。內部同樣需要狀態來存儲匹配濾波器的延遲線。成型濾波器我們選擇簡單的矩形脈沖不整形這樣就沒有碼間串擾也無需狀態便于入門。在后續優化中我們再引入升余弦濾波器及其狀態處理。4. 手把手實現QPSK調制器M語言S函數讓我們從最簡單的開始一個無脈沖成型即使用矩形脈沖的QPSK調制器。它沒有記憶性因此不需要內部狀態。第一步創建S函數模板文件在MATLAB中新建一個.m文件命名為qpsk_modulator_simple.m。開頭是標準的函數定義function [sys,x0,str,ts,simStateCompliance] qpsk_modulator_simple(t,x,u,flag, Ts, EbN0_dB) % QPSK調制器S函數 (簡化版矩形脈沖成型) % 輸入參數: % t: 當前仿真時間 % x: 當前狀態向量本例無狀態為空 % u: 輸入向量本例為單路二進制比特流0/1 % flag: Simulink引擎調用的回調標志 % Ts: 符號周期秒通過S函數參數傳入 % EbN0_dB: 信噪比Eb/N0 (dB)用于計算噪聲功率本例暫未使用 % 輸出: % sys: 根據flag不同返回不同的系統信息 % x0: 初始狀態向量 % str: 保留參數始終設為空矩陣 [] % ts: 采樣時間矩陣 [采樣周期, 時間偏移] % simStateCompliance: 保留參數 switch flag case 0 % flag 0: 初始化回調 (mdlInitializeSizes) [sys,x0,str,ts,simStateCompliance]mdlInitializeSizes(Ts); case 3 % flag 3: 輸出計算回調 (mdlOutputs) sysmdlOutputs(t,x,u); case {1, 2, 4, 9} % flag 1,2,4,9: 本例未使用返回空矩陣 sys []; otherwise error([Unhandled flag ,num2str(flag)]); end第二步實現初始化函數mdlInitializeSizes這個函數定義了模塊的“規格說明書”。function [sys,x0,str,ts,simStateCompliance]mdlInitializeSizes(Ts) % 初始化函數 sizes simsizes; % 獲取一個標準的結構體模板 sizes.NumContStates 0; % 連續狀態數量0 sizes.NumDiscStates 0; % 離散狀態數量0 (無記憶) sizes.NumOutputs 2; % 輸出數量2 (I路和Q路我們輸出兩路實數) sizes.NumInputs 1; % 輸入數量1 (單路比特流) sizes.DirFeedthrough 1; % 直接饋通是。因為輸出直接由當前輸入決定不依賴狀態。 sizes.NumSampleTimes 1; % 采樣時間數量1 sys simsizes(sizes); % 將設置應用到sys x0 []; % 初始狀態為空 str []; % 保留 ts [Ts 0]; % 采樣時間[采樣周期, 時間偏移]。偏移為0表示在整數倍Ts時刻采樣。 simStateCompliance UnknownSimState;這里的關鍵是DirFeedthrough設為1。因為我們的輸出調制符號直接由當前輸入的比特對決定不需要等待下一個采樣點更新狀態。采樣時間ts設置為[Ts, 0]意味著這個模塊以固定的Ts秒為周期運行。第三步實現輸出計算函數mdlOutputs這是核心算法所在。注意輸入u是單個比特0或1。但QPSK每次處理2個比特。我們需要一個機制來緩存比特。由于沒有使用狀態這里有一個關鍵技巧我們假設輸入端口的數據速率是符號速率的兩倍即每個Ts時間內輸入2個比特。但實際上S函數在這個例子中每個Ts被調用一次輸入u只有一個值。因此這個簡單版本實際上有缺陷它每次只能處理1個比特無法完成QPSK映射。這引出了我們第一個重要的設計修正。修正設計我們需要狀態來緩存比特即使是無成型濾波QPSK調制本身也需要記憶2個比特才能映射到一個符號。因此我們必須引入離散狀態來存儲上一個比特。讓我們修改設計sizes.NumDiscStates 1;用于存儲上一個比特。并在mdlUpdate中更新它。同時DirFeedthrough應設為1因為當前輸出由當前輸入u和狀態x上一個比特共同決定。修正后的初始化sizes.NumDiscStates 1; % 離散狀態數量1 (用于緩存上一個比特) sizes.DirFeedthrough 1; % 直接饋通是。輸出依賴于當前輸入和狀態。修正后的輸出函數mdlOutputsfunction sysmdlOutputs(t,x,u) % 輸出計算函數 % x: 狀態這里存儲上一個比特 (0/1) % u: 當前輸入比特 persistent bit_buffer; % 使用持久變量作為另一種緩存方式但更推薦用狀態x % 但為了教學清晰我們這里用狀態x作為緩存。 % 將狀態x(1)解釋為上一個比特 prev_bit x(1); current_bit u(1); % QPSK映射規則: (I, Q) % 比特對 [prev_bit, current_bit] - (I, Q) % 00 - (-1, -1), 01 - (-1, 1), 10 - (1, -1), 11 - (1, 1) % 映射關系可以寫成: I 2*prev_bit - 1, Q 2*current_bit - 1 % 因為prev_bit和current_bit是0或12*bit-1得到-1或1。 I_out 2 * prev_bit - 1; Q_out 2 * current_bit - 1; sys [I_out; Q_out]; % 輸出I和Q兩路信號這里我們假設在仿真開始時狀態x(1)被初始化為0需要在x0中設置。這個映射邏輯將連續的比特流轉換為了I/Q符號。第四步實現狀態更新函數mdlUpdate我們需要在flag2時處理狀態更新。在switch語句中添加case 2 % flag 2: 狀態更新回調 (mdlUpdate) sysmdlUpdate(t,x,u);并實現該函數function sysmdlUpdate(t,x,u) % 狀態更新函數用當前輸入比特更新狀態即緩存當前比特供下一個符號周期使用 % 注意這里有一個重要的時序邏輯。 % 在當前時刻t我們使用狀態x(prev_bit)和輸入u(current_bit)計算輸出。 % 然后我們需要為下一個時刻tTs更新狀態。這個新狀態應該是當前的輸入比特u。 % 因為下一個符號周期當前的u將變成“上一個比特”。 sys u; % 將當前輸入比特存入狀態同時修改初始化中的x0x0 [0]; % 初始狀態假設第一個“上一個比特”為0。一個完整的、可工作的簡單QPSK調制器S函數就完成了。它接收單比特流每兩個比特周期注意由于我們每個Ts處理一個比特對但輸入是每個Ts一個比特所以實際符號周期是2*Ts這里需要仔細思考輸出一個QPSK符號。這里的采樣時間Ts實際是比特周期。符號周期是2*Ts。這是一個需要根據系統設計澄清的點。在更正式的設計中我們可能希望輸入端口接收的是并行的兩路比特I和Q支路或者接收一個整數輸入0,1,2,3代表2個比特。這體現了S函數接口設計的靈活性也反映了實際工程中的權衡。5. 引入脈沖成型讓S函數擁有“記憶”矩形脈沖的頻譜太寬在實際通信中不可行。現在我們引入升余弦滾降濾波器。濾波器是一個有記憶的系統當前輸出是當前輸入和過去多個輸入的加權和。這意味著我們的S函數需要更多的狀態來存儲濾波器的延遲線即過去的一系列輸入符號。我們重新設計調制器假設輸入已經是并行的I、Q兩路符號取值±1。這樣輸入端口數為2每個Ts輸入一個復數符號由I和Q兩個實數構成。成型濾波器作用在I路和Q路上它們是獨立的、相同的實濾波器。狀態設計對于一個長度為L階數為L-1的FIR濾波器我們需要存儲過去的L-1個輸入符號。我們可以用一個長度為L-1的狀態向量來實現這個延遲線。每次新的符號到來我們將其插入延遲線頭部并彈出最舊的符號。修改初始化function [sys,x0,str,ts,simStateCompliance]mdlInitializeSizes(Ts, L) % L: 升余弦濾波器的長度抽頭數 sizes simsizes; sizes.NumContStates 0; sizes.NumDiscStates 2*(L-1); % 狀態數量為I路和Q路分別存儲(L-1)個過去符號 sizes.NumOutputs 2; % 輸出濾波后的I和Q sizes.NumInputs 2; % 輸入當前的I_symbol和Q_symbol (±1) sizes.DirFeedthrough 1; % 是因為輸出與當前輸入有關FIR濾波包含當前輸入 sizes.NumSampleTimes 1; sys simsizes(sizes); x0 zeros(2*(L-1), 1); % 初始狀態設為0表示初始時延遲線內全是0 str []; ts [Ts 0]; simStateCompliance UnknownSimState;修改輸出函數mdlOutputs這里需要實現FIR濾波運算。假設我們已經有了濾波器的系數向量b長度為L可以通過S函數參數傳入。function sysmdlOutputs(t,x,u,b) % u(1): 當前I路輸入符號, u(2): 當前Q路輸入符號 % x: 狀態向量前(L-1)個是I路延遲線后(L-1)個是Q路延遲線 % b: FIR濾波器系數向量 (長度為L) L length(b); % 提取I路延遲線狀態 I_delay_line [u(1); x(1:L-1)]; % 當前輸入 過去狀態 % 提取Q路延遲線狀態 Q_delay_line [u(2); x(L:end)]; % 注意狀態向量的排列 % FIR濾波計算 I_filtered b(:) * I_delay_line(:); % 向量點乘 Q_filtered b(:) * Q_delay_line(:); sys [I_filtered; Q_filtered];修改狀態更新函數mdlUpdate更新延遲線為下一個時刻做準備。function sysmdlUpdate(t,x,u,L) % 更新延遲線狀態將當前輸入移入最舊的值移出 % x: [I_delay(1:L-1); Q_delay(1:L-1)] % 更新后: [u(1); I_delay(1:end-1); u(2); Q_delay(1:end-1)] sys x; % 先復制原狀態 % 更新I路延遲線整體向下移位最舊的被丟棄新的u(1)放在最前面 sys(1:L-1) [u(1); x(1:L-2)]; % 更新Q路延遲線 sys(L:end) [u(2); x(L:end-1)];現在這個S函數實現了一個完整的、帶脈沖成型的QPSK調制器。輸入是符號速率的I/Q符號輸出是經過成型濾波的基帶信號。濾波器系數b和長度L可以作為S函數的參數在Simulink模塊對話框中設置。6. 構建完整鏈路在Simulink中集成與調試有了調制器和解調器解調器是類似的包含匹配濾波和判決判決時需要考慮噪聲影響S函數后我們可以在Simulink中搭建完整鏈路。創建S-Function模塊在Simulink庫瀏覽器中找到“User-Defined Functions”庫拖入一個“S-Function”模塊。雙擊它在“S-function name”框中填入你的M文件名如qpsk_modulator_pulse_shaping在“S-function parameters”框中填入參數如Ts, L, b。這些參數會在仿真時傳遞給S函數。配置輸入源使用“Bernoulli Binary Generator”產生隨機比特流。由于我們的調制器輸入需要的是I/Q符號±1可能需要一個“Bit to Integer”或自定義的映射模塊將每2個比特映射為一個0-3的整數再映射為I/Q符號。更直接的方法是使用兩個“Bernoulli Binary Generator”分別產生I路和Q路的比特0/1然后通過2*u-1的運算轉換為±1。這可以用一個MATLAB Function模塊或Gain/Bias模塊組合實現。添加噪聲使用“AWGN Channel”模塊。需要設置Eb/N0、符號周期、每符號能量等參數。對于QPSK每符號能量Es 2*Eb因為一個符號承載2比特。AWGN模塊的輸入信號功率需要根據你的信號幅度計算。如果I/Q符號為±1則符號能量Es I^2 Q^2 2。在加噪時需仔細換算。連接解調器解調器S函數輸入是含噪的I/Q兩路信號內部先進行匹配濾波系數與發射端成型濾波器相同然后采樣在最佳采樣點通常就是符號周期整數倍時刻最后進行判決。判決邏輯是判斷濾波后的I/Q值最接近哪個星座點(±1, ±1)。觀察與驗證使用“Scope”觀察發射端I/Q信號、接收端濾波前后信號、以及比特誤差。使用“Error Rate Calculation”模塊計算誤比特率BER并與理論QPSK的BER公式0.5*erfc(sqrt(Eb/N0))進行比較。在較高信噪比下仿真結果應與理論值接近。使用“Spectrum Analyzer”觀察成型前后信號的頻譜可以看到升余弦濾波器如何有效地限制帶寬。調試S函數的實用技巧使用disp或fprintf在S函數的各個回調函數中插入打印語句輸出時間t、輸入u、狀態x、輸出sys的值。這是最直接的調試方法。注意在MATLAB Function模塊中也可以這樣做。檢查采樣時間確保所有模塊的采樣時間設置正確且兼容。特別是S函數的ts設置。使用Simulink的“Display - Sample Time - Colors”功能用不同顏色顯示不同采樣時間的信號線。驗證狀態更新邏輯這是最容易出錯的地方。確保在mdlOutputs中使用的狀態是“當前時刻”的狀態而在mdlUpdate中計算的是“下一時刻”的狀態。可以手動計算幾個時間步與仿真輸出對比。處理數據類型Simulink信號有數據類型double, single, int8等。確保你的S函數內部計算能處理輸入的數據類型或者使用mdlInitializeSizes中的DWork向量來指定和轉換數據類型。對于M語言S函數通常使用double類型最簡單。7. 從M語言到C MEX S函數追求極致的仿真速度當你用M語言寫完S函數并驗證功能正確后可能會發現仿真速度較慢尤其是當濾波器階數很高、仿真時間很長時。這是因為M語言是解釋執行的。對于計算密集的部分我們可以將其升級為C MEX S函數。C MEX S函數本質是一個用C語言編寫、并被編譯成MEX文件在Windows上是.mexw64文件的S函數。它被Simulink直接以機器碼調用速度遠超M語言。編寫C MEX S函數更復雜需要熟悉C語言和Simulink C S函數的API。關鍵步驟獲取模板MATLAB提供了模板文件。你可以使用命令edit sfuntmpl_basic.c來查看一個基礎的C S函數模板。更簡單的方法是使用Legacy Code Tool傳統代碼工具將已有的C函數封裝成S函數或者使用S-Function BuilderGUI工具。實現回調函數與M語言類似你需要實現mdlInitializeSizes,mdlInitializeSampleTimes,mdlOutputs,mdlUpdate等函數但使用的是C語言語法和Simulink的宏如ssGetInputPortSignal,ssGetOutputPortSignal,ssGetDiscStates。編譯在MATLAB命令窗口使用mex命令編譯你的.c文件例如mex qpsk_modulator.c sfun.c。這需要你的系統配置了支持的C編譯器如Microsoft Visual C。在Simulink中使用和M語言S函數一樣拖入S-Function模塊名稱填寫編譯后的MEX文件名不含擴展名。一個簡單的C MEXmdlOutputs函數示例概念性static void mdlOutputs(SimStruct *S, int_T tid) { // 獲取輸入、輸出、狀態指針 real_T *u ssGetInputPortRealSignal(S, 0); // 輸入指針 real_T *y ssGetOutputPortRealSignal(S, 0); // 輸出指針 real_T *x ssGetDiscStates(S); // 離散狀態指針 // 從S函數參數獲取濾波器系數和長度 real_T *b mxGetPr(ssGetSFcnParam(S, 0)); // 假設第一個參數是系數向量 int_T L mxGetNumberOfElements(ssGetSFcnParam(S, 0)); // 實現FIR濾波計算 (簡化未考慮延遲線索引細節) real_T I_sum 0.0, Q_sum 0.0; for (int i0; iL; i) { I_sum b[i] * get_I_delay_element(x, i, L); // 需要自定義函數獲取延遲線元素 Q_sum b[i] * get_Q_delay_element(x, i, L); } y[0] I_sum; // I路輸出 y[1] Q_sum; // Q路輸出 }編寫C MEX S函數需要細致的內存管理和指針操作但帶來的性能提升在大型仿真中是非常可觀的。對于通信系統仿真將核心算法如濾波、編碼、解碼用C MEX實現是工業級仿真的常見做法。8. 避坑指南與性能優化實戰心得在實現和集成S函數的過程中我踩過不少坑也總結了一些優化經驗1. 采樣時間同步陷阱問題整個鏈路中比特源、調制器、信道、解調器可能有不同的采樣時間。如果設置不當會導致信號維度不匹配或仿真錯誤。對策明確每個模塊的“時間基準”。通常以符號周期為基準最為清晰。比特源的采樣時間應為Ts/2如果每符號2比特調制器和解調器為Ts成型濾波器內部如果采用多倍采樣上采樣則需要更快的采樣時間。使用“Rate Transition”模塊處理不同速率信號間的轉換或者確保S函數能處理多速率輸入通過設置ts為[0, 0]繼承輸入端口速率。2. 狀態初始化與邊界處理問題濾波器在開始階段延遲線是空的初始化為0。這會導致開始的幾個輸出不正確瞬態效應。在計算誤碼率時通常需要忽略開頭的一段數據。對策在S函數中可以添加一個計數器狀態在濾波器充滿之前不進行有效輸出或者直接輸出0。在誤碼率計算模塊中設置一個“延遲”參數跳過開始的若干符號。3. 復數信號處理問題Simulink默認信號是實數的。處理復基帶信號I jQ時需要拆分成兩路實數信號I路和Q路進行傳輸和處理這使模型連線翻倍容易混亂。對策可以考慮使用real-imag到complex的轉換模塊或者在S函數內部直接處理復數。對于C MEX S函數可以創建COMPLEX類型的端口。但在大多數情況下保持I/Q分離更直觀也便于用Scope觀察。4. 參數化與模塊封裝問題S函數參數如Ts,L,b在模塊對話框中以逗號分隔的字符串形式輸入不易管理和修改尤其是濾波器系數向量b。對策使用MATLAB工作空間變量。在模塊參數框中填寫變量名如Ts_symbol,rc_filter_coeff。在仿真腳本或模型回調函數如PreLoadFcn中計算這些變量。更好的方法是使用模塊封裝Mask。右鍵點擊S-Function模塊選擇“Mask - Create Mask”。在封裝編輯器中可以創建友好的參數輸入控件如編輯框、下拉菜單并為每個參數指定一個變量名。在“Initialization”頁簽你可以寫MATLAB代碼來處理這些輸入例如將字符串轉換成數值向量。封裝后的模塊看起來和Simulink內置模塊一樣專業。5. 性能分析使用Simulink Profiler通過simulink.profiler.start和simulink.profiler.stop來分析仿真時間都花在哪里。如果你發現你的S函數占了大部分時間那就是將其轉換為C MEX S函數的最佳信號。向量化運算在M語言S函數的mdlOutputs中避免使用for循環處理向量。盡量使用MATLAB的向量化操作。例如FIR濾波用dot(b, delay_line)或卷積conv注意處理邊界會比循環快很多。通過這個從原理到實現從簡單到復雜從M語言到C語言的QPSK鏈路S函數開發旅程你獲得的不僅僅是一個能跑的仿真模型。你獲得的是對Simulink仿真內核如何調度一個自定義組件的深刻理解是對數字通信基帶處理流程的親手實踐更是將算法思想轉化為可執行、可集成、可優化仿真模塊的工程能力。下次當Simulink庫模塊無法滿足你的需求時你會自信地知道可以自己動手造出最合適的那塊“積木”。