辐射场指向性分析（Directivity）产生的问题

榫眼

问题：如何用Matlab绘出sysnoise中的辐射场指向性图（Directivity）
Sysnoise中的步骤：计算完模型响应后
1.     Analysis\Directivity
2.    postprocess\directivity diagram  如图1、图2
3.  提取数据 tools\combine 图3  图4
目的：根据第3步提取的数据，用Matlab绘出图2的图形，
我以图4中，A列代表角度，X实部数据画极轴图
却始终绘不出图2那样的图形，不知问题在哪，是否牵涉到公式什么的？谢谢!






[ 本帖最后由 榫眼 于 2009-3-10 19:43 编辑 ]

w89986581

呵呵，肯定是可以的，记得取幅度。

榫眼

谢谢主任的回答！
我想问的是幅度指的是幅值吗？实部与虚部的平方和开根号吗？
这样好像也不对

w89986581

提供一个sysnoise输出结果转存程序(matlab)。计算完指向性以后，自动生成场点，重新计算场点结果（可以多个频率），输出result（格式为txt），然后调用下面这个函数。
仅供参考，自己多动手。

function [] = loadsysnoise(filename)
%%%读取sysnoise计算的指向性结果
% filename = ['C:\MATLAB6p5\work\cylindercone\directivitytest4.txt'];
fid  = fopen(filename,'r');
if fid == -1,
    disp('Invalid datafile')
    return
end
while 1,
    dataline1 = fgetl(fid);
    if ~ischar(dataline1),
        fclose(fid);
        break;
    end
    if isequal(dataline1, 'FREQUENCY'),        
        dataline1 = fgetl(fid);
        filename1 = ['C:\MATLAB6p5\work\cylindercone\directivitytest4',num2str(str2num(dataline1)),'Hz.txt'];
        fid1  = fopen(filename1,'w');
    end   
    if isequal(dataline1, 'PRESSURE VALUES'),
        nodenum = 1;         
        while 1,         
            dataline1 = fgetl(fid);
            if isequal(dataline1, 'VELOCITY VALUES'),
                fclose(fid1);
                break
            end   
            fprintf(fid1,'%s \n',dataline1);
            dataline1 = fgetl(fid);
        end
    end   
end

return

榫眼

谢谢主任！！！！！！！！！！！！！！

oyzy83scut

原帖由 w89986581 于 2009-3-11 13:20 发表

                               
登录/注册后可看大图

提供一个sysnoise输出结果转存程序(matlab)。计算完指向性以后，自动生成场点，重新计算场点结果（可以多个频率），输出result（格式为txt），然后调用下面这个函数。
仅供参考，自己多动手。

function [] = loa ...

请问这个程序是在sysnoise导出指向性数据后，在matlab读入之前导出的数据并运行的吗？

w89986581

是的。

榫眼

0.00000E+00  9.25155E-06 -4.56975E-06
  1.00000E+01  9.25140E-06 -4.56975E-06
  2.00000E+01  9.25097E-06 -4.56975E-06
  3.00000E+01  9.25030E-06 -4.56975E-06
  4.00000E+01  9.24949E-06 -4.56975E-06
  5.00000E+01  9.24862E-06 -4.56975E-06
  6.00000E+01  9.24781E-06 -4.56975E-06
  7.00000E+01  9.24715E-06 -4.56974E-06
  8.00000E+01  9.24671E-06 -4.56974E-06
  9.00000E+01  9.24656E-06 -4.56974E-06
  1.00000E+02  9.24671E-06 -4.56974E-06
  1.10000E+02  9.24715E-06 -4.56974E-06
  1.20000E+02  9.24781E-06 -4.56975E-06
  1.30000E+02  9.24862E-06 -4.56975E-06
  1.40000E+02  9.24949E-06 -4.56975E-06
  1.50000E+02  9.25030E-06 -4.56975E-06
  1.60000E+02  9.25097E-06 -4.56975E-06
  1.70000E+02  9.25140E-06 -4.56975E-06
  1.80000E+02  9.25155E-06 -4.56975E-06
  1.90000E+02  9.25140E-06 -4.56975E-06
  2.00000E+02  9.25097E-06 -4.56975E-06
  2.10000E+02  9.25030E-06 -4.56975E-06
  2.20000E+02  9.24949E-06 -4.56975E-06
  2.30000E+02  9.24862E-06 -4.56975E-06
  2.40000E+02  9.24781E-06 -4.56975E-06
  2.50000E+02  9.24715E-06 -4.56974E-06
  2.60000E+02  9.24671E-06 -4.56974E-06
  2.70000E+02  9.24656E-06 -4.56974E-06
  2.80000E+02  9.24671E-06 -4.56974E-06
  2.90000E+02  9.24715E-06 -4.56974E-06
  3.00000E+02  9.24781E-06 -4.56975E-06
  3.10000E+02  9.24862E-06 -4.56975E-06
  3.20000E+02  9.24949E-06 -4.56975E-06
  3.30000E+02  9.25030E-06 -4.56975E-06
  3.40000E+02  9.25097E-06 -4.56975E-06
  3.50000E+02  9.25140E-06 -4.56975E-06
  3.60000E+02  9.25155E-06 -4.56975E-06
    角度值              声压实部        声压虚部
这是图2的原始数据，是以直角坐标提出的，在MATLAB中用PLOT()可以
画出与SYSNOISE中选项XYPLOT一样的图，可在MATLAB用POLAR()确画不
出SYSNOISE中选项XYPOLAR的极坐标图如图2所示，这个问题已经让我困扰很久了
请大家指点一下

榫眼

polar ('上面数据第一列'，‘上面数据第二列’)，更本就不对

w89986581

别把弧度当成糊涂啦，呵呵。

榫眼

呵呵，急糊涂啦！谢谢主任！