matlab課程設(shè)計——調(diào)幅廣播系統(tǒng)的仿真設(shè)計

1、<p><b>　　課程設(shè)計任務(wù)書</b></p><p>　　課程名稱_____________</p><p>　　院(系) </p><p>　　專業(yè)班級 </p><p>　　姓名 </p>&l

2、t;p>　　學(xué)號 </p><p>　　指導(dǎo)教師 </p><p><b>　　目 錄</b></p><p>　　1、引言……………………………………………………………………………3</p><p>　　1.1課程設(shè)計應(yīng)達(dá)到的目的……

3、………………………………………………………… 3</p><p>　　1.2 課程設(shè)計題目及要求……………………………………………………………………3</p><p>　　2、調(diào)頻廣播系統(tǒng)的模型及仿真環(huán)境………………………………………………… 4</p><p>　　2.1 MATLAB及SIMULINK建模環(huán)境簡介………………………………………………4<

4、/p><p>　　2.2 調(diào)幅廣播系統(tǒng)介紹………………………………………………………………………4</p><p>　　2.3 模型參數(shù)指標(biāo)……………………………………………………………………………4</p><p>　　2.3 仿真參數(shù)設(shè)計…………………………………………………………………………… 5</p><p>　　3、系統(tǒng)的建立與仿真…

5、………………………………………………………………… 6</p><p>　　3.1 仿真參數(shù)設(shè)置……………………………………………………………………………6</p><p>　　3.2 系統(tǒng)中仿真模塊參數(shù)的設(shè)置……………………………………………………………6</p><p>　　3.3 SCOPE端的最終波形圖…………………………………………………………………7&l

6、t;/p><p>　　3.4 調(diào)幅的包絡(luò)檢波和相干解調(diào)性能仿真比較……………………………………………8</p><p>　　3.5腳本程序………………………………………………………………………………… 9</p><p>　　4、總結(jié)與體會………………………………………………………………………………10</p><p>　　5、主要參考文獻(xiàn)………

7、……………………………………………………………………11</p><p><b>　　1引言</b></p><p>　　1.1 設(shè)計目的及任務(wù)要求</p><p>　　1．課程設(shè)計應(yīng)達(dá)到的目的</p><p>　?。?）掌握使用Matlab語言及其工具箱進(jìn)行基本信號分析與處理的方法。</p><p

8、>　　（2）用matlab和simulink設(shè)計一個通信系統(tǒng)，加深對通信原理基本原理和matlab應(yīng)用技術(shù)的理解；學(xué)習(xí)使用計算機(jī)建立通信系統(tǒng)仿真模型的基本方法及基本技能，學(xué)會利用仿真的手段對于實用通信系統(tǒng)的基本理論、基本算法進(jìn)行實際驗證；</p><p>　?。?）提高和挖掘?qū)W生將所學(xué)知識與實際應(yīng)用相結(jié)合的能力，學(xué)習(xí)現(xiàn)有流行通信系統(tǒng)仿真軟件MATLAB的基本使用方法，學(xué)會使用這些軟件解決實際系統(tǒng)出現(xiàn)的問題；

9、</p><p>　?。?）培養(yǎng)學(xué)生的合作精神和獨立分析問題和解決問題的能力；通過系統(tǒng)仿真加深對通信課程理論的理解。</p><p>　?。?）用MATLAB 完成調(diào)幅廣播系統(tǒng)的仿真，提高學(xué)生科技論文的寫作水平。</p><p>　　1.2 課程設(shè)計題目調(diào)幅廣播系統(tǒng)的仿真設(shè)計</p><p><b>　　設(shè)計任務(wù)：</b&g

10、t;</p><p>　　1.采用接收濾波器Analog Filter Design模塊，在同一示波器上觀察調(diào)幅信號在未加入噪聲和加入噪聲后經(jīng)過濾波器后的波形。采用另外兩個相同的接收濾波器模塊，分別對純信號和純噪聲濾波，利用統(tǒng)計模塊計算輸出信號功率和噪聲功率，繼而計算輸出信噪比，用Disply顯示結(jié)果。模型文件保存為ex1_1.mdl。</p><p>　　對中波調(diào)幅廣播傳輸系統(tǒng)進(jìn)行仿真，

