求解方程组

jiezouluxian

本人菜鸟一枚，有一个棘手的问题想请教各位

已知方程组(x-xE1)^2+(y-yE1)^2=0.69*0.69,yE1/xE1=y/x；

给出一组xE1和xE2，即xE1=[x1,x2,x3,x4,……,xn],yE1=[y1,y2,y3,y4,……,yn] 求每一个xE1，yE1值对应的x，y，即求x，y矩阵

不知哪位大侠能帮忙看看怎么弄

求指教，谢谢！！！急

悠悠为谁

1. 对应每一个参数方程组的解可能有两组解，所以得到的矩阵是2*m的。
2. 利用solve求解符号方程时候有时候需要指明未知量，要不然自动寻找的时候可能会出错。
改成这样表达：

1. >> S=solve('(x-xE1)^2+(y-yE1)^2=0.69*0.69','yE1/xE1=y/x','x,y');

复制代码

3. 采用你的参数，发现也能求出来：

1. >> x=subs(S.x);
   
2. >> y=subs(S.y);
   
3. >> x
   
4. 
   
5. x =
   
6. 
   
7.   Columns 1 through 11
   
8. 
   
9.    11.3450   10.6346    9.9322    9.2374    8.5496    7.8684    7.1930    6.5229    5.8573    5.1955    4.5370
   
10.    10.6550    9.9796    9.3125    8.6535    8.0022    7.3580    6.7202    6.0884    5.4620    4.8402    4.2226
    
11. 
    
12.   Columns 12 through 22
    
13. 
    
14.     3.8811    3.2274    2.5752    1.9243    1.2744    0.6251   -0.0237   -0.6720   -1.3199   -1.9673   -2.6140
    
15.     3.6086    2.9978    2.3896    1.7839    1.1802    0.5783   -0.0219   -0.6206   -1.2177   -1.8134   -2.4074
    
16. 
    
17.   Columns 23 through 33
    
18. 
    
19.    -3.2597   -3.9042   -4.5472   -5.1881   -5.8266   -6.4620   -7.0940   -7.7218   -8.3450   -8.9629   -9.5750
    
20.    -2.9998   -3.5904   -4.1790   -4.7654   -5.3494   -5.9306   -6.5088   -7.0837   -7.6550   -8.2223   -8.7852
    
21. 
    
22.   Columns 34 through 44
    
23. 
    
24.   -10.1806  -10.7792  -11.3702  -11.9532  -12.5276  -13.0930  -13.6490  -14.1952  -14.7312  -15.2569  -15.7718
    
25.    -9.3433   -9.8964  -10.4439  -10.9854  -11.5207  -12.0492  -12.5706  -13.0845  -13.5905  -14.0883  -14.5776
    
26. 
    
27.   Columns 45 through 55
    
28. 
    
29.   -16.2759  -16.7688  -17.2504  -17.7205  -18.1791  -18.6260  -19.0610  -19.4839  -19.8946  -20.2928  -20.6781
    
30.   -15.0580  -15.5293  -15.9912  -16.4435  -16.8859  -17.3181  -17.7401  -18.1514  -18.5519  -18.9412  -19.3189
    
31. 
    
32.   Columns 56 through 61
    
33. 
    
34.   -21.0503  -21.4086  -21.7527  -22.0815  -22.3943  -22.6900
    
35.   -19.6847  -20.0379  -20.3778  -20.7038  -21.0149  -21.3100
    
36. 
    
37. >> y
    
38. 
    
39. y =
    
40. 
    
41.   Columns 1 through 11
    
42. 
    
43.    19.6501   19.7189   19.7632   19.7846   19.7844   19.7639   19.7241   19.6661   19.5906   19.4984   19.3903
    
44.    18.4550   18.5042   18.5302   18.5342   18.5177   18.4818   18.4277   18.3563   18.2684   18.1650   18.0466
    
45. 
    
46.   Columns 12 through 22
    
47. 
    
48.    19.2667   19.1281   18.9751   18.8080   18.6271   18.4328   18.2254   18.0051   17.7722   17.5270   17.2697
    
49.    17.9138   17.7674   17.6076   17.4352   17.2504   17.0536   16.8454   16.6261   16.3960   16.1556   15.9052
    
50. 
    
51.   Columns 23 through 33
    
52. 
    
53.    17.0005   16.7197   16.4276   16.1245   15.8107   15.4863   15.1519   14.8077   14.4540   14.0911   13.7195
    
54.    15.6452   15.3759   15.0977   14.8108   14.5158   14.2128   13.9021   13.5841   13.2588   12.9267   12.5879
    
55. 
    
56.   Columns 34 through 44
    
57. 
    
58.    13.3395   12.9513   12.5553   12.1518   11.7410   11.3232   10.8986   10.4672   10.0292    9.5846    9.1332
    
59.    12.2424   11.8906   11.5324   11.1680   10.7973   10.4205   10.0375    9.6482    9.2526    8.8505    8.4417
    
60. 
    
61.   Columns 45 through 55
    
62. 
    
63.     8.6751    8.2099    7.7375    7.2573    6.7691    6.2722    5.7660    5.2498    4.7228    4.1841    3.6327
    
64.     8.0260    7.6031    7.1727    6.7343    6.2875    5.8318    5.3664    4.8907    4.4040    3.9054    3.3939
    
