| 查看: 2893 | 回复: 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矩阵维度不匹配，求高人的指点如何修改？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 矩阵元素这个循环怎么循环，有点晕
已经有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