11、其技術(shù)指標(biāo)為：</p><p>　　1)載波信號：幅度為1的正弦波，設(shè)初相為0，頻率在550~1605Hz內(nèi)可調(diào)；</p><p>　　2)基帶信號：調(diào)制度 (信號最大幅度與載波幅度之比) ma=0.3，頻率在100~600Hz內(nèi)可調(diào)；</p><p>　　3)接收機(jī)選頻濾波器帶寬為12KHz，中心頻率為1000KHz；</p><p>　　

12、4)在信道中加入加性高斯噪聲，需要先計算出信道中應(yīng)該加入噪聲的方差。設(shè)計接收機(jī)選頻濾波器輸出信噪比為20dB。</p><p>　　2.構(gòu)建包絡(luò)解調(diào)和相干解調(diào)電路，用示波器顯示解調(diào)波形。構(gòu)建一個信噪比測試子系統(tǒng)，該系統(tǒng)能使輸入的兩路解調(diào)信號中的信號和噪聲近似分離，以分別計算信號和噪聲分量的功率，進(jìn)而計算信噪比，并用Display顯示，同時將信噪比數(shù)據(jù)送入Workspace。模型文件保存為ex1_2.mdl。<

13、;/p><p>　　3.編寫腳本程序ex1.m，通過選擇不同信噪比，計算加性噪聲的方差送入仿真模型，調(diào)用模型文件執(zhí)行仿真，并通過matlab繪圖得到包絡(luò)解調(diào)和相干解調(diào)后的輸出信噪比與輸入信噪比的關(guān)系曲線。</p><p>　　2 調(diào)幅廣播系統(tǒng)的模型及仿真環(huán)境</p><p>　　MATLAB及Simulink建模環(huán)境簡介</p><p>　　M

14、ATLAB 是美國MathWorks公司出品的商業(yè)數(shù)學(xué)軟件，用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計算的高級技術(shù)計算語言和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分。</p><p>　　Simulink是MATLAB最重要的組件之一，它提供一個動態(tài)系統(tǒng)建模、仿真和綜合分析的集成環(huán)境。在該環(huán)境中，無需大量書寫程序，而只需要通過簡單直觀的鼠標(biāo)操作，就可構(gòu)造出復(fù)雜的系統(tǒng)。Simulink具有適應(yīng)面

15、廣、結(jié)構(gòu)和流程清晰及仿真精細(xì)、貼近實際、效率高、靈活等優(yōu)點，并基于以上優(yōu)點Simulink已被廣泛應(yīng)用于控制理論和數(shù)字信號處理的復(fù)雜仿真和設(shè)計。同時有大量的第三方軟件和硬件可應(yīng)用于或被要求應(yīng)用于Simulink。 </p><p>　　Simulink是MATLAB中的一種可視化仿真工具，是一種基于MATLAB的框圖設(shè)計環(huán)境，是實現(xiàn)動態(tài)系統(tǒng)建模、仿真和分析的一個軟件包，被廣泛應(yīng)用于線性系統(tǒng)、非線性系統(tǒng)、數(shù)字控制及

16、數(shù)字信號處理的建模和仿真中。Simulink可以用連續(xù)采樣時間、離散采樣時間或兩種混合的采樣時間進(jìn)行建模，它也支持多速率系統(tǒng)，也就是系統(tǒng)中的不同部分具有不同的采樣速率。為了創(chuàng)建動態(tài)系統(tǒng)模型，Simulink提供了一個建立模型方塊圖的圖形用戶接口(GUI) ，這個創(chuàng)建過程只需單擊和拖動鼠標(biāo)操作就能完成，它提供了一種更快捷、直接明了的方式，而且用戶可以立即看到系統(tǒng)的仿真結(jié)果.</p><p>　　Simulink&a

17、mp;reg;是用于動態(tài)系統(tǒng)和嵌入式系統(tǒng)的多領(lǐng)域仿真和基于模型的設(shè)計工具。對各種時變系統(tǒng)，包括通訊、控制、信號處理、視頻處理和圖像處理系統(tǒng)，Simulink提供了交互式圖形化環(huán)境和可定制模塊庫來對其進(jìn)行設(shè)計、仿真、執(zhí)行和測試。.</p><p><b>　　調(diào)幅廣播系統(tǒng)介紹</b></p><p>　　模擬幅度調(diào)制是無線電最早的遠(yuǎn)距離傳輸技術(shù)。在幅度調(diào)制中，以聲音信控

