氣動引數的擬合

2021-10-25 13:51:28 字數 1363 閱讀 4505

氣動引數一般有兩種使用方法,插值和擬合,擬合計算較快,寫程式簡單。matlab有自帶函式regress可以進行多元引數最小二乘擬合,可以擬合出想要的形式,得到每一項的係數。

clc;clear;close;

alpha1=xlsread('sr72資料.xlsx','sheet1','c4:c19');

cd1=xlsread('sr72資料.xlsx','sheet1','d4:d19');

plot(alpha1,cd1)

alpha2=xlsread('sr72資料.xlsx','sheet1','e4:e17');

cd2=xlsread('sr72資料.xlsx','sheet1','f4:f17');

hold on

plot(alpha2,cd2)

alpha3=xlsread('sr72資料.xlsx','sheet1','g4:g17');

cd3=xlsread('sr72資料.xlsx','sheet1','h4:h17');

hold on

plot(alpha3,cd3)

hold off

%% regress函式呼叫

ma=[0.3*ones(1,length(alpha1)) 0.4*ones(1,length(alpha2)) 0.5*ones(1,length(alpha3))]';

alpha=[alpha1;alpha2;alpha3];

cd=[cd1;cd2;cd3];

x=[ones(size(cd)) ma.^2 alpha.^2 ma alpha ma.*alpha];

[b,bint,r,rint,stats]=regress(cd,x);

% cd=b(2)ma^2+b(3)alpha^2+b(4)ma+b(5)alpha+b(6)ma*alpha+b(1);

ma=0:0.1:0.5;

alpha=0:30;

[ma alpha]=meshgrid(ma,alpha);

cd=b(2)*ma.^2+b(3)*alpha.^2+b(4)*ma+b(5)*alpha+b(6)*ma.*alpha+b(1);

mesh(ma,alpha,cd)

hold on

plot3(0.3*ones(1,length(alpha1)),alpha1,cd1,'o')

plot3(0.4*ones(1,length(alpha2)),alpha2,cd2,'*')

plot3(0.5*ones(1,length(alpha3)),alpha3,cd3,'s')

xlabel('ma')

ylabel('alpha')

zlabel('cd')

氣動調節的選擇

氣動調節閥隨著科技技術的進步和生活水平的提高廣泛應用於生活中,一般是化工 電力 食品 冶金等等行業使用的居多,是現代工業企業廣泛使用的工業過程控制儀表之一。氣動調節閥通常由氣動執行機構和調節閥連線安裝除錯後形成的組合儀表。可以看出氣動調節閥不是乙個陌生而遙遠的詞彙,它與我們的生活息息相關,有了氣動調...

zookper的安裝與氣動

zookeeper 是apacahe hadoop 的子專案,是乙個樹型的目錄服務,支援變更推送,適合作為 dubbox 服務的註冊中心,工業強度較高,可用於生產環境 1 準備linux環境 沒有linux環境,使用vm,linux環境使用centeros 6.5 演示使用 2 安裝jdk環境 演示...

python 資料擬合求解方程引數

首先引入三件套和scipy import numpy as np import pandas as pd import matplotlib.pyplot as plt from scipy.optimize import curve fit 拿到實驗資料,通過pandas讀取為dataframe ...