vpa编程

A. matlab自变量多个值,经公式运算，输出多个结果，如何编程

1、首先，建立一个自定义函数文件。如ex_fun.m，其文件内容

function [C,R]=ex_fun(X,Y,A)

syms C R

eq1=((C/((1/2)+X+(1/(0.013*R))))*(1-Y))-R;

eq2=((A*(11*(C-40)))/(8*(C+80)))-R;

s=solve(eq1,eq2);

C=vpa(s.C,5);

R=vpa(s.R,5);

end

2、再命令窗口或程序文件中，输入下列命令

X=5;Y=10;A=3;

[C,R]=ex_fun(X,Y,A)

3、计算结果

C = 40.， -6.929， -84.139

R = 0.， -2.649， 123.70

B. 用matlab做，牛顿迭代法

function [ A ] = cal( a,b,v )%a,b表示区间，v是精度

i=1;

x = (a+b)/2;

A=[i x];

t = x-(x^3-x-1)/(3*x^2-1);%迭代函数

while(abs(t-x)>v)

i=i+1;

x = t;

A = [A;i x];

t = x-(x^3-x-1)/(3*x^2-1);%迭代函数

end

A = [A;i+1 t];

end

运行结果：

>> format long;

>> cal(1,2,0.00001)

ans =

1.000000000000000 1.500000000000000

2.000000000000000 1.347826086956522

3.000000000000000 1.325200398950907

4.000000000000000 1.324718173999054

5.000000000000000 1.324717957244790

C. matlab编程求解方程组

如何求解多组二元一次方程组呢？题主的想法是对的，是要用循环来求解，其求解过程：

1、首先将t，h数据赋值给t1，h1，即

t1=tan(t);

h1=h.^2;

2、使用for循环语句，求解t(i)，h(i)对应的a，b值，即

for i=1:6

i

t=t1(i);h=h1(i);

syms a b

eqn1 = a > 0;

eqn2 = b/a==t;

eqn3 = a^2+b^2==h;

eqns = [eqn1 eqn2 eqn3];

S = solve(eqns,[a b]);

a=vpa(S.a)

b=vpa(S.b)

end

3、运行上述代码，可以得到如下结果。