关于arrayfun在复杂匿名函数中的用法请教

bainhome

今天测试arrayfun函数的用法，例如下面的这个函数：

给出其参变量a,b,c的变化规律如下列矩阵，第一列为a，第二列为b，第三列为c参变量的数据：

1. >> data=[2,10,pi/3;4.5,100,5*pi/6;6,1000,2*pi/3]
   
2. data =
   
3. 1.0e+003 *
   
4. 0.0020 0.0100 0.0010
   
5. 0.0045 0.1000 0.0026
   
6. 0.0060 1.0000 0.0021

复制代码

想通过定义符合匿名函数f：

1. >> f=@(a,b,c,x) a*sin(b*x+c)+2*a*sin(b/3*x-2*c)
   
2. f =
   
3. @(a,b,c,x)a*sin(b*x+c)+2*a*sin(b/3*x-2*c)

复制代码

一次在一幅figure中绘制出不同参变量下的f函数图形，如果代入data数据，当然就应当形成三条曲线，每次[a,b,c]对应取data(i,1),data(i,2),data(i,3)形成一组参数绘图。大致上看了看rocwoods大作中所有arrayfun函数的用法，发现老弟给出的命令中没有对多参数条件下的arrayfun绘制图像示例，不知道能否经由嵌套arrayfun来实现？另外所给事关plot和arrayfun结合的示例均是让plot的自变量与arrayfun自身的循环节奏同调，如下例(见rocwoods书page53,例4.2-1续)：

1. ...
   
2. aa=0:.01:2;
   
3. plot(aa,arrayfun(@(a) fzero(f(a),4),aa),'*-')
   
4. ...

复制代码

注：关于用循环实现的方式个人以为就不必讨论了，目的纯粹是研究arrayfun的用法。

qibbxxt

一个arrayfun的确比较难，我用了两个

1. clear;clc;close all
   
2. hold on
   
3. data=[2,10,pi/3;4.5,100,5*pi/6;6,1000,2*pi/3];
   
4. f=@(x,a,b,c) a*sin(b*x+c)+2*a*sin(b/3*x-2*c);
   
5. x=linspace(0,2);
   
6. arrayfun(@(a,b,c)plot(x,arrayfun(@(x)f(x,a,b,c),x),'color',rand(1,3)),data(:,1),data(:,2),data(:,3));
   
7. legend(strcat({'第'},num2str((1:size(data))','%d'),{'条曲线'}),'Location','best');

复制代码

rocwoods

用一个arrayfun也能实现，就是定义匿名函数f的时候，f应该接受向量化输入，这个问题比较简单，不用写成点乘的形式

1. f=@(x,a,b,c) a*sin(b*x+c)+2*a*sin(b/3*x-2*c);

复制代码

已经可以接受x为向量化形式了。所以稍微改动祁老弟的代码如下即可：

1. 
   
2. clear;clc;close all
   
3. hold on
   
4. data=[2,10,pi/3;4.5,100,5*pi/6;6,1000,2*pi/3];
   
5. f=@(x,a,b,c) a*sin(b*x+c)+2*a*sin(b/3*x-2*c);
   
6. x=linspace(0,2);
   
7. arrayfun(@(a,b,c)plot(x,f(x,a,b,c),'color',rand(1,3)),data(:,1),data(:,2),data(:,3));
   
8. legend(strcat({'第'},num2str((1:size(data))','%d'),{'条曲线'}),'Location','best');
   

复制代码

qibbxxt

3# rocwoods
恩，看来对函数的理解还应该再深刻一些
对问题的分析应该精益求精，这样才会有深入的进步，向吴兄学习

bainhome

刚回来看到。
昨天脑子进死循环了，单变量条件下看到示例是写成

1. arrayfun(@(a) f(a),...)

复制代码

想当然就按这个思路往下走，没想到arrayfun中需要把函数写成

1. arrayfun(@() f(x,a,b,c),...)

复制代码

ps:嵌套arrayfun的方式我也很喜欢。