24小时热门版块排行榜

返回列表

csuwt

木虫 (著名写手)

应助: 6 (幼儿园)
金币: 4738.3
散金: 1
红花: 1
帖子: 1189
在线: 95.8小时
虫号: 645925
注册: 2008-11-04
专业: 创新管理

[交流] 求解虚数相位奇变问题

要计算一个系统相频特性，就要用到反正切函数，计算机中反正切函数规定，在一、二象限中的角度为0～pi，三四象限的角度为0～-pi。
若一个角度从0变到2pi，但实际得到的结果是0～pi，再由-pi～0，在w=pi处发生跳变，跳变幅度为2pi，这就叫相位的卷绕。
unwrap（w）就是解卷绕，使相位在pi处不发生跳变，从而反应出真实的相位变化。
----------------------------------------------------------------
%Prepare the data, coyboy 2010-4-1
a = -150:200; b = -53 : 297; y = a + b*i;
% 图1，使用函数angle
y1 = angle(y);
figure(1); subplot(1,3,1);plot(y1);title('angle');
% 图2，单纯使用 atan，也就是楼主的做法
fai = atan(imag(y)./real(y));
subplot(1,3,2); plot(fai);title('atan');
% 图3，使用unwrap 和 atan，注意，我这边乘以10再除以10，为了骗过unwrap函数
y2 = unwrap(10*atan(imag(y)./real(y)),3)/10;
subplot(1,3,3); plot(y2);title('unwrap(atan)');

回复此楼

» 猜你喜欢

1楼 2011-12-19 22:55:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 csuwt 的主题更新

返回列表