18、制高頻率正弦信號的幅度，并以幅度變化的高頻率正選信號放大后經(jīng)過天線發(fā)射出去，成為電磁波輻射。</p><p>　　波動的信號要能夠有效地從天線發(fā)送出去，或者有效地將信號接收過來，需要天線的長度至少達(dá)到波長的四分之一。聲音轉(zhuǎn)換成電信號后其波長為15～15000km之間，實際中不能造出這樣長度的天線進(jìn)行有效的信號收發(fā)。因此需要將象聲音信號這樣的低頻信號搬到較高的頻段上去，以便通過較短的天線發(fā)射出去。例如：移動通信所使

19、用的900MHz頻率段上的電磁波信號長度約為0.33米，其收發(fā)天線的尺寸應(yīng)為波長的四分之一，即約8cm左右。而調(diào)幅廣播中波范圍為550~1605kHz，短波約為3·30MHz，其波長范圍在幾十米到幾百米，相應(yīng)的天線要長一些。</p><p>　　人耳可聞的聲音信號通過話筒轉(zhuǎn)化為波動的電信號，其頻率范圍為20~20kHz。大量實驗發(fā)現(xiàn)，人耳對語音頻率敏感區(qū)域約為300~3400Hz,為了節(jié)約頻帶帶寬資源，

20、國際標(biāo)準(zhǔn)中將電話通信的傳輸頻帶規(guī)定為300~3400Hz。調(diào)幅廣播除了傳輸語音之外，還要播送音樂節(jié)目，這就需要更寬的頻帶。一般而言，調(diào)幅廣播的傳輸頻率范圍約為100~6000Hz。</p><p><b>　　模型參數(shù)指標(biāo)</b></p><p>　?。?）基帶信號：調(diào)制度 (信號最大幅度與載波幅度之比) ma=0.3，頻率在100~600Hz內(nèi)可調(diào)；</p&g

21、t;<p>　　（2）載波信號：幅度為1的正弦波，設(shè)初相為0，頻率在550~1605Hz內(nèi)可調(diào)；</p><p>　?。?）接收機(jī)選頻濾波器帶寬為12kHz，中心頻率為1000 kHz。</p><p>　?。?）信道中加入噪聲。當(dāng)調(diào)制度為0.3時，設(shè)計接收機(jī)選頻濾波器輸出信噪比為20dB。要求設(shè)計信道中應(yīng)加入噪聲的方差，并能夠測量接收選頻濾波器實際輸出信噪比。</p&

22、gt;<p><b>　　仿真參數(shù)設(shè)計</b></p><p>　　系統(tǒng)工作最高頻率為調(diào)幅載波頻率1605kHz，設(shè)計仿真采樣率為最高工作頻率的10倍左右，因此取仿真步長為</p><p>　　相應(yīng)的仿真帶寬為仿真采樣頻率的一半，即</p><p>　　設(shè)基帶測試正弦信號為 ，載波為，則調(diào)制度為為 的調(diào)制輸出信號 為</p&

23、gt;<p><b>　　顯然，的平均功率為</b></p><p>　　設(shè)信道無衰減，其中加入的白噪聲功率譜密度為，那么仿真帶寬內(nèi)噪聲的方差為</p><p>　　設(shè)接收選頻濾波器的功率增益為1，帶寬為B，擇選通濾波器的輸出噪聲功率為</p><p>　　因此，接收選通濾波器輸出信噪比為</p><p>

24、　　故信道中的噪聲方差為</p><p>　　根據(jù)上面的公式，編程計算出噪聲的方差，并將方差值和其它已知值作為仿真系統(tǒng)的參數(shù)。</p><p>　　根據(jù)仿真設(shè)計要求的輸出信噪比SNRout可計算出相應(yīng)信道中應(yīng)加入的噪聲方差值，計算程序和結(jié)果如下：</p><p>　　程序： SNR_dB=20;</p><p>　　SNR=10.^(

25、SNR_dB/10);</p><p><b>　　m_a=0.3;</b></p><p>　　P=0.5+(m_a^2)/4;</p><p>　　W=8025.7e3;</p><p><b>　　B=12e3;</b></p><p>　　sigma2=P/SNR*W

