1、上海電機學院 盧昌鈺 BG0801 10號1.單相半波可控整流電路（1）電阻性負載（R=1歐姆，U2=220V,=30）接線圖電阻性負載二次電壓，輸出電壓，二次電流，輸出電流，晶閘管電壓曲線輸入電壓與輸出電壓波形（2）阻感負載（R=1歐姆，L=0.05H，U2=220V,=30）接線圖阻感負載二次電壓，輸出電壓，二次電流，輸出電流，晶閘管電壓曲線輸入電壓與輸出電壓波形（3）阻感負載+續流二極管（R=1歐姆，L=0.05H，U2=220V,=30）有問題接線圖阻感負載二次電壓，輸出電壓，二次電流，輸出電流，晶閘管電壓曲線輸入與輸出電壓波形2.單相橋式全控整流電路（1）電阻性負載（R=1歐姆，U2

2、=220V,=60）電阻性負載電路圖搭建電阻負載輸入電壓和輸出電壓對比 電阻負載直流電壓和電流波形電阻負載時晶閘管T1的波形電流i2的曲線（2）電感性負載（R=1歐姆，L=0.05H,=60，U2=220V,）阻感負載電路圖搭建阻感負載電壓輸入與輸出波形 阻感負載輸出電流id阻感負載輸出電壓ud阻感負載交變時的電流i2阻感負載交變時的電壓u2阻感負載VT1的電壓波形（3）電感性負載+續流二極管（R=1歐姆，L=0.05H,=60，U2=220V,）電感性負載+續流二極管接線圖輸入和輸出電壓波形負載電流負載電壓二次側電流晶閘管兩端電壓3.單相橋式半空整流電路（1）電阻負載（R=1歐姆，=60，U

3、2=220V,）接線圖二次側電壓，負載電壓，二次側電流，負載電流，晶閘管電壓，二極管電壓，二極管電流波形圖（2）阻感負載（R=1歐姆，L=0.05H，=60，U2=220V,）接線圖二次側電壓，負載電壓，二次側電流，負載電流，晶閘管電壓，二極管電壓，二極管電流波形圖（3）阻感負載+續流二極管（R=1歐姆，L=0.05H，=60，U2=220V,）接線圖二次側電壓，負載電壓，二次側電流，負載電流，晶閘管VT1電壓，二極管VD4電壓，二極管VD4電流波形圖4.三相半波可控整流電路電阻負載接線圖（0）三相輸入電壓 輸出電流和電壓 晶閘管1的電流電壓輸出波形（電阻0）三相輸入電壓 輸出電流和電壓 晶閘

4、管1的電流電壓輸出波形（電阻30）阻感負載接線圖（30）三相輸入電壓 輸出電流和電壓 晶閘管1的電流電壓輸出波形（阻感30）阻感負載+續流二極管接線圖（30）三相輸入電壓 輸出電流和電壓 晶閘管1的電流電壓輸出波形（阻感+續流二極管30）5三相全控整流電路電阻負載接線圖（30導通角）三相輸入輸出電壓對比，晶閘管1電壓，輸出電流電壓圖形（30）阻感負載接線圖（30導通角）三相輸入輸出電壓對比，晶閘管1電壓，輸出電流電壓圖形（30）阻感負載+續流二極管接線圖（30導通角）三相輸入輸出電壓對比，晶閘管1電壓，輸出電流電壓圖形（30）6 降壓BUCK電路降壓斬波電路（電流連續）接線圖BUCK變換器電感

5、電流連續時仿真波形BUCK變換器電感電流斷續時仿真波形7 升壓Boost電路升壓Boost變換器仿真接線圖升壓Boost變換器連續工作升壓Boost變換器斷續工作8 單相全橋方波逆變電路單相全橋方波電阻負載逆變電路接線圖電阻負載逆變器直流側電流，輸出交流電壓電流方波波形單相全橋方波阻感負載逆變電路接線圖阻感負載逆變器直流側電流，輸出交流電壓電流方波波形9 三相方波逆變電路三相方波逆變電路接線圖三相方波逆變電路仿真波形（感性無功=100Var）三相方波逆變電路仿真波形（感性無功=1000Var）10單極性的PWM方式下的單相全橋逆變電路單極性的PWM方式下的單相全橋逆變電路接線圖在下：輸出電壓，

6、電流和直流側電流波形附錄資料：MATLAB函數和命令的用法Binocdf二項式累積分布函數語法格式Y = binocdf(X,N,P)函數功能Y = binocdf(X,N,P) 計算X中每個X(i)的二項式累積分布函數，其中，N中對應的N(i)為試驗數，P中對應的P(i)為每次試驗成功的概率。Y, N, 和 P 的大小類型相同，可以是向量、矩陣或多維數組。輸入的標量將擴展成一個數組，使其大小類型與其它輸入相一致。The values in N must all be positive integers, the values in X must lie on the interval 0,N

