| 24小時(shí)熱門(mén)版塊排行榜 |
| 6 | 1/1 | 返回列表 |
| 查看: 910 | 回復(fù): 5 | ||
xzfgggthn新蟲(chóng) (小有名氣)
|
[求助]
在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; 然后主程序里我首先寫(xiě)上了那個(gè)ask信號(hào)的代碼,就是我帖子里最開(kāi)始那段代碼,然后寫(xiě)上了[t,y]=ode45(@odefun,[0 0.5],[1 1]); 請(qǐng)問(wèn)這樣寫(xiě)哪里有問(wèn)題,為啥程序總報(bào)錯(cuò),感覺(jué)就是方程里的out有問(wèn)題,但是不知道怎么改 |
銀蟲(chóng) (著名寫(xiě)手)
新蟲(chóng) (小有名氣)
![]() ![]() ![]() ![]() |
新蟲(chóng) (初入文壇)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程專碩調(diào)劑 +3 | 204818@lcx 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +4 | yangfz 2026-03-17 | 4/200 |
|
|
[考研] 268求調(diào)劑 +6 | 好運(yùn)連綿不絕 2026-03-12 | 7/350 |
|
|
[考研] 332求調(diào)劑 +5 | Zz版 2026-03-13 | 5/250 |
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 5/250 |
|
|
[考研] 085601求調(diào)劑 +3 | Du.11 2026-03-16 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 快樂(lè)的香蕉 2026-03-11 | 4/200 |
|
|
[考研] 308求調(diào)劑 +5 | 是Lupa啊 2026-03-11 | 5/250 |
|
|
[考研] 求材料調(diào)劑 +5 | 隔壁陳先生 2026-03-12 | 5/250 |
|
|
[考研] 四川大學(xué)085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 274求調(diào)劑0856材料化工 +12 | z2839474511 2026-03-11 | 13/650 |
|
|
[考研] 283求調(diào)劑,材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|