杭电Matlab与仿真-MATLAB基本语法（八）-匿名函数

实验目的

1. 掌握函数的表示与计算；
2. 掌握函数绘制的方法；
3. 认识函数的极点以及零点分析的方法。

实验内容与要求

1. 熟悉内联函数，掌握子函数、匿名函数的表示方法；
2. 掌握fplot()、ezplot()函数的运用方法并能够进行函数图像的绘制；
3. 掌握fminbnd()、fminsearch()、 fzero()等函数的运用，能够对函数的零点和极点进行分析；
4. 能够运用匿名函数解决各种实际问题。

实验程序与结果

1. 对于a=[0,0.01,0.02,…2]，求下列方程响应的x值，并画出a和相应x的图像
   

clear;
clc;
s=0:0.01:2;
f=@(a)@(x)exp(x)+x^a+x^sqrt(x)-100; 
m=arrayfun(@(a)fzero(f(a),3),s);
plot(s,m)

2. 画出下列隐函数中z关于x,y的图形。
   

clear;
clc;
m=zeros(41);
z=@(x,y)fzero(@(z)sin((z*x-0.5).^2+2.*x.*y.^2-z/10).*exp(-((x-0.5-exp(-y+z)).^2+y.^2-z/5+3))-z,1)
for i=-10:0.5:10
    for j=-10:0.5:10
        m(i*2+21,j*2+21)=z(i,j);
    end
end
x=-10:0.5:10;
y=-10:0.5:10;
[X,Y]=meshgrid(x,y);      
mesh(X,Y,m)

3. 生物韵律。生物韵律是100多年前被发现的，1960年左右进入我们的公众文化。现在你仍然可以发现很多能给出个性化生物韵律或兜售能计算生物韵律软件的网站。生物韵律是建立在三个正弦周期就能影响我们生活这样的观念上的。体力周期为23天，感情周期为28天，智力周期为33天。对每个人来说，这些周期的起点是出生的那一天。假设小民出生于1987年2月16日，请绘制出小民2016年4月25日前后8周期间的生物韵曲线，并绘制出自己在以上时间内的生物韵律曲线。假设体力、感情、智力峰值为100。（提示：需要用到datenum函数）.

function swgl(d1,d2)
dt=rem((d2-d1),23);
dq=rem((d2-d1),28); 
dz=rem((d2-d1),33);
x=-56:1:56;
yt=100*sin(2*pi/23*(x+dt));
yq=100*sin(2*pi/28*(x+dq));
yz=100*sin(2*pi/33*(x+dz));
plot(x,yt,'r');
hold on;
plot(x,yq,'g');
hold on;
plot(x,yz,'b');
legend('体力','智力','情感')
hold off;

clear;
clc;
d1=datenum(1987,2,16);
d2=datenum(2016,4,25);
swgl(d1,d2)