7、, and the values in P must lie on the interval 0, 1.The binomial cdf for a given value x and a given pair of parameters n and p is The result, y, is the probability of observing up to x successes in n independent trials, where the probability of success in any given trial is p. The indicator functio

8、n I(0,1,.,n)(i)ensures that x only adopts values of 0,1,.,n.示例若一個棒球隊在一個賽季要比賽162場，每場比賽取勝的機會是50-50，則該隊取勝超過100 場的概率為： 1-binocdf(100,162,0.5)ans = 0.0010433相關函數 HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binofit.html binofit | HYPERLINK jar:file:/C:/Program%20Fil

9、es/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binoinv.html binoinv | HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binopdf.html binopdf | HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binornd.html binornd | HYPERLINK jar:f

10、ile:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binostat.html binostat | HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/cdf.html cdf附：二項式分布(binomial distribution )定義二項分布的概率密度函數為where k is the number of successes in n trials of a Bernoulli pro

11、cess with probability of success p.The binomial distribution is discrete, defined for integers k = 0, 1, 2, . n, where it is nonzero.背景The binomial distribution models the total number of successes in repeated trials from an infinite population under the following conditions:Only two outcomes are po

12、ssible on each of ntrials.The probability of success for each trial is constant.All trials are independent of each other.The binomial distribution is a generalization of the HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/brn2ivz-2.html Bernoulli distribution; it

13、generalizes to the HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/brn2ivz-84.html multinomial distribution.參數Suppose you are collecting data from a widget manufacturing process, and you record the number of widgets within specification in each batch of100. You mi

14、ght be interested in the probability that an individual widget is within specification. Parameter estimation is the process of determining the parameter, p, of the binomial distribution that fits this data best in some sense.One popular criterion of goodness is to maximize the likelihood function. T

15、he likelihood has the same form as the binomial pdf above. But for the pdf, the parameters (nandp) are known constants and the variable isx. The likelihood function reverses the roles of the variables. Here, the sample values (the xs) are already observed. So they are the fixed constants. The variab

16、les are the unknown parameters. MLE involves calculating the value of p that give the highest likelihood given the particular set of data.The function HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binofit.html binofit returns the MLEs and confidence intervals fo

17、r the parameters of the binomial distribution. Here is an example using random numbers from the binomial distribution with n=100 and p=0.9. r = binornd(100,0.9)r = 85 phat, pci = binofit(r,100)phat = 0.85pci = 0.76469 0.91355The MLE for parameterp is0.8800, compared to the true value of0.9. The 95%

18、confidence interval forp goes from 0.7998 to0.9364, which includes the true value. In this made-up example you know the true value ofp. In experimentation you do not.示例The following commands generate a plot of the binomial pdf for n = 10 and p = 1/2.x = 0:10;y = binopdf(x,10,0.5);plot(x,y,+) 相關內容 HY

19、PERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/bqt29ct.html l bqt290a-1 Discrete Distributions附：二項式分布（網上）定義若某事件概率為p，現重復試驗n次，該事件發生k次的概率為:P=C(k,n)pk(1-p)(n-k)C(k,n)表示組合數，即從n個事物中拿出k個的方法數。二項分布的概念考慮只有兩種可能結果的隨機試驗，當成功的概率（）是恒定的，且各次試驗相互獨立，這種試驗在統計學上稱為貝努里試驗（Bernoullitrial）。如果

20、進行n次貝努里試驗，取得成功次數為X（X=0,1,n）的概率可用下面的二項分布概率公式來描述：P=C(X,n) X(1-)(n-X)式中的n為獨立的貝努里試驗次數，為成功的概率，（1-）為失敗的 HYPERLINK /doc/5399415.html t _blank 概率，X為在n次貝努里試驗中出現成功的次數，C(X,n)表示在n次試驗中出現X的各種組合情況，在此稱為二項系數（binomialcoefficient）。內容簡介二項分布，伯努里分布：進行一系列試驗,如果1. 在每次試驗中只有兩種可能的結果,而且是互相對立的；2. 每次實驗是獨立的，與其它各次試驗結果無關；3. 結果事件發生的

