作者hopeandhmm (......)
看板MATLAB
标题[其他] 除法与函数问题
时间Tue Apr 24 16:04:19 2012
请教一下大家
我想了很久
还是不知道哪个地方出错了
先附上我的程式码
x1=linspace(3.4*10^12,3.6*10^12,1024); %%频率轴
sig1=sqrt(1000); %%高斯函数之参数
u1=0; %%高斯函数之参数
y1= 1/(sig1*sqrt(2*pi))*exp(-((((x1/7/10^12)*10-5)-u1).^2/2*sig1^2));%%高斯函数
A=2.*real(y1.*exp(i.*x1.*5));%%高斯函数乘以一个延迟相位的实部之後再乘以两倍
B=2.*imag(y1.*exp(i.*x1.*5));%%高斯函数乘以一个延迟相位的虚部之後再乘以两倍
for v=1:1024
C(v)=B(v)/A(v); %%B阵列除以A阵列
end
theta=atan(C); %%希望获得相位值x1.*5
plot(x1,theta,'-g'),xlabel('Frequency'),ylabel('Phase'),title('Spectral interferograms');
axis tight
我想请教的第一个问题是B阵列除以A阵列的数值
和实际上的数值有一些差距
这个要怎麽处理呢
第二个问题是
相位值theta对频率作图应该是一条直线
因为相位值是x1.*5
而频率值是x1
但结果却不是
想请教要怎麽解决这个问题呢
谢谢大家
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.116.92.238
1F：→ vincent59:实际上的数值有一些差距是指???? 04/25 23:23
2F：→ hopeandhmm:与实际上述直有些差距，是指跟用计算机计算的精确结果 04/26 00:08
3F：→ hopeandhmm:有差异 04/26 00:08