matlab 简单问题请教

千亿星辰

请问这样编程画图可以不？
syms x w;
f=2./((2*pi).^0.5*x)*exp(-(log(x)-log(1.2)).^2*4);
f1=int(x*f,x,0,w);
f2=int(f,x,0,w);
f3=f1/f2;
fplot('f3',[0.1,10],1*exp-4)

mhkmars

1# 千亿星辰
你这个是接着昨天的问题么？感觉还是没有理清思路。我这里只是猜想你想通过几个表达式，积分后，f3随w变化的图线。

1. f1=@(w)quad(@(x)x.*2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4),0,w);
   
2. f2=@(w)quad(@(x)2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4),0,w);
   
3. w=linspace(0.1,10);
   
4. plot(w,arrayfun(@(w)f1(w)./f2(w),w))
   
5. grid;

复制代码

PS：对于这个代码，我在想能够简化不？就是通过两个匿名函数f1，f2，直接形成第三个f3的匿名函数，变量同为w？
请版上的各位朋友看看行不行，我搜索了一下，用匿名函数来生成匿名函数好像讨论的不多

messenger

匿名函数来生成匿名函数，rocwoods版主讨论的较多，你可以搜搜rocwoods版主以前的贴子或者参考rocwoods版主的《MATLAB高效编程技巧与应用：25个案例分析》那本书，比如：
http://forum.simwe.com/thread-885049-1-1.html
http://forum.simwe.com/thread-961191-1-1.html
http://forum.simwe.com/thread-931918-1-1.html

nwcwww

1# 千亿星辰  
你这个是接着昨天的问题么？感觉还是没有理清思路。我这里只是猜想你想通过几个表达式，积分后，f3随w变化的图线。
f1=@(w)quad(@(x)x.*2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4),0,w);
...
mhkmars 发表于 2011-3-31 16:36

用f1f2生成一个匿名函数f3倒是挺容易。。
不过如果直接ezplot(f3,[0.110]),
应该会比你用plot+arrayfun(@(w)f1(w)./f2(w),w)更慢。
而用plot+arrayfun(@(w)f3(w),w)的话和你的方法速度差不多。


用ezplot+f3:

1. 
   
2. 
   
3. clear;figure;tic;
   
4. f1=@(w)quad(@(x)x.*2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4),0,w);
   
5. f2=@(w)quad(@(x)2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4),0,w);
   
6. f3=@(w)f1(w)/f2(w);
   
7. ezplot(f3,[0.1,10]);grid;toc
   
8. 
   
9. 
   
10. Elapsed time is 0.786917 seconds.
    
11. 
    

复制代码

用plot+arrayfun+f3:

1. 
   
2. 
   
3. clear;figure;tic;
   
4. f1=@(w)quad(@(x)x.*2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4),0,w);
   
5. f2=@(w)quad(@(x)2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4),0,w);
   
6. f3=@(w)f1(w)/f2(w);
   
7. w=linspace(0.1,10);plot(w,arrayfun(@(w)f3(w),w));grid;toc
   
8. 
   
9. 
   
10. Elapsed time is 0.167830 seconds.
    
11. 
    

复制代码

用plot+arrayfun+f1+f2,也就是你的代码：

1. 
   
2. 
   
3. clear;figure;tic;
   
4. f1=@(w)quad(@(x)x.*2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4),0,w);
   
5. f2=@(w)quad(@(x)2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4),0,w);
   
6. w=linspace(0.1,10);plot(w,arrayfun(@(w)f1(w)./f2(w),w));grid;toc
   
7. 
   
8. 
   
9. Elapsed time is 0.180164 seconds.
   
10. 
    

复制代码

后两种方法在速度上没啥差别。。

mhkmars

3# messenger
谢谢斑竹，找到了，还是自己有点马虎。下面这样写就可以了，匿名函数生成匿名函数。

1. 
   
2. clc;clear;
   
3. f1=@(w)quad(@(x)x.*2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4),0,w);
   
4. f2=@(w)quad(@(x)2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4),0,w);
   
5. f3=@(w)f1(w)/f2(w)
   
6. fplot(f3,[0.1,10])
   

复制代码

mhkmars

4# nwcwww
学习了
觉得论坛上确实是藏龙卧虎咯，O(∩_∩)O哈哈~。又学习到了ezplot的用法了，还进行了优劣的比较。
我本人也才接触matlab不久，一个个帖子看下来，觉得有意思的就自己help一下，兴趣是好老师，慢慢的，沉浸下来，发现真的能够学到好多东西。
感谢所有帮助过我的朋友，谢谢，O(∩_∩)O哈哈~

qibbxxt

6# mhkmars
还可以这样再改下

1. 
   
2. clc;clear;
   
3. f=@(x)2./((2*pi).^0.5*x).*exp(-(log(x)-log(1.2)).^2*4);
   
4. f3=@(w)quad(@(x)x.*f(x),0,w)./quad(f,0,w);
   
5. fplot(f3,[0.1,10])
   

复制代码

mhkmars

7# qibbxxt
恩，这样写更简洁，学习一下。