21、HYPERLINK /doc/5399415.html t _blank 概率在整個系列試驗中保持不變，則這一系列試驗稱為伯努力試驗。在這試驗中，事件發生的次數為一 HYPERLINK /doc/2273201.html t _blank 隨機事件,它服從二次分布。二項分布可以用于可靠性試驗。可靠性試驗常常是投入n個相同的式樣進行試驗T小時，而只允許k個式樣失敗，應用二項分布可以得到通過試驗的概率。一個事件必然出現，就說它100%要出現。100%=1，所以100%出現的含義就是出現的概率P=1。即必然事件的出現概率為1。若擲一枚硬幣，正面向上的結果的概率為0.5。反面向上的結果的概率也是0.5

22、。則出現正面向上事件或者反面向上事件的概率就是0.5+0.5=1，即二者必居其一。若擲兩次硬幣，由獨立事件的概率乘法定理那么兩次都是正面（反面）向上的概率是0.50.5=0.25。另外第一個是正第二個是反的出現概率也是0.50.5=0.25。同理第一個反第二個正的出現概率也是0.50.5=0.25。于是一正一反的概率是前面兩個情況的和，即0.25+0.25=20.25=0.5。它們的合計值仍然是1。binopdf二項分布的概率密度函數語法格式Y = binopdf(X,N,P)函數功能Y = binopdf(X,N,P) 計算X中每個X(i)的概率密度函數，其中，N中對應的N(i)為試驗數，P

23、中對應的P(i)為每次試驗成功的概率。Y, N, 和 P 的大小類型相同，可以是向量、矩陣或多維數組。輸入的標量將擴展成一個數組，使其大小類型與其它輸入相一致。N中的值必須是正整數，0 P 1 。對于給出的x和兩個參數n和p，二項分布概率密度函數為其中 q = 1 p。 y為n次獨立試驗中成功x次的概率，其中，每次試驗成功的概率為p。指示器函數 I(0,1,.,n)(x) 確保x 取值為 0, 1, ., n。 示例一個質量檢查技術員一天能測試200塊電路板。假設有 2% 的電路板有缺陷，該技術員在一天的測試中沒有發現有缺陷的電路板的概率是多少？binopdf(0,200,0.02)ans =

24、0.0176質量檢查技術員一天中最大可能檢測出有缺陷的電路板是多少塊？defects=0:200;y = binopdf(defects,200,.02);x,i=max(y);defects(i) ans =4相關函數 HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/binocdf.html binocdf | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/binofit.html

25、 binofit | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/binoinv.html binoinv | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/binornd.html binornd | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/binostat.htm

26、l binostat | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/pdf.html pdfdlmread將以ASCII碼分隔的數值數據文件讀入到矩陣中語法格式M = dlmread(filename)M = dlmread(filename, delimiter)M = dlmread(filename, delimiter, R, C)M = dlmread(filename, delimiter, range)函數功能M = dlmread(filename)讀取有分隔符的AS

27、CII數值數據文件filename，并把數據給到矩陣M中。文件名(filename)以字符串形式用單引號括起來。 dlmread 從文件的格式中推斷分隔符。M = dlmread(filename, delimiter)指定分隔符delimiter。使用 t 代表制表符 tab 作為分隔。M = dlmread(filename, delimiter, R, C)R和C指定了數據讀取數據，其左上角的值位于文件中的第R行，第C列。R和C從0開始，所以R=0, C=0 指向文件中的第一個值。M = dlmread(filename, delimiter, range)讀取由range = R1 C

28、1 R2 C2指定的區域塊， (R1,C1)是左上角，(R2,C2)是右下角。也可以使用電子表格表示法來指定，如range = A1.B7。備注All data in the input file must be numeric. dlmread does not read files that contain nonnumeric data, even if the specified rows and columns contain only numeric data.若沒有指定分隔符，當從文件格式中推斷分隔符時，連續的空格符當作一個分隔符對待。若指定了分隔符，則重復的分隔符將分別作為單獨

29、的分隔符對待。If you want to specify an R, C, or range input, but not a delimiter, set the delimiter argument to the empty string, (two consecutive single quotes with no spaces in between, ). For example,M = dlmread(myfile.dat, , 5, 2)In this case, dlmread treats repeated white spaces as a single delimiter

30、.Dlmread將用0填充沒有邊界的區域。有多行的數據文件，若以非空格分隔符結束，例如分號，則導入后會多產生全0的一列在最后。Dlmread在導入任何復數時，將其作為一個整體導入到一個復數單元中。下面是有效的復數格式：i|jExample: 5.7-3.1ii|jExample: -7j嵌入了空格的復數是不正確的格式，空格將被認為是分隔符。示例例1Export a 5-by-8 test matrix M to a file, and read it with dlmread, first with no arguments other than the filename:M = galler