26、/B</p><p>　　運行后結(jié)果：sigma2 =3.4945</p><p>　　3. 系統(tǒng)的建立與仿真</p><p><b>　　仿真參數(shù)設(shè)置</b></p><p>　　按照調(diào)幅廣播系統(tǒng)的物理與數(shù)學(xué)模型建立系統(tǒng)模型。根據(jù)相干方式的原理圖利用MATLAB的Simulink建立系統(tǒng)的模擬仿真圖。如下圖所示：<

27、;/p><p>　　圖3-1-1　中波調(diào)幅廣播傳輸系統(tǒng)仿真參考模型</p><p>　　系統(tǒng)中仿真模塊參數(shù)的設(shè)置</p><p>　　Signal Generator: 信號發(fā)生器，產(chǎn)生基帶信號</p><p>　　Wave form: sine</p><p>　　Amplitude: 0.3</p>&l

28、t;p>　　Frequency: 1000</p><p>　　Signal Generator1: 信號發(fā)生器，產(chǎn)生載波信號</p><p>　　Wave form: sine</p><p>　　Amplitude: 1</p><p>　　Frequency: 1000000</p><p>　　Rando

29、m Number: 隨機(jī)噪聲發(fā)生器，產(chǎn)生高斯正態(tài)分布隨機(jī)信號，這里用來構(gòu)造高斯白噪聲信道</p><p><b>　　Mean: 0</b></p><p>　　Variance: 3.4945</p><p>　　Analog Filter Design: 模擬濾波器設(shè)計，三個模擬濾波器分別用于純信號，純噪聲以及信號和噪聲混合信號的濾波<

30、;/p><p>　　Design method: Butterworth</p><p>　　Filter type: Bandpass</p><p>　　Lower passband edge frequency (rads/sec): 2*pi*(1e6-6e3)</p><p>　　Upper passband edge frequenc

31、y (rads/sec): 2*pi*(1e6+6e3)</p><p>　　Zero-Order Hold: 零界保持器</p><p>　　Sample time: 6.23e-8</p><p>　　Variance: 計算向量的方差</p><p>　　選中Running variance</p><p>　　

32、dB Conversion: 分別對純信號和混合信號做對數(shù)變換</p><p>　　Convert to: dB </p><p>　　Input signal: Power </p><p><b>　　Fun: 運算函數(shù)</b></p><p>　　Expression: u(1)-u(2)</p>&

33、lt;p>　　Display: 顯示SNR的結(jié)果</p><p>　　Format: short</p><p>　　Scope 端的最終波形圖</p><p>　　在系統(tǒng)仿真模型圖中，用加法器和乘法器實現(xiàn)調(diào)幅，用 Random Number 產(chǎn)生噪聲樣值序列，并用加法器實現(xiàn)AWGN通道。為了測量輸出信噪比，以參數(shù)完全相同的三個濾波器模塊分別對純信號，純噪聲以

34、及信號和噪聲混合信號的濾波，最后利用統(tǒng)計模塊計算輸出信號功率和噪聲功率，繼而計算輸出信噪比。</p><p>　　某次仿真執(zhí)行后，測試信噪比為20.11dB，與設(shè)計值 20dB相符。按接收濾波器輸出的調(diào)幅信號以及發(fā)送調(diào)幅信號的波形對比仿真結(jié)果如下圖所示：</p><p>　　圖3-2-1 接收濾波器輸出調(diào)幅信號以及發(fā)送調(diào)幅信號的波形及仿真</p><p>　　3.

35、4 調(diào)幅的包絡(luò)檢波和相干解調(diào)性能仿真比較</p><p>　　實例2：以實例1為傳輸模型，在不同輸入信噪比條件下仿真測量包絡(luò)檢波解調(diào)和同步相干解調(diào)對調(diào)幅波的解調(diào)輸出信噪比，觀察包絡(luò)檢波解調(diào)的門限效應(yīng)。</p><p>　　圖1-2所示的仿真模型(ex2.mdl)用于測量包絡(luò)檢波的門限效應(yīng)，發(fā)送的調(diào)幅波參數(shù)以及仿真步進(jìn)與實例1相同。首先，調(diào)幅信號通過AWGN信道后，分別送入包絡(luò)檢波器和同

