matlab程序作分形蕨,matlab学习第一天笔记

>> phi=(1+sqrt(5))/2

phi =

1.6180

>> format long

>> phi

phi =

1.618033988749895

>> p=[1 -1 -1];   %表示多项式 x^2-x-1%

>> r=roots(p)

r =

-0.618033988749895

1.618033988749895

>> syms x

>> r=solve(1/x==x-1)     %直接求解方程不用转换多项式%

r =

1/2 - 5^(1/2)/2

5^(1/2)/2 + 1/2

>> pretty(r)

/ 1   sqrt(5) \

| - - ------- |

| 2      2    |

|             |

| sqrt(5)   1 |

| ------- + - |

\    2      2 /

>>  phi=r(1)   %获取第一个元素%

phi =

1/2 - 5^(1/2)/2

>> vpa(phi,50)  %变精度%

ans =

-0.61803398874989484820458683436563811772030917980576

>> phi=double(phi)

phi =

-0.618033988749895

>> f=@(x)1./x-(x-1)  %必须用./数组除算符%

f =

@(x)1./x-(x-1)

>> ezplot(f,[0,4])  %easyplot%

>>  phi=fzero(f,1)  %找零点%

phi =

1.618033988749895

>> phi=(1+sqrt(5))/2;

>> x=[0 phi phi 0 0];  %x包含5个元素%

>> y=[0 0 1 1 0];

>> u=[1 1];

>> v=[0 1];

>> plot(x,y,'b',u,v,'b--')   %x-y是蓝色实线，u-v是蓝色虚线%

>> text(phi/2,1.05,'\phi')   %标识字符，不同位置的标识字符%

>> text((1+phi)/2,-.05,'\phi-1')

>> text(-.05,.5,'1')

>> text(.5,-.05,'1')

>> axis equal  %使x，y方向坐标相等%

>> axis off   %使坐标消失%

>> set(gcf,'color','white')  %背景色设为白色%

以下保存在m文件中