| 查看: 919 | 回复: 5 | ||
[求助]
在matlab里用ode45函数，微分方程包括ask信号怎么办 已有2人参与
|
|
大概的意思就是我想用obe45解一个微分方程组，方程里的一项S（t）为ask信号 首先我的输入信号是一个ask信号，代码如下： %=============================1、构造载波，产生8个码元，生成已调信号========% a=randsrc(1,8,[0:1]);%产生8个随机的二进制数 l=linspace(0,2*pi,50);%利用linspace函数创建数组,2pi长度取点50个模拟一个码元 f=sin(2*l);%生成载波 t=linspace(0,10*pi,400);%定义时轴length为10pi,取点400个，代表8个码元的总取样点数 out=1:400;%规定已调信号length b=1:400;% 规定基带信号length w=1:400;%规定载波length %=================生成ASK信号===========% 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作为调制信号输出 w(j+50*(i-1))=f(j); %w作为载波输出 end end （这个ask信号代码是没问题的，我单独运行过，可以生成ask信号） 然后我要把这个ask的输出信号作为微分方程的一项，我觉得他的输出信号应该是out。 于是我创建了个m文件，定义了一个微分方程，微分方程里要包括我刚才的ask信号，代码是 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; 然后主程序里我首先写上了那个ask信号的代码，就是我帖子里最开始那段代码，然后写上了[t,y]=ode45(@odefun,[0 0.5],[1 1]); 请问这样写哪里有问题，为啥程序总报错，感觉就是方程里的out有问题，但是不知道怎么改 |
回复此楼» 猜你喜欢
309分085801求调剂
已经有7人回复
-
生物学308分求调剂（一志愿华东师大）
已经有3人回复
-
求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过
已经有3人回复
-
环境工程 085701，267求调剂
已经有14人回复
-
308求调剂
已经有12人回复
-
求调剂：085600材料与化工，考材科基，总分319
已经有21人回复
-
311求调剂
已经有8人回复
-
294分080500材料科学与工程求调剂
已经有13人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有19人回复
-
288资源与环境专硕求调剂，不限专业，有学上就行
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB求解数值微分方程的问题
已经有3人回复
- matlab 如何求解多元二阶微分方程组
已经有7人回复
- matlab 运行ode45出错，不知道什么原因
已经有4人回复
- 三元二次的三角方程如何用MATLAB求解？？
已经有10人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- 如何用Matlab求时变微分方程的解？
已经有5人回复
- matlab ode15s方程中两个参数怎么调用其他函数带入求解
已经有6人回复
- matlab解偏微分方程求助
已经有12人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab解微分方程
已经有10人回复
- matlab-常微分方程参数估计
已经有12人回复
- matlab解微分方程组
已经有15人回复
- matlab 计算过程中工作空间变量保存及读取 求解决内存不足的方法
已经有5人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】求助用matlab解偏微分方程
已经有15人回复
2楼2015-04-07 06:56:02
3楼2015-04-07 11:17:25
xzfgggthn
兑换贵宾
 |
4楼2015-04-08 12:03:15
a11835
兑换贵宾
- 应助: 1 (幼儿园)
- 金币: 48
- 散金: 2269
- 红花: 3
- 帖子: 142
- 在线: 36.3小时
- 虫号: 2822898
- 注册: 2013-11-24
- 性别: GG
- 专业: 通信理论与系统
5楼2015-04-25 20:27:21
6楼2018-12-04 18:54:34
