24小时热门版块排行榜    

查看: 717  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

bohrzqc

新虫 (正式写手)


[交流] 【求助】dde23方程求解的数据调用问题

求教一个dde23数据调用问题。问题的背景是这样的。要实现一个通信,发射机用solt=dde23解出的方程解做发射机的载波,然后数据yt(2,和yt(3,要经过信道传输,传输后变为yr2和yr3,这两组数据要送入接收机,即solr=dde23中求解。
此时遇到了一个问题,如何将yr2和yr3的数据给接收机呢?

源程序如下:
CODE:
solt=dde23(@sctrt...
   t=1/fs:1/fs:tf;   
[yt,ypt]=deval(solt,t);

e=sqrt(yt(2,:));
  pht=yt(3,:);
  u=e.*exp(i*pht);

p.yr2=fibE.^2;
p.yr3=angle(f);  

solr=dde23(@sctrr...   

子程序是
CODE:
function  xdot=sctrr(t,x,z,p)
                 
xlag=z;  

xt(2)????
xt(3)????如何编程

[ Last edited by xiegangmai on 2011-3-9 at 21:20 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bohrzqc

新虫 (正式写手)



ben_ladeng: 金币+1, 欢迎常来 2013-04-10 21:14:27
引用回帖:
Originally posted by 信彼南山 at 2011-03-07 17:34:10:
fs和tf没定义
子程序里面啥都没有

这叫什么源程序啊?



fs是采样率
tf是结束时间

下面的是子程序
function  xdot=sctrr(t,x,z,p,st,stt)
                 
xlag=z;  

kcf=p.kfr/p.tl;
kcj=p.kj/p.tl;


if t-p.tj>0      

xt(2)=
  xt(3)=

%
thr=p.wr*p.ter+p.th_e+x(3)-xlag(3);
thj=p.wt*p.tj+x(3)-xt(3)+(p.wt-p.wr)*t;
      dnr=p.ir/(p.V*p.q)-1/p.tn*x(1)-x(2)*p.Gn*(x(1)-p.n0)*(1+p.eps*x(2))^-1;
      dsr=1/p.tn*p.gm*p.yt*x(1)-1/p.tpr*x(2)+x(2)*p.Gn*(x(1)-p.n0)*((1+p.eps*x(2))^-1)*p.yt...
                                  +2*kcf*sqrt(xlag(2)*x(2))*cos(thr)+2*kcj*sqrt(xt(2)*x(2))*cos(thj);
      dfr=1/2*p.af*(p.Gn*(x(1)-p.n0)*((1+p.eps*x(2))^-1)*p.yt-1/p.tpr)...
          -kcf*sqrt(xlag(2)/x(2))*sin(thr)-kcj*sqrt(xt(2)/x(2))*sin(thj);

xdot=[dnr;
       dsr;
       dfr];
3楼2011-03-07 18:52:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

信彼南山

木虫 (著名写手)


bohrzqc(金币+1): 2011-03-07 18:50:02
fs和tf没定义
子程序里面啥都没有

这叫什么源程序啊?
2楼2011-03-07 17:34:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见