| 5 | 1/1 | 返回列表 |
| 查看: 915 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 085601英二數(shù)二求調(diào)劑 總分325 +3 | 余航航 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +6 | Faune 2026-03-30 | 6/300 |
|
|
[考研] 一志愿浙江大學(xué)工科動(dòng)力工程370,數(shù)一121,專業(yè)課135,現(xiàn)在能去哪里 +3 | 080700調(diào)劑 2026-03-30 | 4/200 |
|
|
[考研] 289求調(diào)劑 +6 | BrightLL 2026-03-29 | 6/300 |
|
|
[考研] 266分,求材料相關(guān)專業(yè)調(diào)劑 +10 | 哇呼哼呼哼 2026-03-30 | 12/600 |
|
|
[考研] 291求調(diào)劑 +12 | Y-cap 2026-03-29 | 15/750 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +9 | @taotao 2026-03-29 | 9/450 |
|
|
[考研] 285求調(diào)劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +10 | Nadiums 2026-03-25 | 11/550 |
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研] 一志愿上海理工能源動(dòng)力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|