Matlab解方程组结果提示Access to an object's fields is only permitted within its methods.
哪需要这么麻烦? 请看
for n=1:1:9 x = rand(1,10); plot(x); title(num2str(n)); pause;end
没有解析解，就计算数值解吧。
function myode45
[t,xyz]=ode45(@fun,0:1,[1 2 3]) % 'x(0)=1','y(0)=2','z(0)=3');
plot(t,xyz)
function f=fun(t,X)
x=X(1);y=X(2);z=X(3);
dx=10*(y-x);
dy=28*x-y-x*z;
dz=-2.6*x+x*y;
f=[dx;dy;dz];
结果：
t =
0
0.0044
0.0087
0.0131
0.0175
0.0357
0.0538
0.0720
0.0902
0.1120
0.1338
0.1555
0.1773
0.1938
0.2104
0.2269
0.2435
0.2561
0.2686
0.2812
0.2938
0.3064
0.3190
0.3316
0.3442
0.3553
0.3664
0.3775
0.3886
0.3981
0.4075
0.4170
0.4264
0.4388
0.4511
0.4635
0.4758
0.4902
0.5046
0.5191
0.5335
0.5484
0.5634
0.5784
0.5933
0.6078
0.6223
0.6368
0.6513
0.6650
0.6787
0.6923
0.7060
0.7198
0.7336
0.7473
0.7611
0.7748
0.7884
0.8021
0.8158
0.8294
0.8431
0.8567
0.8704
0.8841
0.8978
0.9115
0.9252
0.9391
0.9530
0.9670
0.9809
0.9857
0.9904
0.9952
1.0000
xyz =
1.0000 2.0000 3.0000
1.0449 2.1027 2.9975
1.0924 2.2100 2.9955
1.1426 2.3221 2.9938
1.1955 2.4394 2.9927
1.4479 2.9896 2.9958
1.7594 3.6561 3.0166
2.1420 4.4648 3.0678
2.6099 5.4465 3.1693
3.3059 6.8949 3.4017
4.1859 8.7089 3.8230
5.2910 10.9471 4.5609
6.6626 13.6388 5.8172
7.9038 15.9722 7.2827
9.3295 18.4930 9.3708
10.9296 21.0526 12.2541
12.6610 23.3885 16.0864
14.0190 24.7922 19.7005
15.3492 25.6225 23.8900
16.5748 25.6508 28.5360
17.6016 24.6885 33.4144
18.3310 22.6026 38.2107
18.6596 19.4541 42.5220
18.5140 15.4456 46.0040
17.8737 10.8723 48.4345
16.9118 6.6846 49.6332
15.6188 2.6685 50.0145
14.0581 -0.9609 49.7262
12.3054 -4.0693 48.9687
10.7265 -6.2572 48.1235
9.1119 -8.0333 47.2096
7.5014 -9.4210 46.3193
5.9272 -10.4599 45.5232
3.9640 -11.3679 44.6996
2.1496 -11.8411 44.1651
0.5135 -11.9534 43.9300
-0.9264 -11.7622 43.9718
-2.3470 -11.2093 44.3100
-3.4811 -10.3500 44.8766
-4.3287 -9.2345 45.5781
-4.8978 -7.9131 46.3257
-5.2092 -6.3830 47.0674
-5.2582 -4.7732 47.7115
-5.0759 -3.1644 48.2204
-4.6996 -1.6311 48.5855
-4.1893 -0.2793 48.8156
-3.5779 0.8948 48.9449
-2.9060 1.8607 48.9999
-2.2095 2.6059 49.0107
-1.5593 3.1069 49.0041
-0.9407 3.4197 48.9927
-0.3734 3.5589 48.9846
0.1289 3.5434 48.9831
0.5602 3.3932 48.9873
0.9091 3.1331 48.9940
1.1736 2.7892 48.9994
1.3554 2.3858 48.9995
1.4587 1.9504 48.9911
1.4922 1.5021 48.9734
1.4646 1.0605 48.9467
1.3853 0.6409 48.9124
1.2645 0.2557 48.8734
1.1129 -0.0848 48.8325
0.9407 -0.3737 48.7926
0.7572 -0.6077 48.7566
0.5701 -0.7865 48.7267
0.3877 -0.9105 48.7041
0.2164 -0.9828 48.6894
0.0609 -1.0080 48.6827
-0.0774 -0.9909 48.6832
-0.1931 -0.9378 48.6900
-0.2848 -0.8559 48.7018
-0.3523 -0.7523 48.7170
-0.3700 -0.7131 48.7227
-0.3851 -0.6724 48.7287
-0.3975 -0.6305 48.7348
-0.4074 -0.5878 48.7409
用matlab解微分方程组
解析解:[x,y]=dsolve('Dx=x-y-x*(x^2+y^2)','Dy=x+y-y*(x^2+y^2)','x(0)=2','y(0)=1')得到的结果是解析解没有找到。用数值解。在Matlab下输入:edit，然后将下面两行百分号之间的内容，复制进去，保存 function y=zhidao_rk4_5(t,x)x,y变量分别用x(1),x(2)表示 y...
如何用matlab求解微分方程组的解析解
如何用matlab求解微分方程组的解析解?为了说明问题。特列举一个例子。例：x"(t)+2x'(t)=x(t)+2y(t)-exp(-t)，y'(t)=4x(t)+3y(t)+4exp(-t)对于线性微分方程组，一般可以用dsolve()函数求解直接得到其解析解，其计算过程如下：syms t x(t) y(t) %声明变量 D2x=diff(x,2); %求...
matlab解三元一次方程组,解全为0怎么办?
你的问题不属于求解三元一次方程组，因为所求的变量是l、m、n，而给出的方程有四个，其中l²+m²+n²=1是约束条件方程，所以你的问题属于非线性规划问题。无法用solve函数来求解。该问题可以运用fmincon函数来求解。其求解方法如下：1、创建目标函数 2、创建约束函数 3、初定l、m...
matlab求矩阵方程组
这相当于线性方程组的第三行。所以相当于：\/\/ 代码如下 A = [A00 - 1,A10;A01,R*A2 - 1;e,*inv(I-R)*e];b = [0;0;1];x = A\\b;\/\/ Matlab自动求解线性方程组的语句，没记错的话，对于这种超定方程组，应该也能自动搞定，算出最小二乘解 \/\/ 代码结束 最后，求出来的x是个2...
matlab 如何解非线性方程组
使用solve函数。举个例子，解非线性方程组 x^2+y^3=10 x^3-y^2=1 其中x,y为方程组的未知量 在Matlab的命名窗口中输入：syms x y [x y]=solve('x^2+y=10','x^2-y^2=1','x','y')即可 输出计算结果为：x = (37^(1\/2)\/2 + 21\/2)^(1\/2)(21\/2 - 37^(1\/2)\/2)...
非线性方程组的解法matlab
用matlab求解非线性方程组方法，可以用下列方法来实现：方法一，使用solve函数求解 x = optimvar('x');y = optimvar('y');prob = optimproblem;prob.Objective = -x - y\/3;prob.Constraints.cons1 = x + y <= 2;prob.Constraints.cons2 = x + y\/4 <= 1;prob.Constraints.cons3 = x ...
matlab有变量的方程组求解,写出具体编程
matlab有变量的方程组求解方法。利用for循环语句，有vpasolve()函数解出x,y值。求解过程：a=0:0.1:(50*3.14\/180); x1=2720*cos(a)-280*sin(a); y1=2720*sin(a)+280*cos(a);for i=1:length(a)syms x y [x,y]=vpasolve((x-2312)^2+(y-100)^2-1450^2==0,(x-x1(i))^2...
MATLAB解带三角函数的方程组,如何编程
1、对于简单的三角方程，可以用solve()来求解，得到解析解。如 其求解代码和结果为 2、对于复杂的三角方程，可以用Broyden（拟牛顿法）等方法来求解，得到数值解。如 其求解过程如下 ①建立自定义函数文件，funm。m function y = funm(x)x=x(1,1)，y=x(2,1)y(1,1) = x(1,1) - 0.7...
matlab中solve解方程组
一。用matlab 中的solve函数 >>syms x y; %定义两个符号变量；>>[x ,y]=solve('y=2*x+3','y=3*x-7');%定义一个 2x1 的数组，存放x，y >>x >>x=10.0000 >>y >>y=23.0000 二。用matlab 中的反向斜线运算符（backward slash）分析：方程组可化为 2*x-y=-3；3*x-y...
怎么用MATLAB解这个六元方程组呀?急求求大佬解答
用solve（）可以得到，其求解过程如下：syms xA yA xB yB xC yC eq1=25-(xA^2+yA^2);eq2=25-((xB-5)^2+yB^2);eq3=9-(xC^2+(yC-6)^2);eq4=9-((xB-xC)^2+(yB-yC)^2);eq5=18-((xA-xC)^2+(yA-yC)^2);eq6=9-((xA-xB)^2+(yA-yB)^2);s=solve(eq1,eq2,...
都甄青可： 我的版本是2012a.(注:solve求解得到表达式) 结果如下: b = -log(cos((pi*F)/2))-log(cos((pi*F)/2))F = -a(1 - F*a)^(1/2)/(b^(1/2)*(b + 1)^(1/2))a(1 - F*a)^(1/2)/(b^(1/2)*(b + 1)^(1/2))f = atan((b^(1/2)*(a - 1))/(a(1 - F*a)^(1/2)*(b + 1)^(1/2)))...
新华区18453111778： 如何用matlab解非齐次线性方程组,其中方程的个数小于未知量的个数用matlab判断方程组否有解,如果有,求其通解:例如:a+b - 3c - d=13a - b - 3c+4d=4a+5... - ？
都甄青可：[答案] clear all A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; b=[1 4 0]' %输入矩阵A,b A;b; %输入矩阵A,b [m,n]=size(A); R=rank(A); B=[A b]; Rr=rank(B); format rat if R==Rr&R==n % n为未知数的个数,判断是否有唯一解 x=A\b; elseif R==Rr&R
新华区18453111778： matlab使用fsolve解一个非线性方程组,出现下面结果怎么办?我改大了MaxFunEvals的值,结果还是一样. - ？
都甄青可： ....上面都说了是它已经超出函数估值的极限了,你还把估值过程量再改大,试试改小吧.
新华区18453111778： Matrix is close to singular or badly scaled 出错提示求助 - ？
都甄青可： Matlab程序运行出现问题:Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 4.777592e-017. 分析:警告:矩阵接近奇异,即det(Matrix)接近0,解方程组得到的结果可能不正确.RCOND是条件数的倒数,RCOND越小,其越接近奇异!就是矩阵接近奇异啦,迭代结果可能不收敛到真实值了,建议换个初始迭代值试试.随着迭代的进行,矩阵中的数值太小,容易出现这种情况.
新华区18453111778： 遇到一个二元方程组,直接用matlab求解出现如下信息 - ？
都甄青可： 函数值接近0,与收敛容差相近,求解完成;梯度分析,问题解可能是奇异解;估计方程给的不好,稳定性很差,是自己随便写的方程吧...这个求解不可靠,需谨慎使用...
新华区18453111778： matlab出现Matrix is close to singular or badly scaled是怎么回事?谢谢 - ？
都甄青可： 矩阵奇异,解方程组得到的结果可能不正确
新华区18453111778： 数值计算方法 - ？
都甄青可： 占个位,明天下午再看看. %%%%%%%%%%%%%%%%%%%%%%%%%%%%% 一题:(你的题目中精度没有说清楚,应当是公式复制过,丢失信息了) 你改一下精度和初始值吧(自己设计的迭代法的收敛与初值关系比较大)f=inline('(x^2+2-...
新华区18453111778： 用matlab解方程组提示Explicit solution could not be found. - ？
都甄青可： 试试fsolve
新华区18453111778： matlab解方程中出现 is not a valid expression or equation.该怎么改 - ？
都甄青可： 方程表达式有误,中间的那个2cos(x),少了一个乘号,应写成2*cos(x);运算结果为x=0.30455140531179438887495863392766
新华区18453111778： matlab在求解微分方程时遇到的问题？
都甄青可： 我的也是2008a,遇到同样的问题.Root0f其实是根本无法解,只好用这种形式来表式一种结果,是Maple中的叫法.>> y=dsolve('(Dy)^2+y^2=1','y(0)=0')y =-1/2*RootOf(_Z^2-1)*(-1+exp(RootOf(_Z^2+1)*t)^2)*RootOf(_Z^2+1)/exp(RootOf(_Z^...
