Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#1 |
初级会员
注册日期: 2009-12-07
年龄: 38
帖子: 7
声望力: 0 ![]() |
![]()
本人是新手,想实现高斯分布的图,可是N=1时,应该会出现高斯分布,可是只出现最大值到最小值的一半的图形,望高手指教
clc,clear,close all; global N t0=0;tf=200; y1=0.66e-6;%设φm=0.66,φ(0,0)=φm*e-4 y2=0.0001;%设初始A=0 %options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]); [T,Y]=ode45('lie',[t0,tf],[y1,y2]); plot(T,Y(:,1),'r-') lie.m文件的内容如下 function dy=lie(t,y) N=2; dy=zeros(2,1); dy(1)=y(1)*N*exp(-y(2))-y(1); dy(2)=y(1); |
![]() |
![]() |