查看单个帖子
旧 2015-06-03, 22:41   #1
小媛爱蓝色
初级会员
 
注册日期: 2015-06-03
帖子: 1
声望力: 0
小媛爱蓝色 正向着好的方向发展
问题 【求助各位】新人不知如何调用.m文件

问题描述见附件,编写函数如下:
代码:
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
请问应如何调用,在指令窗应该输入什么?
上传的图像
文件类型: jpg dsg.jpg (17.3 KB, 4 次查看)
小媛爱蓝色 当前离线   回复时引用此帖