| 5 | 1/1 | 返回列表 |
| 查看: 914 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
xzfgggthn新蟲 (小有名氣)
|
[求助]
在matlab里用ode45函數(shù),微分方程包括ask信號(hào)怎么辦 已有2人參與
|
|
|
大概的意思就是我想用obe45解一個(gè)微分方程組,方程里的一項(xiàng)S(t)為ask信號(hào) 首先我的輸入信號(hào)是一個(gè)ask信號(hào),代碼如下: %=============================1、構(gòu)造載波,產(chǎn)生8個(gè)碼元,生成已調(diào)信號(hào)========% a=randsrc(1,8,[0:1]);%產(chǎn)生8個(gè)隨機(jī)的二進(jìn)制數(shù) l=linspace(0,2*pi,50);%利用linspace函數(shù)創(chuàng)建數(shù)組,2pi長(zhǎng)度取點(diǎn)50個(gè)模擬一個(gè)碼元 f=sin(2*l);%生成載波 t=linspace(0,10*pi,400);%定義時(shí)軸length為10pi,取點(diǎn)400個(gè),代表8個(gè)碼元的總?cè)狱c(diǎn)數(shù) out=1:400;%規(guī)定已調(diào)信號(hào)length b=1:400;% 規(guī)定基帶信號(hào)length w=1:400;%規(guī)定載波length %=================生成ASK信號(hào)===========% for i=1:8 if a(i)==0 for j=1:50 out(j+50*(i-1))=0; %若碼元為0則將0輸出 end else for j=1:50 out(j+50*(i-1))=5*f(j); %若碼元為1則將載波輸出 end end end for i=1:8 for j=1:50 b(j+50*(i-1))=a(i); %b作為調(diào)制信號(hào)輸出 w(j+50*(i-1))=f(j); %w作為載波輸出 end end (這個(gè)ask信號(hào)代碼是沒(méi)問(wèn)題的,我單獨(dú)運(yùn)行過(guò),可以生成ask信號(hào)) 然后我要把這個(gè)ask的輸出信號(hào)作為微分方程的一項(xiàng),我覺(jué)得他的輸出信號(hào)應(yīng)該是out。 于是我創(chuàng)建了個(gè)m文件,定義了一個(gè)微分方程,微分方程里要包括我剛才的ask信號(hào),代碼是 function dy=odefun(t,y)dy=zeros(2,1); dy(1)=y(2); dy(2)=y(2)+y(1)-y(1)^3+cos(w)+out; 然后主程序里我首先寫上了那個(gè)ask信號(hào)的代碼,就是我帖子里最開始那段代碼,然后寫上了[t,y]=ode45(@odefun,[0 0.5],[1 1]); 請(qǐng)問(wèn)這樣寫哪里有問(wèn)題,為啥程序總報(bào)錯(cuò),感覺(jué)就是方程里的out有問(wèn)題,但是不知道怎么改 |
新蟲 (小有名氣)
銀蟲 (著名寫手)
![]() ![]() ![]() ![]() |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +12 | 辰沐5211314 2026-03-26 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)材料與化工(085600)296求調(diào)劑 +24 | 稻妻小編 2026-03-26 | 24/1200 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
[考研] 329求調(diào)劑,一志愿西北工業(yè)大學(xué),材料工程(085601) +5 | 小小機(jī)靈蟲 2026-03-29 | 11/550 |
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +6 | 深V宿舍吧 2026-03-29 | 6/300 |
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑。。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-25 | 7/350 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|