PDA

查看完整版本 : [MATLAB基础] 【求助各位】新人不知如何调用.m文件


小媛爱蓝色
2015-06-03, 22:41
问题描述见附件,编写函数如下:

function F=myfun(x)
%%v=x(1);sigma(v)=x(2);
rou=0.95;
B=10885223;
E=1105358;
A=0.124118;
r=0.0271;
a1=0.3194;a2=-0.3566;a3=1.7815;a4=-1.8213;a5=1.3303;
F=[x(1)*(1-(1/sqrt(2*pi))*(exp(((log(x(1)/rou*B))+((x(2)^2)/2)/(x(2)))^2/(-2)))*(a1/(1+r*((log(x(1)/(rou*B))+((x(2))^2)/2)/(x(2))))+a2/(1+r*((log(x(1)/(rou*B))+((x(2))^2)/2)/(x(2)))^2)+a3/(1+r*((log(x(1)/(rou*B))+((x(2))^2)/2)/(x(2)))^3)+a4/(1+r*((log(x(1)/(rou*B))+((x(2))^2)/2)/(x(2)))^4)+a5/(1+r*((log(x(1)/(rou*B))+((x(2))^2)/2)/(x(2)))^5)))-rou*B*(1-(1/sqrt(2*pi))*(exp(((log(x(1)/rou*B))-((x(2))^2)/2)/(x(2)))^2)/(-2)))*(a1/(1+r*((log(x(1)/(rou*B))-((x(2))^2)/2)/(x(2))))+a2/(1+r*((log(x(1)/(rou*B))-((x(2))^2)/2)/(x(2)))^2)+a3/(1+r*((log(x(1)/(rou*B))-((x(2))^2)/2)/(x(2)))^3)+a4/(1+r*((log(x(1)/(rou*B))-((x(2))^2)/2)/(x(2)))^4)+a2/(1+r*((log(x(1)/(rou*B))-((x(2))^2)/2)/(x(2)))^5)))-E;
x(2)-A/(x(1)*(1-(1/sqrt(2*pi))*(exp(((log(x(1)/rou*B))+((x(2))^2)/2)/(x(2)))^2/(-2)))*(a1/(1+r*((log(x(1)/(rou*B))+((x(2))^2)/2)/(x(2))))+a2/(1+r*((log(x(1)/(rou*B))+((x(2))^2)/2)/(x(2)))^2)+a3/(1+r*((log(x(1)/(rou*B))+((x(2))^2)/2)/(x(2)))^3)+a4/(1+r*((log(x(1)/(rou*B))+((x(2))^2)/2)/(x(2)))^4)+a5/(1+r*((log(x(1)/(rou*B))+((x(2))^2)/2)/(x(2)))^5))))];
return

请问应如何调用,在指令窗应该输入什么?