MATLAB 板


LINE

最後一篇, 这篇会提出一个简单的 MEX file 的 Framework 废话不多说, 直接看 example, 程式码才是最好的说明书 这个 example 实际上是 vol.2 和 vol.3 的合体 如果前面的程式都看懂了, 这个程式应该不会有任何困难 建议把 code 复制到 Matlab 的 Editor 来看, 比较不会伤眼睛 :p 范例程式接受一个 2-D 的 input, 把每个元素 +1 以後输出 output = input + 1; 在程式码里, 我把程式区分了很多块 以一个标准的 MEX 档来说, 大概会有这几个部分 * 取得输入参数的资料 (eg. dimension) * 配置输出参数 (Output Allocation) * 根据输入, 计算输出 (Data Processing) 把下面的程式看懂, 碰到要写 MEX 的时候 直接套下面这个范例, 应该可以省不少时间 :) #include "mex.h" #include <math.h> #include <stdio.h> // Program framework for input, output and processing // usage: // input = zeros(2, 10) // input(1, :) = 1:10 // input(2, :) = 10:-1:1 // mex test3.c // a = test3(input) // note: type all commands above in Matlab Command Window void mexFunction( int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { int i, j, k; // input int input_dim_x; int input_dim_y; double *in; // pointer to process content of the input // output int output_dimension[2]; // note: you should specify output dimension int output_dim_x; int output_dim_y; double *out; /* -------------------------------- */ /* NECESSARY (input processing) */ in = mxGetPr(prhs[0]); // get data pointer // mx: Matrix // P: Pointer // r: real // mxGetPr() has a counterprt mxGetPi(); input_dim_x = mxGetM(prhs[0]); input_dim_y = mxGetN(prhs[0]); // print some message about input data printf("nrhs: %d\n", nrhs); printf("mxGetM(prhs[0]): %d\n", mxGetM(prhs[0])); printf("mxGetN(prhs[0]): %d\n", mxGetN(prhs[0])); // print input content for(i=0; i<input_dim_x; i++) // x for(j=0; j<input_dim_y; j++) // y // notice: data type is "float", you shall use "%f" insted of "%d" printf("%f\n", in[i + j*input_dim_x]); /* -------------------------------- */ /* NECESSARY (output allocation) */ // specify output matrix's dimension output_dim_x = input_dim_x; output_dim_y = input_dim_y; output_dimension[0] = output_dim_x; // mxCreateNumericArray() required output_dimension[1] = output_dim_y; // Allocate the output matrix plhs[0] = mxCreateNumericArray( 2, output_dimension, mxDOUBLE_CLASS, mxREAL); // Get output matrix's data pointer out = mxGetPr(plhs[0]); // just point at the start /* -------------------------------- */ /* NECESSARY (data processing) */ // Copy and add 1 for(i=0; i<output_dim_x; i++) for(j=0; j<output_dim_y; j++) { int temp = i + j*input_dim_x; out[temp] = in[temp]*2; } } --



※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.128.237
1F:推 kaishi:推荐这篇文章 我也在 140.113.128.XXX 哈 05/17 11:24
2F:推 Evanny:范例似乎是将输入乘上两倍 而不是加一喔 考我有没有认真看 05/17 23:19
3F:→ Evanny:真巧妙 刚好过了一年 05/17 23:20







like.gif 您可能会有兴趣的文章
icon.png[问题/行为] 猫晚上进房间会不会有憋尿问题
icon.pngRe: [闲聊] 选了错误的女孩成为魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一张
icon.png[心得] EMS高领长版毛衣.墨小楼MC1002
icon.png[分享] 丹龙隔热纸GE55+33+22
icon.png[问题] 清洗洗衣机
icon.png[寻物] 窗台下的空间
icon.png[闲聊] 双极の女神1 木魔爵
icon.png[售车] 新竹 1997 march 1297cc 白色 四门
icon.png[讨论] 能从照片感受到摄影者心情吗
icon.png[狂贺] 贺贺贺贺 贺!岛村卯月!总选举NO.1
icon.png[难过] 羡慕白皮肤的女生
icon.png阅读文章
icon.png[黑特]
icon.png[问题] SBK S1安装於安全帽位置
icon.png[分享] 旧woo100绝版开箱!!
icon.pngRe: [无言] 关於小包卫生纸
icon.png[开箱] E5-2683V3 RX480Strix 快睿C1 简单测试
icon.png[心得] 苍の海贼龙 地狱 执行者16PT
icon.png[售车] 1999年Virage iO 1.8EXi
icon.png[心得] 挑战33 LV10 狮子座pt solo
icon.png[闲聊] 手把手教你不被桶之新手主购教学
icon.png[分享] Civic Type R 量产版官方照无预警流出
icon.png[售车] Golf 4 2.0 银色 自排
icon.png[出售] Graco提篮汽座(有底座)2000元诚可议
icon.png[问题] 请问补牙材质掉了还能再补吗?(台中半年内
icon.png[问题] 44th 单曲 生写竟然都给重复的啊啊!
icon.png[心得] 华南红卡/icash 核卡
icon.png[问题] 拔牙矫正这样正常吗
icon.png[赠送] 老莫高业 初业 102年版
icon.png[情报] 三大行动支付 本季掀战火
icon.png[宝宝] 博客来Amos水蜡笔5/1特价五折
icon.pngRe: [心得] 新鲜人一些面试分享
icon.png[心得] 苍の海贼龙 地狱 麒麟25PT
icon.pngRe: [闲聊] (君の名は。雷慎入) 君名二创漫画翻译
icon.pngRe: [闲聊] OGN中场影片:失踪人口局 (英文字幕)
icon.png[问题] 台湾大哥大4G讯号差
icon.png[出售] [全国]全新千寻侘草LED灯, 水草

请输入看板名称,例如:Soft_Job站内搜寻

TOP