a=[0 0 0 1 1 0 1 0 1 1];
K=length(a)
a3=a;
a2=a;
for j=[2:1:K]
if a(j)==a(j-1)
a2(j)=0;
else
a2(j)=1;
end
end
a
a=a2;
a
initial_phase=pi;
f1=1;
N=200;
i=[0:1:N-1];
sin2=sin(2*pi*f1*i/N);
sin1=sin(2*pi*f1*i/N+pi);
for j=[1:1:K]
for i=[1:1:N]
yout(N*(j-1)+i)=a(j)*sin1(i)+(1-a(j))*sin2(i);
end
end
sin22=sin2;
sin11=sin1;
for j=[1:1:K-1]
sin22=[sin22 sin2];
sin11=[sin11 sin1];
end
a=a3;
for j=[1:1:K]
for i=[1:1:N]
a1(N*(j-1)+i)=a(j);
end
end
figure(1)
plot(a1,'--r')
hold on
plot(yout)
grid on

mgs_sgm
- 粉丝: 0
最新资源
- 软件测试工程师的疑惑.doc
- 供电系统安全管理.doc
- 烟草专卖局(公司)年度信息系统安全检查工作自查报告.doc
- 网络舆情管理信息系统技术方案.doc
- 电子商务专业英语词汇表.doc
- assembly_learning-汇编语言资源
- 网站需求说明书软件工程课程设计.doc
- 服务热线网络管理平台(DOC页).docx
- 动态规划算法原理与的应用.doc
- 2023年MSoffice计算机二级考点.docx
- 用友软件食品行业烘焙细分行业ERP信息化解决方案.doc
- 网络营销调研培训教材.pptx
- 信息安全与计算机病毒的防范教材.pptx
- 供应链网络设计.ppt
- 自考数据库系统原理04735真题模拟含答案.doc
- 北京交通大学微机原理与接口技术作业答案.docx
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈


