




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、1第第6章章 醫(yī)用薄膜滲透率的確定醫(yī)用薄膜滲透率的確定曲線擬合曲線擬合 2概述概述p利用數(shù)據(jù)確定方程的參數(shù)利用數(shù)據(jù)確定方程的參數(shù)n如:剎車距離與剎車速度的有如下關(guān)系如:剎車距離與剎車速度的有如下關(guān)系 d=Cv2 其中,其中,d表示剎車距離,表示剎車距離,v表示剎車時(shí)的速度,需要確定表示剎車時(shí)的速度,需要確定常數(shù)常數(shù)C。p根據(jù)一組二維數(shù)據(jù),即平面上的若干點(diǎn),要求確定根據(jù)一組二維數(shù)據(jù),即平面上的若干點(diǎn),要求確定一個(gè)一元函數(shù)一個(gè)一元函數(shù)y=f(x),使這些點(diǎn)與曲線總體來說,使這些點(diǎn)與曲線總體來說盡量接近,這就是數(shù)據(jù)擬合曲線的基本思想,簡(jiǎn)稱盡量接近,這就是數(shù)據(jù)擬合曲線的基本思想,簡(jiǎn)稱曲線擬合曲線擬合(
2、fitting a curve)。 36.1 醫(yī)用薄膜的滲透率醫(yī)用薄膜的滲透率測(cè)定薄膜滲透率測(cè)定薄膜滲透率VAVBS用面積為S的薄膜將容器分成體積分別為VA,VB的兩部分,在兩部分中分別注滿兩種不同濃度的溶液,此時(shí)物質(zhì)分子會(huì)從高濃度溶液穿過薄膜向低濃度溶液擴(kuò)散,通過單位面積薄膜分子擴(kuò)散的速度與薄膜兩側(cè)溶液的濃度差成正比,試確定薄膜滲透率K。 46.2 確定醫(yī)用薄膜滲透率的數(shù)學(xué)模型確定醫(yī)用薄膜滲透率的數(shù)學(xué)模型p問題假設(shè)問題假設(shè)n薄膜兩側(cè)的溶液始終是均勻的,即在任何時(shí)刻薄膜兩側(cè)薄膜兩側(cè)的溶液始終是均勻的,即在任何時(shí)刻薄膜兩側(cè)的每一處溶液的濃度都是相同的的每一處溶液的濃度都是相同的n當(dāng)薄膜兩側(cè)的溶液
3、濃度不一致時(shí),物質(zhì)的分子穿透薄膜當(dāng)薄膜兩側(cè)的溶液濃度不一致時(shí),物質(zhì)的分子穿透薄膜總是從高濃度溶液向低濃度溶液擴(kuò)散總是從高濃度溶液向低濃度溶液擴(kuò)散n通過單位面積薄膜分子擴(kuò)散的速度與薄膜兩側(cè)溶液的濃通過單位面積薄膜分子擴(kuò)散的速度與薄膜兩側(cè)溶液的濃度差成正比。度差成正比。n薄膜是雙向同性的,即物質(zhì)從薄膜的任何的一側(cè)向另一薄膜是雙向同性的,即物質(zhì)從薄膜的任何的一側(cè)向另一側(cè)滲透的性能是相同的側(cè)滲透的性能是相同的 5p符號(hào)說明符號(hào)說明nCA(t),CB(t)表示表示t時(shí)刻薄膜兩側(cè)溶液的濃度時(shí)刻薄膜兩側(cè)溶液的濃度nA 、B表示初始時(shí)刻兩側(cè)溶液的濃度(表示初始時(shí)刻兩側(cè)溶液的濃度(mg/cm3)nK表示薄膜滲透
4、率表示薄膜滲透率nVA ,VB表示由薄膜阻隔的容器兩側(cè)的體積表示由薄膜阻隔的容器兩側(cè)的體積6.2 確定醫(yī)用薄膜滲透率的數(shù)學(xué)模型確定醫(yī)用薄膜滲透率的數(shù)學(xué)模型 66.2 確定醫(yī)用薄膜滲透率的數(shù)學(xué)模型確定醫(yī)用薄膜滲透率的數(shù)學(xué)模型問題分析問題分析解決問題的思路:首先通過機(jī)理分析尋找某一側(cè)濃度隨時(shí)間變化的函數(shù)關(guān)系CA(t)或CB(t),其中可能含有待定參數(shù),如薄膜滲透率K,然后根據(jù)一組測(cè)量值(ti,Ci),i=1,2,n去確定模型中待定參數(shù)。t,t+t,t+tt時(shí)段容器時(shí)段容器A A側(cè)的物質(zhì)質(zhì)量增加量側(cè)的物質(zhì)質(zhì)量增加量從從B側(cè)滲透至側(cè)滲透至A側(cè)的質(zhì)量為側(cè)的質(zhì)量為 76.2 確定醫(yī)用薄膜滲透率的數(shù)學(xué)模型確
5、定醫(yī)用薄膜滲透率的數(shù)學(xué)模型由質(zhì)量守恒定律,有由質(zhì)量守恒定律,有=兩邊除以兩邊除以t,令,令t 0并整理得并整理得整個(gè)容器的溶液中物質(zhì)的質(zhì)量不變,有整個(gè)容器的溶液中物質(zhì)的質(zhì)量不變,有(6.1) 8兩邊除以兩邊除以V VA A,并整理得,并整理得6.2 確定醫(yī)用薄膜滲透率的數(shù)學(xué)模型確定醫(yī)用薄膜滲透率的數(shù)學(xué)模型代入代入(6.1)得得再利用初始條件再利用初始條件解出解出 96.2 確定醫(yī)用薄膜滲透率的數(shù)學(xué)模型確定醫(yī)用薄膜滲透率的數(shù)學(xué)模型數(shù)學(xué)模型數(shù)學(xué)模型問題歸結(jié)為利用CB在時(shí)刻tj的測(cè)量數(shù)據(jù)Cj(j=1,2,N)來辨識(shí)參數(shù)K和A,B,對(duì)應(yīng)的數(shù)學(xué)模型變?yōu)槭购瘮?shù)達(dá)到最小令問題轉(zhuǎn)化為如下最優(yōu)化問題 10最小二
6、乘法最小二乘法給定平面上的點(diǎn)給定平面上的點(diǎn)(xi,yi),i=1,2,n,求,求f(x),使,使達(dá)到最小。達(dá)到最小。其中i為點(diǎn)(xi,yi)與曲線y=f(x)的距離。曲線擬合的實(shí)際含義是尋求一個(gè)函數(shù)y=f(x),使f(x)在某種準(zhǔn)則下與所有數(shù)據(jù)點(diǎn)最為接近,即曲線擬合得最好。6.3 一元最小二乘法簡(jiǎn)介一元最小二乘法簡(jiǎn)介 116.3 一元最小二乘法簡(jiǎn)介一元最小二乘法簡(jiǎn)介設(shè)擬合函數(shù)可由一些簡(jiǎn)單的“基函數(shù)”(如冪函數(shù),三角函數(shù)等) 來線性表示:現(xiàn)在要確定系數(shù)c0,c1,cm,使達(dá)到最小 126.3 一元最小二乘法簡(jiǎn)介一元最小二乘法簡(jiǎn)介曲線擬合示意圖 136.3 一元最小二乘法簡(jiǎn)介一元最小二乘法簡(jiǎn)介已知
7、一組數(shù)據(jù)已知一組數(shù)據(jù)(xi,yi),i=1,2,n,如何選擇函數(shù),如何選擇函數(shù)f(x)?一是根據(jù)機(jī)理分析來確定函數(shù)形式一是根據(jù)機(jī)理分析來確定函數(shù)形式二是根據(jù)散點(diǎn)圖直觀判斷函數(shù)二是根據(jù)散點(diǎn)圖直觀判斷函數(shù)f(x)的形式,常用的一的形式,常用的一元曲線擬合函數(shù)有元曲線擬合函數(shù)有雙曲線雙曲線指數(shù)曲線指數(shù)曲線多項(xiàng)式多項(xiàng)式f(x)=a+b/xf(x)=aebx 146.4 用曲線擬合方法確定醫(yī)用薄膜滲透率用曲線擬合方法確定醫(yī)用薄膜滲透率醫(yī)用薄膜滲透率最小二乘法辨識(shí)模型假設(shè):VA=VB=1000cm3,S=10cm2,B部分溶液濃度測(cè)試結(jié)果如下表:(6.2) 156.4 用曲線擬合方法確定醫(yī)用薄膜滲透率用曲
8、線擬合方法確定醫(yī)用薄膜滲透率101202. 0),(minjjKtCbeabaKEj將6.2簡(jiǎn)化為1)編寫M文件(nongdu.m)function f=nongdu(x,tdata)f=x(1)+x(2)*exp(-0.02*x(3)*tdata);其中x(1)=a x(2)=b x(3)=K 166.4 用曲線擬合方法確定醫(yī)用薄膜滲透率用曲線擬合方法確定醫(yī)用薄膜滲透率2)編寫程序(nihe2.m)tdata=linspace(100,1000,10);cdata=4.54,4.99,5.35,5.65,5.9,6.1,6.26,6.39,6.5,6.59*10(-3);x0=0.2,0.0
9、5,0.05;x,resnorm,residual=lsqcurvefit(nongdu,x0,tdata,cdata)t=linspace(100,1000,100);c=nongdu(x,t);plot(tdata,cdata,o,t,c)c1=nongdu(x,tdata);e=c1-cdata; e1=sum(e.*e)輸出結(jié)果:x =0.0063 -0.0034 0.2542e1=3.5604e-0076.4 用曲線擬合方法確定醫(yī)用薄膜滲透率用曲線擬合方法確定醫(yī)用薄膜滲透率曲線擬合效果圖 186.4 用曲線擬合方法確定醫(yī)用薄膜滲透率用曲線擬合方法確定醫(yī)用薄膜滲透率多項(xiàng)式擬合語句a=p
10、olyfit(xdata,ydata,n)其中n表示多項(xiàng)式的最高階數(shù),xdata,ydata為將要擬合的數(shù)據(jù),輸出參數(shù)a為多項(xiàng)式的系數(shù)計(jì)算多項(xiàng)式的擬合值 y=polyval(a,x) 196.4 用曲線擬合方法確定醫(yī)用薄膜滲透率用曲線擬合方法確定醫(yī)用薄膜滲透率汽車剎車距離與剎車速度的觀測(cè)值根據(jù)所給的數(shù)據(jù)畫散點(diǎn)圖plot(v,d,ro) 206.4 用曲線擬合方法確定醫(yī)用薄膜滲透率用曲線擬合方法確定醫(yī)用薄膜滲透率實(shí)測(cè)數(shù)據(jù)d,v的散點(diǎn)圖 216.4 用曲線擬合方法確定醫(yī)用薄膜滲透率用曲線擬合方法確定醫(yī)用薄膜滲透率設(shè)距離d與速度v之間的函數(shù)結(jié)構(gòu)為二次多項(xiàng)式:d=a2v2+a1v+a0,求參數(shù)a=a2
11、,a1,a0的matlab程序如下v=20 25 30 35 40 45 50 55 60 65 70 75 80;d=42 56 73.5 91.5 116 142.5 173 209.5 248 292.5 343 401 464;a=polyfit(v,d,2)dl=polyval(a,v);plot(v,d,ro,v,dl,b.);輸出結(jié)果:a=0.0886 -1.9701 50.0594擬合的多項(xiàng)式為 d=0.0886v2-1.9701v+50.0594模型一 226.4 用曲線擬合方法確定醫(yī)用薄膜滲透率用曲線擬合方法確定醫(yī)用薄膜滲透率擬合數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)比較 236.4 用曲線擬合方
12、法確定醫(yī)用薄膜滲透率用曲線擬合方法確定醫(yī)用薄膜滲透率v=20 25 30 35 40 45 50 55 60 65 70 75 80;d=42 56 73.5 91.5 116 142.5 173 209.5 248 292.5 343 401 464;v=sqrt(v);a=polyfit(v,d,1)dl=polyval(a,v);plot(v,d,ro,v,dl,b.); box off輸出結(jié)果:a=91.4428 -430.1865模型二通過觀測(cè)散點(diǎn)圖,還可以假設(shè)剎車距離d與v剎車速度之間的函數(shù)關(guān)系為 ,下面為編寫求未知參數(shù)的a1、a0的matlab程序01avad 246.4 用曲線
13、擬合方法確定醫(yī)用薄膜滲透率用曲線擬合方法確定醫(yī)用薄膜滲透率擬合函數(shù)為:1865.4304428.91vd擬合數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)比較 25簡(jiǎn)介曲面擬合簡(jiǎn)介曲面擬合多元最小二乘法已知m個(gè)變量(x1,.,xm)和一個(gè)因變量y的一組觀測(cè)值(x1i,.,xmi,yi),i=1,2,.,n,要確定函數(shù)y=f(x1,.,xm),使得niimiiiyxxxfJ1221),.,(min通過機(jī)理分析或數(shù)據(jù)的直觀判斷,去確定函數(shù)f(x1,.,xm)的結(jié)構(gòu),假定函數(shù)中含有未知參數(shù)a1,.,ak,再通過最小二乘原理具體確定參數(shù)a1,.,ak 。 26簡(jiǎn)介曲面擬合簡(jiǎn)介曲面擬合經(jīng)濟(jì)增長(zhǎng)模型在經(jīng)濟(jì)學(xué)中,有一個(gè)著名的Cobb-Do
14、uglas生產(chǎn)函數(shù):1,0 ,),(LaKLKQQ,K,L分別表示產(chǎn)值、資金、勞動(dòng)力,式中的,a要由經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)確定。現(xiàn)有美國(guó)馬薩諸塞州1900-1926年上述三個(gè)經(jīng)濟(jì)指數(shù)的統(tǒng)計(jì)數(shù)據(jù),試用最小二乘法,求出,a 27簡(jiǎn)介曲面擬合簡(jiǎn)介曲面擬合最小二乘的準(zhǔn)則2712)(iiiiQLaKmin1)建立M函數(shù)文件jingjizz.mfunction Q=jingjizz(x,y)Q=x(1)*(y(1,:).x(2).*(y(2,:).x(3);其中 x(1)=a;x(2)=;x(3)= 28簡(jiǎn)介曲面擬合簡(jiǎn)介曲面擬合2)建立運(yùn)行文件qumiannihe.mQ=1.05 1.18 1.29 1.30 1.
15、3 1.42 1.50 1.52 1.46 1.6 1.69 1.81 1.93 1.95 2.01 2.0 2.09 1.96 2.20 2.12 2.16 2.08 2.24 2.56 2.34 2.45 2.58;y=1.04 1.06 1.16 1.22 1.27 1.37 1.44 1.53 1.57 2.05 2.51 2.63 2.74 2.82 3.24 3.24 3.61 4.1 4.36 4.77 4.75 4.54 4.54 4.58 4.58 4.58 4.54;1.05 1.08 1.18 1.22 1.17 1.3 1.39 1.47 1.31 1.43 1.58 1.59 1.66 1.68 1.65 1.62 1.86 1.93 1.96 1.95 1.9 1.58 1.67 1.82 1.6 1.61 1.64;x0=0.1, 0.1, 0.2;x=lsqcur
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 跨國(guó)公司出口業(yè)務(wù)委托管理三方合同
- 倉(cāng)儲(chǔ)物流場(chǎng)地租賃及一體化服務(wù)合同
- 貝類養(yǎng)殖產(chǎn)業(yè)技術(shù)評(píng)價(jià)合同
- 碑刻與政治史研究合同
- 疏港工程拆除與港口設(shè)施重建合同
- 成都寫字樓租賃合同示范文本
- 美術(shù)課件制作技能
- 美術(shù)外寫生課件
- 安全風(fēng)險(xiǎn)管控的內(nèi)容
- 檢驗(yàn)崗位的安全職責(zé)
- 閱讀療法創(chuàng)新發(fā)展
- 第二十章手術(shù)減肥及體形塑造美容手術(shù)美容外科學(xué)概論講解
- 2024年消防文員理論考試題庫(kù)(濃縮400題)
- 山東省煙臺(tái)市招遠(yuǎn)市(五四學(xué)制)2024-2025學(xué)年八年級(jí)上學(xué)期期末考試物理試卷(含答案)
- 2024年05月陜西秦農(nóng)農(nóng)村商業(yè)銀行股份有限公司2024年招考金融科技勞務(wù)派遣人員筆試歷年參考題庫(kù)附帶答案詳解
- 工程施工質(zhì)量檢驗(yàn)試驗(yàn)計(jì)劃
- 物業(yè)管理委托合同變更協(xié)議
- 產(chǎn)品痛點(diǎn)及解決方案
- 防災(zāi)減災(zāi)培訓(xùn)(安全行業(yè)講座培訓(xùn)課件)
- 電力安全知識(shí)相關(guān)內(nèi)容培訓(xùn)
- 2025年醫(yī)療質(zhì)量安全管理與持續(xù)改進(jìn)考核標(biāo)準(zhǔn)、考核辦法、改進(jìn)措施
評(píng)論
0/150
提交評(píng)論