声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3496|回复: 10

[Virtual.Lab] 辐射场指向性分析(Directivity)产生的问题

[复制链接]
发表于 2009-3-10 19:41 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
问题:如何用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 编辑 ]

图1

图1

图2

图2

图3

图3

图4

图4
回复
分享到:

使用道具 举报

发表于 2009-3-10 21:47 | 显示全部楼层

回复 楼主 榫眼 的帖子

呵呵,肯定是可以的,记得取幅度。
 楼主| 发表于 2009-3-11 11:10 | 显示全部楼层
谢谢主任的回答!
我想问的是幅度指的是幅值吗?实部与虚部的平方和开根号吗?
这样好像也不对
发表于 2009-3-11 13:20 | 显示全部楼层

回复 板凳 榫眼 的帖子

提供一个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
 楼主| 发表于 2009-3-11 21:10 | 显示全部楼层
谢谢主任!!!!!!!!!!!!!!
发表于 2009-3-13 11:45 | 显示全部楼层
原帖由 w89986581 于 2009-3-11 13:20 发表
提供一个sysnoise输出结果转存程序(matlab)。计算完指向性以后,自动生成场点,重新计算场点结果(可以多个频率),输出result(格式为txt),然后调用下面这个函数。
仅供参考,自己多动手。

function [] = loa ...


请问这个程序是在sysnoise导出指向性数据后,在matlab读入之前导出的数据并运行的吗?
发表于 2009-3-13 12:19 | 显示全部楼层

回复 6楼 oyzy83scut 的帖子

是的。
 楼主| 发表于 2009-3-14 09:36 | 显示全部楼层
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所示,这个问题已经让我困扰很久了
请大家指点一下
 楼主| 发表于 2009-3-16 16:26 | 显示全部楼层
polar ('上面数据第一列',‘上面数据第二列’),更本就不对
发表于 2009-3-16 16:35 | 显示全部楼层

回复 9楼 榫眼 的帖子

别把弧度当成糊涂啦,呵呵。
 楼主| 发表于 2009-3-16 19:22 | 显示全部楼层
呵呵,急糊涂啦!谢谢主任!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 19:31 , Processed in 0.083623 second(s), 27 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表