clear?all
clc
N?=?1;
M?=?10;
j?=?sqrt（-1）;
theta1?=?（-40*pi）/180;
theta2?=?（40*pi）/180;
theta?=?linspace（theta1theta21000）;
atheta0?=?ones（M1）;
for?k=1:1000
????
????????atheta?=?zeros（M1）;??????%導(dǎo)引矢量
????????????????for?i=1:M
????????????????????atheta（i1）=exp（-j*（i-1）*pi*sin（theta（k）））;
????????????????end
????????b?=?ones（M1）;
????????c1?=?（N/M）*（abs（atheta0‘*b）.^2）*（abs（atheta‘*atheta0）.^2）;????????%coherent?transmit?signals收發(fā)形式
????????d1?=?（N/（M））*（abs（atheta‘*atheta0）.^4）;???????????????????????????????%orthogonal?transmit?signals收發(fā)形式
????????c1=c1/1000;
????????d1=d1/1000;
????????c（1k）?=?10*log10（c1）;
????????d（1k）?=?10*log10（d1）;
end
theta?=?theta*（180/pi）;
figure
plot（thetac‘r--‘thetad‘b‘）
xlabel（‘θ[deg]‘）ylabel（‘G（θ）[db]‘）
legend（‘β=1‘‘β=0‘4）
title（‘Beam?pattern?for?orthogonal?and?coherent?signals‘）

savepath1=‘E:\matlab\fangzhen?jieguo‘;
savename=‘firstfigure‘;%圖像的文件名

?%print（2‘-djpeg‘strcat（savepath1‘\‘savename‘.jpeg‘））;?%將句柄為2的圖形保存為jpeg/jpg格式的圖片，文件名為‘C:\abc.jpeg‘。
saveas（gcfsavename‘fig‘）;

?
set（gcf‘visible‘‘off‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????1158??2009-04-17?17:16??beampattern.m

-----------?---------??----------?-----??----

?????????????????1158????????????????????1