31、y(integerdata, 100, 5 8, 0); dlmwrite(myfile.txt, M, delimiter, t)dlmread(myfile.txt)ans = 96 77 62 41 6 21 2 42 24 46 80 94 36 20 75 85 61 2 93 92 82 61 45 53 49 83 74 42 1 28 94 21 90 45 18 90 14 20 47 68Now read a portion of the matrix by specifying the row and column of the upper left corner:dlm

32、read(myfile.txt, t, 2, 3)ans = 92 82 61 45 53 42 1 28 94 21 90 14 20 47 68This time, read a different part of the matrix using a range specifier:dlmread(myfile.txt, t, C1.G4)ans = 62 41 6 21 2 80 94 36 20 75 93 92 82 61 45 74 42 1 28 94例2Export matrix M to a file, and then append an additional matri

33、x to the file that is offset one row below the first:M = magic(3);dlmwrite(myfile.txt, M*5 M/5, )dlmwrite(myfile.txt, M/3, -append, . roffset, 1, delimiter, )type myfile.txt40 5 30 1.6 0.2 1.215 25 35 0.6 1 1.420 45 10 0.8 1.8 0.4 2.6667 0.33333 21 1.6667 2.33331.3333 3 0.66667When dlmread imports t

34、hese two matrices from the file, it pads the smaller matrix with zeros:dlmread(myfile.txt) 40.0000 5.0000 30.0000 1.6000 0.2000 1.2000 15.0000 25.0000 35.0000 0.6000 1.0000 1.4000 20.0000 45.0000 10.0000 0.8000 1.8000 0.4000 2.6667 0.3333 2.0000 0 0 0 1.0000 1.6667 2.3333 0 0 0 1.3333 3.0000 0.6667

35、0 0 0替代作為dlmread的替代，可使用導入向導。選擇菜單 “File | Import Data” 激活導入向導。相關函數 HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/techdoc/help.jar%21/ref/dlmwrite.html dlmwrite | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/techdoc/help.jar%21/ref/textscan.html textscanfprintf寫數據到文本文件或輸出到命令窗口語法格式fprintf(fi

36、leID, format, A, .)fprintf(format, A, .)count = fprintf(.)函數功能fprintf(fileID, format, A, .) fileID為文件句柄，指定要寫入的文件；format是用來控制所寫數據格式的格式符；A是用來存放數據的矩陣。applies the format to all elements of array A and any additional array arguments in column order, and writes the data to a text file. fprintf uses the en

37、coding scheme specified in the call to HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/techdoc/help.jar%21/ref/fopen.html fopen.fprintf(format, A, .) 將A中的數據按格式format輸出到命令窗口。count = fprintf(.) 返回fprintf寫的字節數。輸入參數fileID取值為下面中的一種：從文件打開時所得到的文件句柄，為一個整數值。1 為標準輸出設備（屏幕）。2 for standard error.fileID缺省時，輸出到

38、屏幕，即1為默認值。format單引號括起來的是 “轉換控制字符串 ”，指定輸出的格式。由下面給出的各部分的組合而成：百分號后緊跟的格式字符，如%s 。輸出區域的寬度，精度，和其它選項（附加格式說明符）。要原樣輸出的文本（普通字符）。轉義字符，包括：單引號% 百分號反斜線a響鈴b退格符fForm feedn換行r回車t水平制表符v垂直制表符xN十六進制數, NN八進制數, N轉換字符和可選的操作符按以下順序排列（包括空格）：下表列出可用的轉換字符和限定類型字符。數據類型轉換說 明Integer, signed%d 或 %i十進制整數%ld 或 %li64-bit base 10 values%

39、hd 或 %hi16-bit base 10 valuesInteger, unsigned%u十進制整數%o八進制整數(octal)%x十六進制整數 (hexadecimal), 使用字母 af%X十六進制整數 (hexadecimal), 使用字母AF%lu%lo%lx or %lX64-bit values, 十、八、十六進制整數%hu%ho%hx or %hX16-bit values, 十、八、十六進制整數Floating-point number%f定點表示法%e指數表示法，如 3.141593e+00%E同 %e，但用E，如 3.141593E+00%g%e 或 %f的緊湊格式，

40、沒有尾部的0%G%E 或 %f的緊湊格式，沒有尾部的0%bx or %bX%bo%bu雙精度的十六、八、十進制數值如：%bx prints pi as 400921fb54442d18%tx or %tX%to%tu單精度的十六、八、十進制數值如：%tx prints pi as 40490fdbCharacters%c單個字符%s字符串附加格式說明符，包括：字段寬度Minimum number of characters to print. Can be a number, or an asterisk (*) to refer to an argument in the input lis

