MATLAB里函数命令与用法实例.docVIP

MATLAB各类函数命令和用法实例（全） 1.符号函数 Ft=@(t)0*t;%初始化 ft=@(t)0*t; Ft=@(t)1-exp(-(t/theta).^beta);%括号内为自变量 ft=@(t)(beta/theta)*((t/theta).^(beta-1)).*(exp(-(t/theta).^beta)); %符号函数的求解过程 F=@(beta)1/beta+sum(sigma111.*log(fault111))/r-sum((fault111.^(beta)).*log(fault111))/(sum((fault111.^(beta)))); beta=fsolve(F,1); 2.递归运算示例 建立control555.m的M文件 function y=control555(x) if x==1 y=1; else y=x*control555(x-1); end 然后在命令窗口中输入 control555(5) ans =120 3.fslove数值法解方程组，程序示例 myfun.m文件程序如下： function F = myfun(x) t00=[1350 1 502.4 1 836.5 1 596.8 1 1323.4 1 729.3 1 295.4 0 676.2 1 939.6 1 306.58 0 12.68 0 394.6 0 202.1 0 88.2 1 344.7 1 296.6 1 56.4 1 128.2 0 42.8 1 163.3 1 227.7 1 18.6 0 31.1 1 144.5 1 55.6 1 102.8 1 87 1 46.94 1 ]; for i=1:length(t00) t000(i)=sum(t00(1:i,1)); end % t=t0(:,1); t=t000; c=t00(:,2); tt=[0;t]; tt(end)=[]; %下面求M，N for i=1:length(t) m(i)=sum(c(1:i)); n(i)=i-m(i); end m=m;n=n; t4=t(2:length(t)); t44=tt(2:length(tt)); m4=m(2:length(t)); n4=n(2:length(t)); %方程组 x=(a,p,sigma,b) F = [length(c)/x(1)-(1/(x(4)+1))*sum((x(2).^m).*(x(3).^n).*(t.^(x(4)+1))-(x(2).^m).*(x(3).^n).*(tt.^(x(4)+1))); (1/x(2))*sum(c.*m)-(x(1)/(x(4)+1))*sum(m.*(x(2).^m).*(x(3).^n).*(t.^(x(4)+1))-m.*(x(2).^m).*(x(3).^n).*(tt.^(x(4)+1))); (1/x(3))*sum(c.*n)-(x(1)/(x(4)+1))*sum(n.*(x(2).^m).*(x(3).^n).*(t.^(x(4)+1))-n.*(x(2).^m).*(x(3).^n).*(tt.^(x(4)+1))); (1/x(4))*sum(c.*log(t))-(x(1)/((x(4)+1).^2))*(sum((x(2).^m4).*(x(3).^n4).*(t4.^(x(4)+1)).*((x(4)+1).*log(t4)-1)-(x(2).^m4).*(x(3).^n4).*(t44.^(x(4)+1)).*((x(4)+1).*log(t44)-1))+sum((x(2).^m(1)).*(x(3).^n(1)).*(t(1).^(x(4)+1)).*((x(4)+1).*log(t(1))-1))) ]; main.m程序如下： clear,clc x0 = [1 1 1 1]; % Make a starting guess at the solution options=optimset(Display,iter); % Option to display output [x,fval] = fsolve(@myfun,x0,options) % Call optimizer 4.符号函数求定积分 syms x int(x^7,0,1) %定积分 int(1/x,1,2)%定积分 5.求极限 syms x; limit((x-2)/(x^2-4),2)% 求函数在某点极限 limit(1/x,x,0,left); % 求函数在x趋于0的左极限 6.符号函数求根、作图 syms x; roots=solve(x^2