65. 
    
66.   Columns 56 through 61
    
67. 
    
68.     3.0677    2.4880    1.8925    1.2801    0.6496    0.0000
    
69.     2.8687    2.3287    1.7729    1.2002    0.6096    0.0000
    
70. 
    
71. >>

复制代码

悠悠为谁

试试这样：

1. >> S=solve('(x-a)^2+(y-b)^2=0.69^2','b/a=y/x');
   
2. >> a=1:10;
   
3. >> b=1:10;
   
4. >> x=subs(S.x);
   
5. >> y=subs(S.y);
   
6. >> x
   
7. 
   
8. x =
   
9. 
   
10.     1.4879    2.4879    3.4879    4.4879    5.4879    6.4879    7.4879    8.4879    9.4879   10.4879
    
11.     0.5121    1.5121    2.5121    3.5121    4.5121    5.5121    6.5121    7.5121    8.5121    9.5121
    
12. 
    
13. >> y
    
14. 
    
15. y =
    
16. 
    
17.     1.4879    2.4879    3.4879    4.4879    5.4879    6.4879    7.4879    8.4879    9.4879   10.4879
    
18.     0.5121    1.5121    2.5121    3.5121    4.5121    5.5121    6.5121    7.5121    8.5121    9.5121
    
19. 
    
20. >>

复制代码

jiezouluxian

悠悠为谁 发表于 2013-11-27 10:01
试试这样：

你好，首先谢谢您帮忙！我再问一下x，y应该是1行10的，怎么你给出的答案是2行十列的，不是每一组a，b对应一组x，y吗？

jiezouluxian

悠悠为谁 发表于 2013-11-27 10:01
试试这样：

并且我计算的时候怎么提示我的计算有错呢？
S=solve('(x-xE1)^2+(y-yE1)^2=0.69*0.69','yE1/xE1=y/x');
x=subs(S.x);y=subs(S.y);
??? Error using ==> mupadmex
Error in MuPAD command: dimensions do not match
[(Dom::Matrix(Dom::ExpressionField()))::_plus]

Error in ==> sym.subs>mupadsubs at 152
G = mupadmex('symobj::fullsubs',F.s,X2,Y2);

Error in ==> sym.subs at 127
    G = mupadsubs(F,X,Y);

jiezouluxian

悠悠为谁 发表于 2013-11-27 10:01
试试这样：

我的xE1=[11.0000000000000        10.3070850870268        9.62235497020440        8.94545204049313        8.27590739186780        7.61317008500258        6.95663451407019        6.30566540357414        5.65962010795945        5.01786802517944        4.37980705998986        3.74487717897531        3.11257118465870        2.48244289983295        1.85411299637148        1.22727272727273        0.601685829378287        -0.0228111397649641        -0.646309780853740        -1.26883191744072        -1.89033253981744        -2.51070344722481        -3.12977743290430        -3.74733288072536        -4.36309866426483        -4.97675925875500        -5.58795999299088        -6.19631238201803        -6.80139949225759        -7.40278129879194        -8.00000000000000        -8.59258525777001        -9.18005933230370        -9.76194207920853        -10.3377557732897        -10.9070297183152        -11.4693045951481        -12.0241364921359        -12.5711005516691        -13.1097941555761        -13.6398395598006        -14.1608858760329        -14.6726102852198        -15.1747183559470        -15.6669433305853        -16.1490442351134        -16.6208026662027        -17.0820181132606        -17.5325016856298        -17.9720681380295        -18.4005261224947        -18.8176666440566        -19.2232497611353        -19.6169896500581        -19.9985382450551        -20.3674677677774        -20.7232525695504        -21.0652508192881        -21.3926866730426        -21.7046336493994        -22];
yE1=[19.0525588832577        19.1115538107306        19.1466775864553        19.1593973185243        19.1510535385457        19.1228601989119        19.0759080690109        19.0111708921291        18.9295136557403        18.8317023498381        18.7184146350223        18.5902509074078        18.4477453241973        18.2913764355198        18.1215771493965        17.9387438331358        17.7432444230249        17.5354254730494        17.3156181217281        17.0841429940731        16.8413140838299        16.5874416805954        16.3228344184229        16.0478005284448        15.7626483791673        15.4676863855912        15.1632223632249        14.8495623962303        14.5270092810857        14.1958605988261        13.8564064605510        13.5089269628272        13.1536893821245        12.7909451307268        12.4209264908737        12.0438431393968        11.6598784720550        11.2691857353673        10.8718839742978        10.4680538069636        10.0577330429972        9.64091217067274        9.21752974979774        8.78746776305042        8.35054699817373        7.90652255737790        7.45507961835369        6.99582960306025        6.52830694509680        6.05196668265868        5.56618313988550        5.07024999227602        4.56338203860804        4.04471901880470        3.51333182047705        2.96823140251200        2.40838072854794        1.83270994385977        1.24013494479917        0.629579381032349        2.24333027298356e-15];

jiezouluxian

悠悠为谁 发表于 2013-11-26 20:58
1. 对应每一个参数方程组的解可能有两组解，所以得到的矩阵是2*m的。
2. 利用solve求解符号方程时候有时候 ...

麻烦再问一下，我之前按你的方法计算，为什么会报错呢？谢谢

mallocroon

路过，学习了。。。