41、t. For example, the input list (%12d, intmax) is equivalent to (%*d, 12, intmax).精度For %f, %e, or %E: Number of digits to the right of the decimal point.Example: %6.4f prints pi as 3.1416For %g or %G Number of significant digits.Example: %6.4g prints pi as 3.142Can be a number, or an asterisk (*) to

42、 refer to an argument in the input list. For example, the input list (%6.4f, pi) is equivalent to (%*.*f, 6, 4, pi).FlagsActionFlagExampleLeft-justify.%-5.2fPrint sign character (+ or ).+%+5.2fInsert a space before the value. % 5.2fPad with zeros.0%05.2fModify selected numeric conversions:For %o, %x

43、, or %X, print 0, 0 x, or 0X prefix.For %f, %e, or %E, print decimal point even when precision is 0.For %g or %G, do not remove trailing zeros or decimal point.#%#5.0fIdentifierOrder for processing inputs. Use the syntax n$, where n represents the position of the value in the input list.For example,

44、 %3$s %2$s %1$s %2$s prints inputs A, B, C as follows: C B A B.The following limitations apply to conversions:Numeric conversions print only the real component of complex numbers.If you apply an integer or string conversion to a numeric value that contains a fraction, MATLAB overrides the specified

45、conversion, and uses %e.If you apply a string conversion (%s) to integer values, MATLAB converts values that correspond to valid character codes to characters. For example, %s converts 65 66 67 to ABC.Different platforms display exponential notation (such as %e) with a different number of digits in

46、the exponent.環境示例Windows1.23e+004UNIX1.23e+04不同的操作系統環境顯示的負零(-0) 不相同。轉換字符環境%e 或 %E%f%g 或 %GWindows0.000000e+0000.0000000Others-0.000000e+00-0.000000-0A數值或字符數組。示例例1 輸出多個數值和文本到屏幕。 B = 8.8 7.7 ; 8800 7700; fprintf(X is %4.2f meters or %8.3f mmn, 9.9, 9900, B)X is 9.90 meters or 9900.000 mmX is 8.80 mete

47、rs or 8800.000 mmX is 7.70 meters or 7700.000 mm注：對B是按列處理。例2 對雙精度數四舍五入取整，并輸出到屏幕。 a = 1.02 3.04 5.06; fprintf(%dn, round(a);135例3 寫一組數據到文本文件中。x = 0:.1:1;y = x; exp(x); % 打開文件exp.txt為可寫fid = fopen(exp.txt, w);fprintf(fid, %6.2f %12.8fn, y);fclose(fid); % 查看文件內容type exp.txtMATLAB導入功能,所有UNIX應用程序,和微軟Word

48、和寫字板都把 n 作為一個換行符指示器。然而,如果你要用微軟記事本讀取文件,請使用“ r n”移動到新行再寫。將之前的調用改寫如下：fprintf(fid, %6.2f %12.8frn, y);輸出結果不變。例4 在Windows系統上，將PC形式的指數表示法（三位數的指數）轉換到unix形式的表示法（兩位數字）,并打印數據到一個文件:a = 0.06 0.1 5 300 % 用sprintf轉換數值數據到文本，使用%ea_str = sprintf(%et,a) % use strrep to replace exponent prefix with shorter versiona_st

49、r = strrep(a_str,e+0,e+);a_str = strrep(a_str,e-0,e-); %調用fprintf打印更新的文本字符串fid = fopen(newfile.txt,w);fprintf(fid, %s, a_str);fclose(fid); % 查看文件內容type newfile.txt例5 在屏幕上顯示一個超鏈接。例6 創建一個字符矩陣并存入磁盤，再讀出賦值給另一個矩陣。 a=string; fid=fopen(d:char1.txt,w); fprintf(fid,%s,a); fclose(fid); fid1=fopen(d:char1.txt,r

50、t); b=fscanf(fid1,%s)b =string例7 用fprintf輸出一個數值矩陣。i = 1:10;square_root = sqrt(i);square = i.2;cube = i.3;% 拼接成一個大矩陣，行向量轉置成列向量out = i square_root square cube ;% 輸出數據for i = 1:10fprintf ( %2d %11.4f %6d %8dn,out(i,:);end%第1列整數輸出，占2格%第2列實數輸出，占11格，4位小數i = 1:10;square_root = sqrt(i);square = i.2;cube = i.3;% 拼接成一個大矩陣out = i; square_root; square; cube;% 輸出數據fprintf ( %2d %11.4f %6d %8dn,out);%將第1列輸出為第1行，格式對應；%將第2列輸出為第2行，格式對應；參考文獻1 Kernighan, B. W., and D. M. Ritchie, The C Programming Language, Second Edition, Prent