搜搜考试网matlab绘制隐函数图象
来源：学生作业帮 编辑：搜搜考试网作业帮 分类：综合作业 时间：2024/05/02 08:43:38
matlab绘制隐函数图象
char a b f;
syms x y
f=exp((1/(5*sqrt(2*pi)))*exp(-(y-100)^2/50))；
ezplot('int(int(f,y,0,969600*b-(20*a+20)*(2000-x)),x,500,2000)=500')
char a b f;
syms x y
f=exp((1/(5*sqrt(2*pi)))*exp(-(y-100)^2/50))；
ezplot('int(int(f,y,0,969600*b-(20*a+20)*(2000-x)),x,500,2000)=500')
F = @(x,y) exp((1/(5*sqrt(2*pi)))*exp(-(y-100).^2/50));
z = @(a,b) quadl(@(x) arrayfun(@(x) quadl(@(y) F(x,y),0,969600*b-(20*a+20).*(2000-x)),x),500,2000);
[A,B] = meshgrid(0.5:.2:2, 0:.05:0.5);
Z = arrayfun(@(a,b) z(a,b),A,B);
surf(A,B,Z)
xlabel( 'a' );
ylabel( 'b' );
再问： 运行不出来啊！
再答： 你用的Matlab是什么版本？
再问： 加你QQ行吗 想仔细了解一下 拜谢 我的是295118639
再答： 7.1版以上可用
再问： o 我写的那个可以吗 错在哪啊
再答： 1 你的问题是曲面，不是曲线。 2 该积分没有符号解，只能求数值解。 F = @(x,y) exp((1/(5*sqrt(2*pi)))*exp(-(y-100).^2/50)); z = @(a,b) quadgk(@(x) arrayfun(@(x) quadgk(@(y) F(x,y),0,969600*b-(20*a+20).*(2000-x)),x),500,2000) - 500; [A,B] = meshgrid(0.5:.2:2, 0:.05:0.5); Z = arrayfun(@(a,b) z(a,b),A,B); surf(A,B,Z) xlabel( 'a' ); ylabel( 'b' );
再问： 这是第二题的答案？
再答： 上面的代码是画出积分函数的曲面图。 而求出积分值为500时，a与b的关系，则可用下面的代码。 对于计算中出现的Warning，不必深究。 F = @(x,y) exp((1/(5*sqrt(2*pi)))*exp(-(y-100).^2/50)); ymax = @(x,a,b) 969600*b-(20*a+20).*(2000-x); Q = @(a,b) quad2d(F,500,2000,0,@(x) ymax(x,a,b)) - 500 A = 0:0.1:1; n = length(A); B = zeros(1,n); b0 = 0.01; % Make a starting guess at the solution for i = 1:length(A) a = A(i); [b,fval] = fsolve(@(b) Q(a,b),b0) B(i) = b; % b0 = b; end figure;plot(A,B,'k'); xlabel( 'a' ); ylabel( 'b' ); .......................... 对于b μ σ之间的三维图，你可参考上面的代码，多加一重循环即可。 .......... 请到我的百度空间， 标题为 MATLAB积分方程的求解
再问： 高人，我的电脑装不了7.0以上版本，你能帮忙吧这加了条件的 代码组织一下给我吗 非常感谢 题目在追问里
z = @(a,b) quadl(@(x) arrayfun(@(x) quadl(@(y) F(x,y),0,969600*b-(20*a+20).*(2000-x)),x),500,2000);
[A,B] = meshgrid(0.5:.2:2, 0:.05:0.5);
Z = arrayfun(@(a,b) z(a,b),A,B);
surf(A,B,Z)
xlabel( 'a' );
ylabel( 'b' );
再问： 运行不出来啊！
再答： 你用的Matlab是什么版本？
再问： 加你QQ行吗 想仔细了解一下 拜谢 我的是295118639
再答： 7.1版以上可用
再问： o 我写的那个可以吗 错在哪啊
再答： 1 你的问题是曲面，不是曲线。 2 该积分没有符号解，只能求数值解。 F = @(x,y) exp((1/(5*sqrt(2*pi)))*exp(-(y-100).^2/50)); z = @(a,b) quadgk(@(x) arrayfun(@(x) quadgk(@(y) F(x,y),0,969600*b-(20*a+20).*(2000-x)),x),500,2000) - 500; [A,B] = meshgrid(0.5:.2:2, 0:.05:0.5); Z = arrayfun(@(a,b) z(a,b),A,B); surf(A,B,Z) xlabel( 'a' ); ylabel( 'b' );
再问： 这是第二题的答案？
再答： 上面的代码是画出积分函数的曲面图。 而求出积分值为500时，a与b的关系，则可用下面的代码。 对于计算中出现的Warning，不必深究。 F = @(x,y) exp((1/(5*sqrt(2*pi)))*exp(-(y-100).^2/50)); ymax = @(x,a,b) 969600*b-(20*a+20).*(2000-x); Q = @(a,b) quad2d(F,500,2000,0,@(x) ymax(x,a,b)) - 500 A = 0:0.1:1; n = length(A); B = zeros(1,n); b0 = 0.01; % Make a starting guess at the solution for i = 1:length(A) a = A(i); [b,fval] = fsolve(@(b) Q(a,b),b0) B(i) = b; % b0 = b; end figure;plot(A,B,'k'); xlabel( 'a' ); ylabel( 'b' ); .......................... 对于b μ σ之间的三维图，你可参考上面的代码，多加一重循环即可。 .......... 请到我的百度空间， 标题为 MATLAB积分方程的求解
再问： 高人，我的电脑装不了7.0以上版本，你能帮忙吧这加了条件的 代码组织一下给我吗 非常感谢 题目在追问里