How do I get a outer loop in a polar plot in MATLAB

2017 11 6
1 个回答
回答已采纳
更新时间：2017 11 8
12 次查看（30 天）

0 个投票

To display the multiple microphone directive pattern, I wrote a simple MATLAB program.
% Super-cardiod = 0.57
% Hyper-Cardiod = 0.33
A_sen = 20.0; % sensitivity
c = 343; % speed of sound 343 m/s
Beta_1 = 0.33; % beta value for the shape of the directivity pattern
Beta_2 = 0.57; % beta value for the shape of the directivity pattern
d = 0.09; % distance between two microphones
Ang_freq = 90; % Angular frequency of sound
theta = linspace(-2*pi,2*pi); % Angular ref to microphone
U_total_1 = (A_sen *(Ang_freq*d)/c)*(Beta_1+cos(theta));
U_total_2 = (A_sen *(Ang_freq*d)/c)*(Beta_2+cos(theta));
subplot(2,1,1),polar(theta,U_total_1)
subplot(2,1,2),polar(theta,U_total_2)
The MATLAB program works to display the directivity pattern as below,
However, I'd like the inner loop to be an outer loop as shown in the link below.
https://music.stackexchange.com/questions/29895/microphone-pick-up-patterns-for-live-performance-advantages-disadvantages-of-c
I suppose this is due to the axes limits. But I'm not sure. How do I get the inner loop to show as an outer loop. Thank you..
Note: I am using R2014b. Therefore I don't have the polarplot function. So I will have to stick to the polar function in MATLAB.

4 个评论
显示 2更早的评论 隐藏 2更早的评论

KSSV
KSSV 2017-11-6
But there are no loops in your code..?
Big heads
Big heads 2017-11-6
If you see the figure or if you run the code, you will see that there is a big cardioid curve (big loop) and then there is a smaller loop (smaller cardioid) within the bigger cardioid. There is no loop/curve in the code. The loop/curve is in the polar plot.
Walter Roberson
Walter Roberson 2017-11-6
You are trying to create a hypercardioid ?

请先登录，再进行评论。

请先登录，再回答此问题。

 采纳的回答

David Goodmanson
David Goodmanson 2017-11-8

0 个投票

Hi Union Olive,
If you replace
U_total_1 = (A_sen *(Ang_freq*d)/c)*(Beta_1+cos(theta));
with
U_total_1 = abs(A_sen *(Ang_freq*d)/c)*(Beta_1+cos(theta)));
and similarly with U_total_2 it should work.

4 个评论
显示 2更早的评论 隐藏 2更早的评论

Big heads
Big heads 2017-11-8
编辑：Big heads 2017-11-8
Thank you. That worked. I did.
subplot(2,1,1),polar(theta,abs(U_total_1))
If I can get greedy here, do you know if there is another way to change the axes limits, I'd like to be able to display the negative values as well. I know it is possible using the new 'polarplot' function, but is there a way to do so with the 'polar' function?**
David Goodmanson
David Goodmanson 2017-11-8
I had some doubts about how that pattern would look so I ran the hypercardioid case using polarplot and the limits [-.4, .7]. It doesn't look very helpful, since now U_total_1 = 0 is not at the origin any more and at first glance the pattern looks like a subcardioid. You lose a lot by not having zero output at the origin, and I believe the abs option is the most physically meaningful plot (although I have not looked into the possibility that abs^2 is what is wanted).
I would post a picture of the plot except I can't figure out how to copy it into this window, although almost everyone else in Matlabland can do it. (I know I could attach it as a png file or whatever but I would like to have the window option as well).
Walter Roberson
Walter Roberson 2017-11-8
The picture-Frame button between the link icon and the paper clip is to insert images. Just make sure you have an empty line before and after the coding of the link. The system uses <<URL>> as the indication to put the image in-line.
Big heads
Big heads 2017-11-8
Thanks David. I'll stick with the abs option. Since this is just for simulation, I'm fine with 'abs', I just wanted to see if there was anything else that could be done. Thank you again.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Polar Plots 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by