举个简单的例子:
a = [1 2i];
x = zeros(1,length(a));
for n=1:length(a)
x(n) = isreal(a(n));
end
在尝试对代码进行矢量化时,我尝试了:
y = arrayfun(@isreal,a);
但结果并不相同:
x =
1 0
y =
0 0
我做错了什么?
答案 0 :(得分:8)
这肯定是一个错误,但这是一个解决方法:
>> y = arrayfun(@(x) isreal(x(1)),a)
ans =
1 0
为什么这样做?我不是完全肯定,但是当你在调用ISREAL之前对变量执行索引操作时,它会从中移除“复杂”属性。数组元素,如果虚数组件为零。在命令窗口中尝试:
>> a = [1 2i]; %# A complex array
>> b = a(1); %# Indexing element 1 removes the complex attribute...
>> c = complex(a(1)); %# ...but we can put that attribute back
>> whos
Name Size Bytes Class Attributes
a 1x2 32 double complex
b 1x1 8 double %# Not complex
c 1x1 16 double complex %# Still complex
显然,ARRAYFUN必须在内部维护它传递给ISREAL的数组元素的“复杂”属性,因此即使虚数组件为零,也将它们视为复数。
答案 1 :(得分:8)
知道MATLAB分别存储矩阵的实部/复杂部分可能会有所帮助。请尝试以下方法:
>> format debug
>> a = [1 2i];
>> disp(a)
Structure address = 17bbc5b0
m = 1
n = 2
pr = 1c6f18a0
pi = 1c6f0420
1.0000 0 + 2.0000i
其中pr
是指向包含所有值的实部的内存块的指针,pi
指向矩阵中所有值的复杂部分的指针。由于所有元素都存储在一起,因此在这种情况下它们都有一个复杂的部分。
现在比较这两种方法:
>> arrayfun(@(x)disp(x),a)
Structure address = 17bbcff8
m = 1
n = 1
pr = 1bb8a8d0
pi = 1bb874d0
1
Structure address = 17c19aa8
m = 1
n = 1
pr = 1c17b5d0
pi = 1c176470
0 + 2.0000i
与
>> for n=1:2, disp(a(n)), end
Structure address = 17bbc930
m = 1
n = 1
pr = 1bb874d0
pi = 0
1
Structure address = 17bbd180
m = 1
n = 1
pr = 1bb874d0
pi = 1bb88310
0 + 2.0000i
因此,当您在for循环中访问a(1)
时,返回的值(在ans
变量中)具有零复数部分(null pi
),因此被认为是真实的。
另一方面,ARRAYFUN似乎直接访问矩阵的值(不在ANS变量中返回它们),因此它可以访问非{null}的pr
和pi
指针因此,所有元素都被认为是非实际的。
请记住这只是我的解释,我可能会弄错......
答案 2 :(得分:3)
对此问题的回答非常晚...... MATLAB函数ISREAL以一种非常反直觉的方式运行,用于多种用途。它告诉你一个给定的数组作为一个整体没有复杂的部分 - 它告诉你关于存储,它并没有真正告诉你任何关于<数组中的em> values 。在这方面,它有点像ISSPARSE功能。所以,例如
isreal(complex(1)) % returns FALSE
在MATLAB中你会发现某些操作会自动修剪任何全零虚部。所以,例如
x = complex(1);
isreal(x); % FALSE, we just forced there to be an imaginary part
isreal(x(1)); % TRUE - indexing realised it could drop the zero imaginary part
isreal(x(:)); % FALSE - "(:)" indexing is just a reshape, not real indexing
简而言之,MATLAB确实需要一个函数来回答问题“这个值是否具有零虚部”,在元素上以元素方式表示。