【正文】 出新的定點(diǎn)，進(jìn)行細(xì)化。也就是引進(jìn)新的頂點(diǎn)。具體可以去看蘭州大學(xué)信息工程學(xué)院的學(xué)報(bào)。首先解釋一下什么叫做徑向插值 徑向基函數(shù)主要是考慮多維空間的數(shù)據(jù)插值問題，徑向基函數(shù)在三維圖形的變形中常常被采用，用其來變分三維人臉的時(shí)候還應(yīng)結(jié)合人臉的固有特征。插值方法：假設(shè)大多數(shù)人臉的形狀都可以由一個(gè)拓?fù)湓妥兓脕?，那么，通過調(diào)整一個(gè)一般模型的構(gòu)造參數(shù)可以建立不同的面部模型。在離散數(shù)據(jù)的多變量插值問題方面，徑向基函數(shù)（radial basis function，RBF）插值方法是一個(gè)行之有效的工具，所以也適用于類似人臉這樣高維曲面的近似或平滑插值。這種方法的優(yōu)點(diǎn)在于：（1）通過插值可以得到丟失的數(shù)據(jù)點(diǎn)，所以源網(wǎng)格和目標(biāo)網(wǎng)格不需要相同數(shù)目的結(jié)點(diǎn)；（2）如果選擇了合適的匹配點(diǎn)，數(shù)學(xué)上可以保證能夠?qū)⒃淳W(wǎng)格變形到目標(biāo)網(wǎng)格。這種方法的關(guān)鍵是找到一個(gè)很好的核函數(shù)來計(jì)算出新的頂點(diǎn)。然后而且在網(wǎng)上找到了部分代碼。這是我校大四一個(gè)同學(xué)畢設(shè)中的一個(gè)東西。第三次遇到的問題是我在數(shù)學(xué)建模課上遇到的一個(gè)關(guān)于矩陣LU分解的問題：Matlab作出的結(jié)果和手算的結(jié)果竟然不一樣。249。233。249。11一個(gè)矩陣 A=234。ee11U=[epusino 1。L=[1 0。L*U ans =0.***1.***1.***0 而直接用lu函數(shù)對(duì)矩陣A分解的結(jié)果為A=[epusino 1。翻譯過來是L是一個(gè)心理上的下三角矩陣，其實(shí)是下三角矩陣和置換矩陣也就是最后的結(jié)論。由此可見這并不是什么精度的問題導(dǎo)致的結(jié)果只是LU分解的矩陣意義和書本上的不同而已。首先的理解是利用鄰接矩陣，先隨機(jī)生成一個(gè)矩陣，該矩陣為一個(gè)對(duì)稱矩陣。隨機(jī)拓?fù)鋱D10026516398191 ***0430 ***065%產(chǎn)生數(shù)組A用來存放表示兩點(diǎn)之間權(quán)值的矩陣A，也就是臨接矩陣，那么兩點(diǎn)之間權(quán)值不為零元素的個(gè)數(shù)即為該點(diǎn)的度數(shù)DEF=5。%產(chǎn)生DEF*DEF的隨機(jī)矩陣 for i=1:DEF A(i,i)=0%將對(duì)角線上的數(shù)置為0 end A=10*A。%向下去整 for i=1:DEF for j=1:i A(j,i)=A(i,j)%將A矩陣變?yōu)橐粋€(gè)上三角或者下三角矩陣 end end x=100*rand(1,DEF)。%產(chǎn)生10個(gè)隨機(jī)的點(diǎn) plot(x,y,39。)。%將A中的權(quán)值轉(zhuǎn)化為字符型if c~=39。%不顯示為0的值 因?yàn)锳矩陣為零代表兩點(diǎn)不相連 text((x(i)+x(j))/2,(y(i)+y(j))/2,c,39。,18)。line([x(i)x(a(j))],[y(i)y(a(j))])。隨機(jī)拓?fù)鋱D39。e=num2str(DEF)。%左上角顯示節(jié)點(diǎn)的個(gè)數(shù) for m=1:DEF A(m,m)=m。hold on。Fontsize39。%將權(quán)值顯示在兩點(diǎn)連線中間 end 接下來是一個(gè)在信息安全課上寫的一個(gè)關(guān)于256色圖分層的程序，一副圖像有m*n個(gè)像素然后每個(gè)像素是一個(gè)8bit的二進(jìn)制數(shù)據(jù)換算為十進(jìn)制是0255之間。I=imread(39。)。%將圖像轉(zhuǎn)換為double類型便于MATLAB中的計(jì)算 %確定圖像的長(zhǎng)寬 M=size(I,1)。%寬 l=M*N。%剔除每層 for i=1:Mfor j=1:NB=numdec2bin(I(i,j),8)。%刪除相應(yīng)位置的元素y(i,j)=numbin2dec(B)。%將10進(jìn)制的矩陣轉(zhuǎn)換為無符號(hào)整形subplot(4,2,9n)。%顯示這個(gè)圖像title(strcat(39。,num2str(9n),39。))。%從函數(shù)將十進(jìn)制數(shù)x轉(zhuǎn)化為n位二進(jìn)制 y=zeros(1,n)。i=n。a=floor(a/2)。end function y=numbin2dec(x)%將二進(jìn)制數(shù)轉(zhuǎn)換為十進(jìn)制數(shù) a=0。for i=1:lx a=a+x(i)*2^(lxi)。去除第1層后的結(jié)果去除第2層后的結(jié)果去除第3層后的結(jié)果去除第4層后的結(jié)果去除第5層后的結(jié)果去除第6層后的結(jié)果去除第7層后的結(jié)果去除第8層后的結(jié)果Matlab的RGB轉(zhuǎn)換為YCbCr之后轉(zhuǎn)換回來和原圖不符這個(gè)是做JPEG圖像壓縮的時(shí)候遇到的一個(gè)問題，當(dāng)時(shí)的代碼是用C寫的但是結(jié)果不對(duì)，后來想到用matlab進(jìn)行驗(yàn)證結(jié)果發(fā)現(xiàn)通過公式轉(zhuǎn)換得到的圖是錯(cuò)的。Y= R + G + B Cb = + B + 128 Cr = B + 128 然后再進(jìn)行量化，DCT，編碼等等步驟。R = Y + ((Cb128)128)B = Y + (Cb128)但是通過書上給的公式和網(wǎng)上大部分公式卻發(fā)現(xiàn)根本轉(zhuǎn)化不會(huì)來。clc,clear。39。%讀入原始RGB圖像figure(1)。imshow(Source):title(39。)。%計(jì)算圖像大小 %計(jì)算紅色分量并顯示分解圖% R(:,:,1)=Source(:,:,1)。R(:,:,3)=zeros(r,c)。whos。subplot(1,3,1)。Red Component39。%計(jì)算綠色分量并顯示分解圖% G(:,:,2)=Source(:,:,2)。G(:,:,3)=zeros(r,c)。figure(2)。imshow(G)title(39。)。B(:,:,1)=zeros(r,c)。B=uint8(B)。subplot(1,3,3)imshow(B)title(39。)。Comp(:,:,2)=G(:,:,2)。figure(1)。imshow(Comp):title(39。)。Cb=***G+128。red=Y+*(Cr128)。blue=Y+*(Cb128)。Comp2(:,:,2)=green(:,:,2)。figure(3)。title(39。)。GD=G(:,:,2)green(:,:,2)。figure(4)。imshow(RD)。紅色分量差異39。subplot(1,3,2)。title(39。)。imshow(BD)。藍(lán)色分量差異39。然而可以通過figure(4)的圖像清楚看到R 和G分量轉(zhuǎn)換后恢復(fù) 與以前的差別很大。很是不解實(shí)驗(yàn)得到的圖為: original imageposition image將源圖像分為RGB三個(gè)部分。紅色分量差異綠色分量差異藍(lán)色分量差異 公式轉(zhuǎn)換的y分量公式轉(zhuǎn)換的cb分量公式轉(zhuǎn)換的Cr分量matlab函數(shù)轉(zhuǎn)換的y分量matlab函數(shù)轉(zhuǎn)換的Cb分量matlab函數(shù)轉(zhuǎn)換的Cr分量從最后的結(jié)果可以看出原圖轉(zhuǎn)換到Y(jié)CbCr空間之后又轉(zhuǎn)化回來得到的圖像與原圖不符….RGB轉(zhuǎn)換為YCrCb后又轉(zhuǎn)換為RGB的圖像以上只是我自己用Matlab做過的部分問題，還有一些比如場(chǎng)聲源定位Music算法仿真，DES加密算法。二 Simulink學(xué)習(xí)心得與編程實(shí)踐SIMILINK模塊庫(kù)按功能進(jìn)行分類，包括以下8類子庫(kù)： Continuous（連續(xù)模塊）Discrete（離散模塊）Functionamp。Systems（信號(hào)和系統(tǒng)模塊）Sinks（接收器模塊）Sources（輸入源模塊）對(duì)于simulink來說其實(shí)沒什么好說的，需要什么就把什么拖到Model里面，連線之后設(shè)定參數(shù)。但是重要的是對(duì)于系統(tǒng)的設(shè)計(jì)與仿真思想才是最重要的。Simulink的工作方式：（1）模塊內(nèi)的參數(shù)值首先會(huì)送到Matlab中進(jìn)行計(jì)算，得到的參數(shù)值會(huì)用來當(dāng)做以后需要調(diào)用的參數(shù)。（3）模塊按被處理的順序排列，此時(shí)代數(shù)回路結(jié)構(gòu)也被檢查出來，此種排列產(chǎn)生一個(gè)列表，以確保具有代數(shù)回路的模塊驅(qū)動(dòng)輸入的模塊被更新后才更新?，F(xiàn)在可以準(zhǔn)備執(zhí)行仿真操作，仿真時(shí)使用數(shù)值迭代求的的結(jié)果，每種數(shù)值積分模型提供的連續(xù)狀態(tài)的微分能力。定義了一個(gè)模型以后，就可以通過Simulink的菜單或者在Matlab的Command中輸入命令進(jìn)行仿真。1：觀察一個(gè)信號(hào)與積分之后的區(qū)別，目的是了解Scope的功能和用法。用XYGraph看到的對(duì)比Rossler吸引子產(chǎn)生仿真：Rossler 吸引子產(chǎn)生是服從下面這個(gè)方程的 1212 331仿真圖 amp。=x+axxamp。仿真結(jié)果圖形：1：Spectrum Scope的輸出圖形2：Vector Scope的輸出結(jié)果3：Scope的輸出結(jié)果圖形第三篇：MATLAB與控制系統(tǒng)仿真實(shí)驗(yàn)報(bào)告《MATLAB與控制系統(tǒng)仿真》實(shí)驗(yàn)報(bào)告20132014學(xué)年 第 1 學(xué)期專業(yè)： 班級(jí)： 學(xué)號(hào)： 姓名：實(shí)驗(yàn)三 MATLAB圖形系統(tǒng)一、實(shí)驗(yàn)?zāi)康模?。?）繪制多根二維曲線 plot(x,y)當(dāng)x是向量，y是有一維與x同維的矩陣時(shí)，則繪制多根不同顏色的曲線。（3）含有多個(gè)輸入?yún)?shù)的plot函數(shù) plot(x1,y1,x2,y2,…,xn,yn)（4）具有兩個(gè)縱坐標(biāo)標(biāo)度的圖形 plotyy(x1,y1,x2,y2)2,圖形標(biāo)注與坐標(biāo)控制 1）title(圖形名稱)； 2）xlabel（x軸說明）3）ylabel（y軸說明）4）text（x，y圖形說明）5）legend（圖例1，圖例2，…）6）axis（[xmin xmax ymin ymax zmin zmax]）3, 圖形窗口的分割 subplot（m,n,p）4,三維曲線plot3（x1,y1,z1,選項(xiàng)1，x2,y2，選項(xiàng)2，…,xn,yn,zn,選項(xiàng)n）5，三維曲面mesh(x,y,z,c)與surf(x,y,z,c)。X，y是網(wǎng)格坐標(biāo)矩陣，z是網(wǎng)格點(diǎn)上的高度矩陣，c用于指定在不同高度下的顏色范圍。2）image和imagesc函數(shù) 這兩個(gè)函數(shù)用于圖象顯示。三、實(shí)驗(yàn)儀器和設(shè)備：計(jì)算機(jī)一臺(tái)（）。2．復(fù)習(xí)圖形輔助操作。21+x2，已知y1=x2,y2=cos(2x),y3=y1*y2,完成下列操作：（1）在同一坐標(biāo)系下用不同的顏色和線型繪制三條曲線；（2）分別用條形圖、階梯圖、桿圖和填充圖繪制三條曲線。x+p,239。02239。239。238。z=x2+y2+sin(xy)的曲面及等高線圖。+a25a2四、心得體會(huì)：通過這次實(shí)驗(yàn)我能熟練掌握二維和三維圖以及其他特殊圖形的制作，弄清楚了基本的圖形操作規(guī)則，大大加深了我對(duì)matlab的興趣。三、實(shí)驗(yàn)內(nèi)容+bx+c=0的根的函數(shù)（這個(gè)方程不一定為一元二次方程，因a、b、c的不同取值而定），這里應(yīng)根據(jù)a、b、c的不同取值分別處理，有輸入?yún)?shù)提示，當(dāng)a=0,b=0,c~=0時(shí)應(yīng)提示“為恒不等式!”。clear,clc a=input(39。)。請(qǐng)輸入一個(gè)數(shù)b=39。c=input(39。)。if a==0if b==039。end endm=b^24*a*c。不存在正實(shí)根39。其中100分為A+，90分～99分為A，80分～89分為B，70分～79分為C，60分～69分為D，60分以下為E。（2）輸入百分制成績(jī)后要判斷該成績(jī)的合理性，對(duì)不合理的成績(jī)應(yīng)輸出出錯(cuò)信息。b(k)={79+k}。d(k)={59+k}。A(1,:)={39。,39。,39。,39。,39。,39。}。for k=1:6switch A{2,k}case 100r=39。case ar=39。case br=39。case cr=39。case dr=39。otherwiser=39。endA(3,k)={r}。a39。b39。c39。B39。C39。A39。d39。e39。A+39。E39。f39。D39。clear,clc n=input(39。)。輸入錯(cuò)誤39。n!=039。n!=139。for i=1:1:ny=y*i。endy end 請(qǐng)輸入一個(gè)正數(shù)n=89ans =輸入錯(cuò)誤 請(qǐng)輸入一個(gè)正數(shù)n=0ans =n!=0 請(qǐng)輸入一個(gè)正數(shù)n=1ans =n!=1 請(qǐng)輸入一個(gè)正數(shù)n=3y =6 請(qǐng)輸入一個(gè)正數(shù)n=10y =3628800四、實(shí)驗(yàn)心得體會(huì)：通過本次實(shí)驗(yàn)課，我能熟練運(yùn)用for循環(huán)語句，switch條件語句以及if條件語句的新用法，和在C中的區(qū)別。第四篇：金陵石化電力系統(tǒng)故障分析報(bào)告金陵石化2013年1月5日電力系統(tǒng)故障匯報(bào)一、電力系統(tǒng)故障過程及影響故障過程2013年1月5日故障發(fā)生前，金陵分公司供電系統(tǒng)處正常運(yùn)行，即兩路地區(qū)電網(wǎng)聯(lián)絡(luò)線和金陵分公司三臺(tái)發(fā)電機(jī)均正常運(yùn)行。致使南京東北區(qū)域電網(wǎng)220kV系統(tǒng)電壓波動(dòng)，錄波監(jiān)控顯示，自05：06：18起到05：05：22，金陵石化熱電變電站110kV兩母線電壓多次波動(dòng)，最低電壓低于40%額定電壓，最后一次波動(dòng)持續(xù)時(shí)間約1秒。堯化門變電站2主變跳停后，金陵石化熱電Ⅲ發(fā)電機(jī)組（100MW）與南京市電網(wǎng)解列，孤網(wǎng)運(yùn)行（因堯金Ⅰ線721未停，小網(wǎng)包括堯化門變部分負(fù)荷），因小網(wǎng)負(fù)荷嚴(yán)重超過機(jī)組負(fù)載能力，Ⅲ發(fā)電機(jī)組強(qiáng)勵(lì)動(dòng)作引起勵(lì)磁變過流，在05：06：28秒停機(jī)。電氣專業(yè)檢查發(fā)現(xiàn)6kVⅤ、Ⅵ段電源63636斷路器跳閘，備用電源60606斷路器未能自投，6kVⅤ、Ⅵ段及400VⅤ、Ⅵ段失電。鍋爐專業(yè)檢查發(fā)現(xiàn)Ⅴ、Ⅵ爐跳停，Ⅰ、Ⅱ、Ⅲ、Ⅳ爐從觀火孔觀察到爐膛內(nèi)仍有火光，通過Ⅰ～Ⅳ爐盤前電接點(diǎn)水位計(jì)觀察到水位上升很快，令外操就地手操各爐03門，但水位依然上升至+300，于是按照運(yùn)行規(guī)程，按下停爐按鈕將Ⅰ～Ⅳ爐停爐。此時(shí)熱電運(yùn)行部六爐三機(jī)均停，金陵分公司供電系統(tǒng)靠南京市電網(wǎng)堯金Ⅱ線724供電。(3機(jī)跳停后1秒，)2)備自投未自動(dòng)投入的裝置變電所有8個(gè)，其中煉油區(qū)域6個(gè)，煤化工區(qū)域2個(gè)。1)烷一車間112變電所高壓Ⅰ段、Ⅱ段進(jìn)線開關(guān)及11變低壓400V進(jìn)線開關(guān)未受影響，但部分高壓電機(jī)欠壓保護(hù)動(dòng)作。2)烷二車間的13變高壓Ⅰ段、Ⅱ段進(jìn)線開關(guān)未受影響，但部分電機(jī)受晃電影響而停機(jī)。10個(gè)低壓進(jìn)線開關(guān)全部跳閘，母聯(lián)開關(guān)均未自投，低壓設(shè)備全部停電。此次分公司電力系統(tǒng)事故，各部門在事故處理上符合預(yù)案要求，無次生事故發(fā)生，特別是分公司總負(fù)荷在事故后單電源運(yùn)行時(shí)控制在80MW左右，避免了系統(tǒng)故障影響進(jìn)一步擴(kuò)大現(xiàn)象。此時(shí)，金陵公司范圍內(nèi)大部分低壓電機(jī)均已跳停，大機(jī)組同時(shí)停運(yùn)，裝置停工。高壓電機(jī)按使用特點(diǎn)和機(jī)組特性作適當(dāng)延時(shí)跳停,~2秒。熱電廠用電6kVⅤ、Ⅵ失電及Ⅴ、Ⅵ爐