【正文】 九、實驗改進想法：無。subplot(2,1,2),plot(t,u2)。figure(2)。u23=abs(U(3)*H(3))*cos(101*t+angle(U(3)*H(3)))。u21=abs(U(1)*H(1))*cos(99*t+angle(U(1)*H(1)))。subplot(2,1,1),plot(w,abs(H)),grid。H=polyval(b,j*w)./polyval(a,j*w)。a=[1,2,10001]。U=[,1,]。t=linspace(0,2*pi,1001)。已知帶通濾波器的系統(tǒng)函數(shù)為H(s)=五、解題分析：需要知道求頻率響應的方法，并繪制圖形。九、實驗改進想法：無。七、實驗結(jié)果：八、實驗心得體會：通過對做信號的離散序列以及經(jīng)FFT和IFFT的變換，了解了相關特性。subplot(2,1,2)。plot(t,x)。y=real(ifft(fft(x)))。六、實驗程序：t=0:1/255:1。由于計算機內(nèi)寸的限制，MATLAB無法表示一個無窮長的序列。之后畫出序列波形。用MATLAB表示一離散序列，可用兩個向量來表示。連續(xù)時間函數(shù)與離散時間函數(shù)在編程中的區(qū)別主要體現(xiàn)在如下兩個方面：第一，自變量的取值范圍不同，離散時間函數(shù)的自變量是整數(shù)，而連續(xù)時間函數(shù)的自變量為一定范圍內(nèi)的實數(shù)；第二，繪圖所用的函數(shù)不同，連續(xù)函數(shù)圖形的繪制不止一個。實驗四一、實驗者姓名：二、實驗時間：三、實驗地點：四、實驗題目：將信號x(t)=sin(240pt)做離散序列，比較原序列與經(jīng)過FFT和IFFT變換后的序列，并做出說明。通過本次實驗，對幅值譜和相位譜有了更深的了解，并與課程《信號與系統(tǒng)》里的一些相關知識連接到了一起，使得學到的只是更加深刻、有意義。title(‘phase’)。subplot(3,1,3)。plot(q,abs(y))。title(‘source signal’)。subplot(3,1,1)。x=2*sin(4*pi*t)+5*cos(8*pi*t)。t=d*n。N=64。clear。ifft(X,N)計算N點IDFT。若MN，則將原序列截短為N點序列，再計算其N點DFT；若Mifft(X)計算M點的IDFT。M是序列x的長度，即M=length(x)。五、解題分析：在MATLAB信號處理工具箱中，MATLAB提供了4個內(nèi)部函數(shù)用于計算DFT和IDFT，它們分別是：fft(x)，fft(x,N)，ifft(X)，ifft(X,N)。九、實驗改進想法：無。p rad/sample)八、實驗心得體會：通過實驗我們知道了使用Matlab來繪出出一個線性移不變系統(tǒng)的幅頻和相頻曲線。freqz(b,a,w)七、實驗結(jié)果：沖擊響應圖及階躍響應圖:Impulse (sec)Step (sec)3035404550 100Magnitude(dB) Frequency(180。dstep(b,a,50)w=[0:1:512]*pi/512。subplot(2,1,1)。a=[1 ]。close。六、實驗程序:clc。在抽樣點w上的值。若X（ejW）可表示為b0+b1ejW+...bMejWMB(ejW)X(e)==jWjWjWNA(e)a0+a1e+...+aNe則freqz的調(diào)用形式為 X=freqz(b,a,w)，其中的b和 a分別是表示前一個jW式子中分子多項式和分母多項式系數(shù)的向量，即a=[a0,a1,...,aN]，w為抽樣的頻率點，向量w的長度至少為2。當序列的DTFT可寫成ejW的有理多項式時，可用MATLAB信號處理工具箱提供的freqz函數(shù)計算DTFT的抽樣值。X表示輸入序列，y表示輸出序列。filter函數(shù)調(diào)用形式為 y=filter(b,a,x)。(an/a0)y[kn]+229。bx[kn]nnn=0n=0NN其中x[k]、y[k]分別系統(tǒng)的輸入和輸出。許多離散LTI都可用如下的線性常系數(shù)的差分方程描述229。MATLAB函數(shù)conv的返回值y中只有卷積的結(jié)果，沒有y的取值范圍。五、解題分析：離散卷積是數(shù)字信號處理中的一個基本運算，MTLAB提供的計算兩個離散序列卷積的函數(shù)是conv，其調(diào)用方式為 y=conv(x,h)。九、實驗改進想法：無。了解到了N階系統(tǒng)的單位階躍響應的計算方法，和系統(tǒng)的響應波形圖的函數(shù)應用和繪制方法。step(num,den,10)。六、試驗程序：num=[5 25 30]。第五篇：信號與系統(tǒng)仿真實驗報告信號與系統(tǒng)仿真實 驗 報 告班級： 學號： 姓名： 學院：實驗一一、實驗者姓名：二、實驗時間：三、實驗地點：四、實驗題目：5(s2+5s+6)求三階系統(tǒng)H(s)=3的單位階躍響應，并繪制響應波形圖。Matlab提供了強大的用戶圖形界面，以幫助用戶不必編寫底層程序而直接在軟件包基礎上進行自行開發(fā)，這點在諸多軟件中都有所體現(xiàn)。本次試驗收獲很多，學習到了很多知識，首先是熟悉了SIMULINK由簡至繁搭建系統(tǒng)的過程，學習了利用VR建立虛擬模型，并在SIMULINK中連接。6實驗總結(jié)與心得此次MATLAB實驗綜合了SIMULINK、GUI和VR場景等多個部分，對四旋翼飛行器運動進行了仿真模擬。使用默認數(shù)據(jù)，此時xdes=3，ydes=4，zdes=5，開始仿真，可以得到運動軌跡x、y、z的響應函數(shù)，同時可以得到在xyz坐標中的空間運動軌跡。通過俯仰角和滾轉(zhuǎn)角控制飛行器在XW和YW平面上的運動，通過控制偏航角，通過控制飛行器在ZB軸上的運動?？赏ㄟ^PID控制器計算控制信號：是目標懸停位置是我們的目標懸停位置(i=1,2,3)，是期望加速度，即控制信號。外回路：位置控制回路外回路采用以下控制方式：通過位置偏差計算控制信號（加速度）；建立控制信號與姿態(tài)角之間的幾何關系；得到期望姿態(tài)角，作為內(nèi)回路的輸入。因此，這里首先對轉(zhuǎn)速產(chǎn)生的作用進行分析。RigidBodyDynamics是被控對象，模擬四旋翼飛行器的運動特性。內(nèi)回路由AttitudeControl模塊實現(xiàn)，輸入為期望姿態(tài)角，輸出為期望轉(zhuǎn)速。外回路由PositionControl模塊實現(xiàn)。其中，L為旋翼中心建立飛行器質(zhì)心的距離，I為慣量矩陣。四旋翼飛行器所受力矩來源于兩個方面：1）旋翼升力作用于質(zhì)心產(chǎn)生的力矩；2）旋翼旋轉(zhuǎn)產(chǎn)生的扭轉(zhuǎn)力矩。線運動由合外力引起，符合牛頓第二定律：r為飛機的位置矢量。期望轉(zhuǎn)速則需要限制在電機的最小轉(zhuǎn)速和最大轉(zhuǎn)速之間，范圍可分取[1200rpm，7800rpm]。當給定期望轉(zhuǎn)速后，電機的實際轉(zhuǎn)速需要經(jīng)過一段時間才能達到。旋翼旋轉(zhuǎn)產(chǎn)生旋轉(zhuǎn)力矩Mi(i=1,2,3,4)，力矩Mi的旋向依據(jù)右手定則確定。力模型為：，旋翼通過螺旋槳產(chǎn)生升力。四旋翼飛行器受力分析,如圖所示圖2四旋翼飛行器受力分析示意圖旋翼機體所受外力和力矩為：重力mg,機體受到重力沿方向；四個旋翼旋轉(zhuǎn)所產(chǎn)生的升力(i=1,2,3,4)，旋翼升力沿方向；旋翼旋轉(zhuǎn)會產(chǎn)生扭轉(zhuǎn)力矩(i=1,2,3,4)。圖1四旋翼飛行器旋轉(zhuǎn)方向示意圖在圖中，前端旋翼和后端旋翼逆時針旋轉(zhuǎn)，而左端旋翼和右端的旋翼順時針旋轉(zhuǎn)，以平衡旋翼旋轉(zhuǎn)所產(chǎn)生的反扭轉(zhuǎn)矩。四個旋翼位于一個幾何對稱的十字支架前，后，左，右四端，如圖所示。工具軟件：操作系統(tǒng)：Windows系列；軟件工具：MATLAB及simulink。2實驗目的通過在Matlab環(huán)境中對四旋翼飛行器進行系統(tǒng)建模，使掌握以下內(nèi)容：四旋翼飛行器的建模和控制方法在Matlab下快速建立虛擬可視化環(huán)境的方法。我也基本掌握了Comsol Multiphysics這款軟件，我相信在今后我會將我對本課程的學習運用到實際中。通過對仿真實驗課程的學習，及本次試驗，我體會到仿真技術對于實驗的幫助非常巨大，使得實驗室進行的許多實驗可以通過計算機模擬直接完成，節(jié)省了資源消耗，并極大地提高了實驗效率。我在樣品上樣時一直未能獲得良好的上樣結(jié)果，發(fā)現(xiàn)對瞬態(tài)求解器的時間比例進行修改，可以獲得良好上樣結(jié)果，同時，在樣品分離改變chefk左進口濃度時發(fā)現(xiàn)修改數(shù)值導致結(jié)果錯誤，遂未修改濃度，得到了正確結(jié)果。求解模式全部為當前解，對電動流求解，結(jié)果用濃度表示。求解模式為初始值表達式和當前解，對傳到介質(zhì)DC求解，結(jié)果用電位能表示。（2）改變cmdc進口，7邊界“電位能/10”，1邊界“電位能/3”。再求兩次解，完成上樣。（2）改變chekf進口，7邊界改為“通量/向內(nèi)通量nmflux_c_chekf”；求解域中x速度和y速度改為0去除載流作用；求解器設置改為瞬態(tài)求解器，時間改為“0::”。樣品上樣（1）改變emdc進口，邊界7電位能由10改為3。（4）再求解電場，改變求解模式，點變量值不可解和線性化/當前解，對傳導介質(zhì)DC求解，觀察求解結(jié)果，用電位能表示。（2）求解器管理器設置求解模式：初始值/初始值表達式，點變量值不可解和線性化/從初始值使用設定。設置邊界條件mmglf—入口1和7邊界“進口/層流流進/”出口5和12邊界“出口/壓力，粘滯應力/0”；emdc—入口1和7邊界“電位能/10V”出口5和12邊界“接地”其余邊界“電絕緣”；chekf—入口1“濃度/1”，7“濃度/0”出口5和12“通量/向內(nèi)通量nmflux_c_chekf”其余邊界“絕緣/對稱”。網(wǎng)格劃分菜單欄，網(wǎng)格，自由網(wǎng)格參數(shù)，通常網(wǎng)格尺寸，最大單元尺寸：4e7。5和9兩端點取圓角，半徑1e6。復制該矩形，旋轉(zhuǎn)90176。Better Defaults第三篇：仿真實驗報告仿真軟件實驗實驗名稱：基于電滲流的微通道門進樣的數(shù)值模擬實驗日期：、實驗目的對建模及仿真技術初步了解學習并掌握Comsol Multiphysics的使用方法了解電滲進樣原理并進行數(shù)值模擬運用Comsol Multiphysics建立多場耦合模型，加深對多耦合場的認識二、實驗設備實驗室計算機，Comsol Multiphysics 。Quicker Defaultsl Better Defaults做法第一步：選擇Better Defaults參數(shù)第二步：分析結(jié)果。l Quicker Defaults做法第一步：選擇Quicker Defaults參數(shù),并開始運行第二步：分析結(jié)果。Better Defaults：將庫存檢查間隔ei克隆拖入Optimizer模塊上，在訂貨模型create模塊后加隊列防止其堵塞，ei的范圍設置如下。：選擇Quicker Defaults參數(shù)，（s，S）優(yōu)化結(jié)果為（27,43）。分別用Quicker Defaults參數(shù)和Better Defaults參數(shù)運行兩次比較結(jié)果。運行結(jié)束后選擇Better Default按鈕設置規(guī)模更大的優(yōu)化參數(shù)（耗時長精度高），并比較兩次的結(jié)果。 Init模塊初始化最小最大庫存(s，S) 從value庫中將Opt