36、步相干解調(diào)器。包絡(luò)檢波器由Saturation模塊來模擬具有單向?qū)ㄐ阅艿臋z波二極管，模塊的上下門限分別設(shè)置為inf和0。兩解調(diào)器后接的低通濾波器相同。解調(diào)后的兩路信號送到示波器顯示，同時送入信噪比測試模塊，即圖中的子系統(tǒng)SNR Detection，其內(nèi)部如圖1-3所示。在SNR Detection模塊中，輸入的兩路解調(diào)信號通過濾波器將信號和噪聲近似分離，以分別計算信號和噪聲分量的功率，進(jìn)而計算信噪比。兩個帶通濾波器參數(shù)相同，其中心頻率

37、為1000Hz，帶寬為200Hz，對應(yīng)于發(fā)送基帶測試信號頻率，其輸出近似視為純信號分量。兩個帶阻濾波器參數(shù)也相同，其中心頻率為1000Hz，帶寬為200Hz，其輸出可近似為信號中的噪聲分量。之后，通過零階保持模塊將信號離散化，再由buffer模塊和方差模塊計算出信號和噪聲的功率，buffer緩沖區(qū)長設(shè)置為1.6051e+005個樣值，這樣將在0.01s內(nèi)進(jìn)行一次統(tǒng)計計算。最后，由分貝轉(zhuǎn)換模塊</p><p>　　

38、圖3-4-1 包絡(luò)檢波和相干解調(diào)性能測試仿真模型</p><p>　　系統(tǒng)中仿真模塊參數(shù)的設(shè)置</p><p>　　Signal Generator: 信號發(fā)生器，產(chǎn)生基帶信號</p><p>　　Wave form: sine</p><p>　　Amplitude: 0.3</p><p>　　Frequency:

39、 1000</p><p>　　Signal Generator1: 信號發(fā)生器，產(chǎn)生載波信號</p><p>　　Wave form: sine</p><p>　　Amplitude: 1</p><p>　　Frequency: 1000000</p><p>　　Random Number: 隨機(jī)噪聲發(fā)生器，產(chǎn)

40、生高斯正態(tài)分布隨機(jī)信號，這里用來構(gòu)造高斯白噪聲信道</p><p><b>　　Mean: 0</b></p><p>　　Variance:1</p><p>　　Analog Filter Design: </p><p>　　Design method: Butterworth</p><p&g

41、t;　　Filter type: Lowpass</p><p>　　Filter Order:2</p><p>　　passband edge frequency (rads/sec): 2*pi*6000</p><p>　　Display: 顯示SNR的結(jié)果</p><p>　　Format: short</p><

42、;p><b>　　Scope</b></p><p>　　Time range:0.01</p><p>　　圖3-4-2解調(diào)輸出信噪比近似于測量子系統(tǒng)SNR Detection的內(nèi)部結(jié)構(gòu)</p><p>　　系統(tǒng)中仿真模塊參數(shù)的設(shè)置</p><p>　　Analog Filter Design: </p&g

43、t;<p>　　Filter type: Bandstop</p><p>　　Lower passband edge frequency (rads/sec): 2*pi*900</p><p>　　Upper passband edge frequency (rads/sec): 2*pi*1100</p><p>　　Analog Filter

44、Design: </p><p>　　Filter type: Bandpass</p><p>　　Lower passband edge frequency (rads/sec): 2*pi*900</p><p>　　Upper passband edge frequency (rads/sec): 2*pi*1100</p><p>

45、　　Zero-Order Hold: 零界保持器</p><p>　　Sample time: 6.23e-8</p><p>　　Variance: 計算向量的方差</p><p>　　選中Running variance</p><p><b>　　Buffer</b></p><p>　　O

46、utput buffersize:1.6051e+005</p><p>　　dB Conversion: 分別對純信號和混合信號做對數(shù)變換</p><p>　　Convert to: dB </p><p>　　Input signal: Power </p><p><b>　　Fun: 運算函數(shù)</b></p

47、><p>　　Expression: u(3)-u(1)</p><p>　　Expression: u(4)-u(2)</p><p>　　Display: 顯示SNR的結(jié)果</p><p>　　Format: short</p><p>　　圖3-4-3 執(zhí)行仿真所得到的解調(diào)信號波形噪聲方差設(shè)為1</p>

