| 查看: 2577 | 回复: 5 | ||
[求助]
怎么将此for循环改为矩阵形式？
|
|
M=2; % M进制调制 BW=25e+3; delta_f=1e+3; BW1=(BW-2*delta_f)/2 Td=1/BW1 % 码元周期 fd=1/Td fs=15.5e6; % 采样频率 f1=1e+6; f2=f1+delta_f; BW=2*BW1+2*delta_f fc=[f1,f2]; % 载频 counter=Td*fs; T=20e-3; % 模拟信号长度 N=T*fs % 总的样本点数 Ns=floor(N/counter) % 码元数 code0=rand(1,Ns); code=fix(code0*M); % 产生随机码元 step = 1/fs; t= 0: step: Td*Ns-step; length(t) for m=1:length(t) a=floor(t(m)/(Td)); y(m)=sin(2*pi*fc(code(a+1)+1)*t(m)); phase_shift=exp(-i*2*pi*fc(code(a+1)+1)); end; % m=[1:length(t)]; % y(m)=sin(2*pi*fc(code(floor(t(m)/(Td))+1)+1)*t(m)); % phase_shift=exp(-i*2*pi*fc(code(floor(t(m)/(Td))+1)+1)); 程序如上，可正常运行，但考虑到for循环耗时太久，所以讲for循环改为下面注释的部分，改之后提示inner矩阵维度不匹配，求高人的指点如何修改？ |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 矩阵元素这个循环怎么循环，有点晕
已经有56人回复
- 如何用for或while循环给矩阵赋值
已经有3人回复
- matlab如何实现矩阵多行数据的循环调用
已经有5人回复
- matlab中如何通过for语句得到已有矩阵的部分数据
已经有3人回复
- 用fortran存储稀疏矩阵
已经有4人回复
- fortran能求带字母矩阵的本征值吗？（矩阵元素带有字母，能得到形式解吗）
已经有5人回复
- 【分享】Direct MEthods for Sparse Linear Systems （稀疏矩阵，基本算法）
已经有76人回复
- 怎样将for循环的结果存入数组
已经有6人回复
- fortran程序怎样调用lapack库进行矩阵对角化？
已经有7人回复
- 【求助】fortran求解矩阵
已经有7人回复
- 【求助】fortran：如何求解 复数矩阵的秩【已完结】
已经有5人回复
- 【求助】用fortran解矩阵问题【已解决】
已经有5人回复
zkw666666
木虫 (正式写手)
益虫。。。。
- 应助: 5 (幼儿园)
- 贵宾: 1
- 金币: 1728.3
- 红花: 1
- 帖子: 732
- 在线: 230.7小时
- 虫号: 21704
- 注册: 2003-08-27
- 性别: GG
- 专业: 计算机应用技术
2楼2012-11-25 22:07:11
zkw666666
木虫 (正式写手)
益虫。。。。
- 应助: 5 (幼儿园)
- 贵宾: 1
- 金币: 1728.3
- 红花: 1
- 帖子: 732
- 在线: 230.7小时
- 虫号: 21704
- 注册: 2003-08-27
- 性别: GG
- 专业: 计算机应用技术
3楼2012-11-25 22:09:25
mei3014
金虫 (正式写手)
木有头屑
- 应助: 14 (小学生)
- 金币: 2957.5
- 散金: 274
- 帖子: 541
- 在线: 335小时
- 虫号: 1333109
- 注册: 2011-06-28
- 性别: GG
- 专业: 信号理论与信号处理
4楼2012-11-27 12:01:05
mei3014
金虫 (正式写手)
木有头屑
- 应助: 14 (小学生)
- 金币: 2957.5
- 散金: 274
- 帖子: 541
- 在线: 335小时
- 虫号: 1333109
- 注册: 2011-06-28
- 性别: GG
- 专业: 信号理论与信号处理
【答案】应助回帖
★
xiegangmai: 金币+1, 辛苦了！ 2012-11-27 22:51:07
xiegangmai: 金币+1, 辛苦了！ 2012-11-27 22:51:07
|
% for m=1:length(t) % a=floor(t(m)/(Td)); % y(m)=sin(2*pi*fc(code(a+1)+1)*t(m)); % phase_shift=exp(-i*2*pi*fc(code(a+1)+1)); % end; m=[1:length(t)]; y(m)=sin(2*pi*fc(code(floor(t(m)/(Td))+1)+1).*t(m)); phase_shift=exp(-i*2*pi*fc(code(floor(t(m)/(Td))+1)+1)); |
» 本帖已获得的红花（最新10朵）
化名为哀5楼2012-11-27 12:02:04
6楼2012-11-27 14:10:47
