- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我想计算最后一个的平均值,例如matlab中矩阵的每一列中的3个非零数字。这些列的末尾用零填充以创建相同长度的向量。
示例矩阵:
A = [5 6 3 5 6 8 9;
1 2 3 5 4 7 6;
0 1 2 3 4 5 6;
0 0 1 2 3 4 5;
0 0 0 1 2 3 4;
0 0 0 0 2 3 4;
0 0 0 0 2 3 4;
0 0 0 0 0 0 3]
最佳答案
可能有更有效的解决方案,但一种方法是使用 sum
来查找给定列中非零行的数量。然后通过使用 arrayfun
遍历所有列并对列中零之前的 N
行进行平均来获取 A 的平均值。
%// Number of elements to average
N = 3;
%// Last non-zero row in each column
lastrow = sum(A ~= 0, 1);
%// Ensure that we don't have any indices less than 1
startrow = max(lastrow - N + 1, 1);
%// Compute the mean for each column using the specified rows
means = arrayfun(@(k)mean(A(startrow(k):lastrow(k),k)), 1:size(A, 2));
示例
对于您的示例数据,这将产生:
3.0000 3.0000 2.0000 2.0000 2.0000 3.0000 3.6667
另一种方法是使用卷积来实际为您解决这个问题。您可以使用卷积核计算均值。如果您想要一个矩阵的所有 3 行组合的平均值,您的内核将是:
kernel = [1; 1; 1] ./ 3;
当与感兴趣的矩阵进行卷积时,这将计算输入矩阵中所有 3 行组合的平均值。
B = [1 2 3;
4 5 6;
7 8 9];
conv2(B, kernel)
0.3333 0.6667 1.0000
1.6667 2.3333 3.0000
4.0000 5.0000 6.0000
3.6667 4.3333 5.0000
2.3333 2.6667 3.0000
在下面的示例中,我执行此操作然后仅返回我们关心的区域的值(其中平均值仅由每列中的最后一个 N
非零值组成)
%// Find the last non-zero entry in each column
lastrow = sum(A ~= 0, 1);
%// Use convolution to compute the mean for every N rows
%// This will be applied to ALL of A
convmean = conv2(A, ones(N, 1)./N);
%// Select only the means that we care about
%// Because of the padding of CONV2, these will live at the rows
%// stored in LASTROW
means = convmean(sub2ind(size(convmean), lastrow, 1:size(A, 2)));
%// Now correct for cases where fewer than N samples were averaged
means = (means * N) ./ min(lastrow, N);
同样,输出是一样的
3.0000 3.0000 2.0000 2.0000 2.0000 3.0000 3.6667
我运行了一个快速测试脚本来比较这两种方法的性能。很明显,基于卷积的方法要快得多。
这是完整的测试脚本。
function benchmark()
dims = round(linspace(1, 1000, 100));
times1 = zeros(size(dims));
times2 = zeros(size(dims));
N = 3;
for k = 1:numel(dims)
A = triu(rand(dims(k)));
times1(k) = timeit(@()test_arrayfun(N, A));
A = triu(rand(dims(k)));
times2(k) = timeit(@()test_convolution(N, A));
end
figure;
plot(dims, times1);
hold on
plot(dims, times2);
legend({'arrayfun', 'convolution'})
xlabel('Dimension of A')
ylabel('Execution Time (seconds)')
end
function test_arrayfun(N, A)
%// Last non-zero row in each column
lastrow = sum(A ~= 0, 1);
%// Ensure that we don't have any indices less than 1
startrow = max(lastrow - N + 1, 1);
%// Compute the mean for each column using the specified rows
means = arrayfun(@(k)mean(A(startrow(k):lastrow(k),k)), 1:size(A, 2));
end
function test_convolution(N, A)
%// Find the last non-zero entry in each column
lastrow = sum(A ~= 0, 1);
%// Use convolution to compute the mean for every N rows
%// This will be applied to ALL of A
convmean = conv2(A, ones(N, 1)./N);
%// Select only the means that we care about
%// Because of the padding of CONV2, these will live at the rows
%// stored in LASTROW
means = convmean(sub2ind(size(convmean), lastrow, 1:size(A, 2)));
%// Now correct for cases where fewer than N samples were averaged
means = (means * N) ./ min(lastrow, N);
end
关于matlab - 如何确定最后一个的平均值matlab矩阵每列100个非零数,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/36288101/
在 Windows 世界中,什么是正确的名称。具有导出函数的老式 C++ DLL?不是 COM DLL,也不是 .NET DLL。我们以前通过调用 LoadLibrary() 和 GetProcAdd
目前我正在使用javaEE7,我有一个场景如下。在我的 JSF Web 应用程序中,我有一个事件监听器(不是 JSF 事件),当事件调用时,它会执行某些操作,然后将这些信息更新到我的 Web 应用程序
这不是 AJAX 请求/响应回调问题... 我正在使用 Dojo 1.5 构建网格。我正在尝试 dojo.connect具有功能的扩展/收缩按钮。我的问题是 grid.startup()在创建实际 D
非 Webkit Opera 是 very specific在某些功能中,因此通常通过 JavaScript 检测到 the following way . 但是,Opera Next 几乎是 Goo
我已查看以下链接中给出的所有日志,但未能找到 IP 地址: https://developer.couchbase.com/documentation/server/3.x/admin/Misc/Tr
我有一个命令行程序,它根据一组源文件生成一个我想在我的 Android gradle 构建 (A) 中使用的 jar 文件。这个命令行程序只是将一个 jar 文件存储在磁盘上的一个目录中。 我如何创建
下面的 htaccess 命令将所有非 www 转移到 http www RewriteEngine On RewriteCond %{HTTP_HOST} !^www\. RewriteRule ^
我正在使用自定义链接器脚本将内核镜像分为两部分。第一个是普通代码和数据,第二个是初始化代码和不再需要时将被丢弃的数据。初始化部分也不像内核本身那样在地址空间之间共享,因此如果 fork() 仍然存在(
这个问题在这里已经有了答案: Several unary operators in C and C++ (3 个答案) What is the "-->" operator in C++? (29
假设我有一个类设置如下: class A { public: virtual void foo() { printf("default implementation\n"); } }; c
#include using namespace std; int main(int argc, char *argv[]) { int i=-5; while(~(i)) {
近期,百度搜索引擎变化无常,很多企业站、行业站、门户站、论坛等站点遭到了降权,特别是比比贴分类信息网直接遭到了拔毛,这对于广大站长来说是一种打击,也是各个企业、行业的打击。 至今,很多网站已经恢复
我现在正在使用 IBM TPM v1332 + IBM TSS v1470 并尝试将一些基本关键字/密码存储到 TPM 上的非 volatile 内存中。我找到了两种方法。一种是创建一个密封对象并使用
我的 PHP 脚本中有一个正则表达式,如下所示: /(\b$term|$term\b)(?!([^)/iu 这与 $term 中包含的单词匹配,只要前后有单词边界并且它不在 HTML 标记内即可。 但
我想显示用户名称地址(请参阅 www.ipchicken.com ),但我唯一能找到的是 IP 地址。我尝试了反向查找,但也没有用: IPAddress ip = IPAddress.Parse(th
只有 UI 线程能够显示到屏幕上,还是其他线程也可以这样做? 最佳答案 不,您只能直接从 UI 线程访问 UI,但您可以编码来自其他线程的结果,例如使用 Control.Invoke 或 contro
我正在使用现代 Excel 滚动条(不是旧的 ActiveX 类型,即开发人员 > 插入 > 表单控件 > 滚动条)并且想检测它的值何时更改。我找不到有关此类对象的更改事件的任何信息。您可以在单击时分
当我使用这段代码时 IE 6 确实正确使用了指定的样式表,但所有其他浏览器在应该使用基本上声明的样式表时会忽略这两种样式表,如果您不是 IE,请使用此样式表。 有什么想法吗? 最佳答案 n
我想指定 2 mssql 表之间的关系。 付款类别和付款。 paymentcategory.id 加入 payout.category 列。 在 payout.json 模型中 我指定为外键:id,
我正在尝试制作非 volatile UDF,但似乎不可能。因此,这是我非常简单的test-UDF: Option Explicit Dim i As Integer Sub Main() i = 0
我是一名优秀的程序员,十分优秀!