48、<p>　　3.5任務(wù)三matlab源程序ex1.m</p><p><b>　　clc;</b></p><p>　　clear all;</p><p>　　SNR_in_dB=0:2:6;</p><p>　　SNR_in=10.^(SNR_in_dB./10); % 信道信噪比</p>

49、<p>　　m_a=0.3; % 調(diào)制度</p><p>　　P=0.5+(m_a^2)/4; % 信號功率</p><p>　　for k=1:length(SNR_in)</p><p>　　sigma2=P/SNR_in(k);</p><p>　　% 計算信道噪聲

50、方差并送入仿真模型</p><p>　　sim('ex12.mdl');% 執(zhí)行仿真</p><p>　　SNRdemod(k,:)=SNR_out(:,1); % 記錄仿真結(jié)果</p><p><b>　　end</b></p><p>　　plot(SNR_in_dB, SNRdemod(:,size

51、(SNRdemod)));</p><p>　　xlabel('輸入信噪比 dB');</p><p>　　ylabel('解調(diào)輸出信噪比 dB');</p><p>　　legend('包絡(luò)檢波','相干解調(diào)');</p><p>　　將圖3—4—1中的Random Numbe

52、r: 設(shè)置為</p><p><b>　　Mean: 0</b></p><p>　　Variance: sigma2</p><p>　　通過matlab繪圖得到包絡(luò)解調(diào)和相干解調(diào)后的輸出信噪比與輸入信噪比的關(guān)系曲線：</p><p><b>　　4. 總結(jié)與體會</b></p>&

53、lt;p>　　通過這次的課程設(shè)計，我們對信息和通信系統(tǒng)有了更進(jìn)一步的認(rèn)識，尤其是在系統(tǒng)設(shè)計方面，盡管是非?；A(chǔ)的調(diào)幅廣播從系統(tǒng)的仿真，也是經(jīng)過若干設(shè)備協(xié)同工作，才能保證信號有效傳輸，而小到僅僅是一個參數(shù)，都有可能導(dǎo)致整個系統(tǒng)無法正常運行。兩周的課程設(shè)計，讓我們領(lǐng)教了MATLAB矩陣實驗室強(qiáng)大的功能和實力。通過在Simulink環(huán)境下對系統(tǒng)進(jìn)行模塊化設(shè)計與仿真，使我們獲得兩方面具體經(jīng)驗，第一是MATLAB中Simulink功能模塊的

54、使用方法，第二是圖形化和結(jié)構(gòu)化的系統(tǒng)設(shè)計方法。</p><p>　　在整個課程設(shè)計過程中也遇到很多現(xiàn)實的問題，比如各版本MATLAB軟件并不完全兼容，許多復(fù)雜模塊參數(shù)深奧難以正確設(shè)置，，比如參數(shù)設(shè)置的不理想因此總是會出現(xiàn)波形失真的現(xiàn)象等問題。但是通過上網(wǎng)查找資料和查詢參考書能夠讓我更好的完成此次設(shè)計。同時這次設(shè)計也讓我能夠更好的對應(yīng)用工具M(jìn)ATLAB有一個進(jìn)一步的了解和應(yīng)用。</p><p&g

55、t;　　在學(xué)習(xí)MATLAB理論基礎(chǔ)后，我們又在此基礎(chǔ)上通過利用MATLAB仿真真正的看到了通信中傳輸信息的一系列的問題。比如說要使信號不失真的能夠傳輸?shù)浇邮斩司鸵紤]很多的因數(shù)。在發(fā)送端要注意噪聲的加入，盡量的減少噪聲進(jìn)入信道中，以免在接收端使信號失真度過大而不能夠恢復(fù)成原來的信號。而在接收端，采用哪種解調(diào)方式能夠更好的恢復(fù)出原來的信號。</p><p><b>　　5、主要參考文獻(xiàn)</b>

56、</p><p><b>　　參考資料：</b></p><p>　　【1】 邵玉斌. Matlab/Simulink通信系統(tǒng)建模與仿真實例分析. 北京:清華大學(xué)出版社, 2008</p><p>　　【2】張化光, 劉鑫蕊, 孫秋野. MATLAB/SIMULINK實用教程. 北京:人民郵電出版社, 2009</p><p