im trying to solve this system of 10 equations

3 ビュー (過去 30 日間)
joshua payne
joshua payne 2023 年 2 月 6 日
回答済み: Walter Roberson 2023 年 2 月 6 日
clc
clear all
M=10
M = 10
deltax=2*pi/M
deltax = 0.6283
syms f1 f2 f3 f4 f5 f6 f7 f8 f9 f10
U1=f1*exp(1i*deltax)+f2*exp(1i*2*deltax)+f3*exp(1i*3*deltax)+f4*exp(1i*4*deltax)+f5*exp(1i*5*deltax)+f6*exp(1i*6*deltax)+f7*exp(1i*7*deltax)+f8*exp(1i*8*deltax)+f9*exp(1i*9*deltax)+f10*exp(1i*10*deltax)==cos(deltax)
U1 = 
U2=f1*exp(2*1i*deltax)+f2*exp(2*1i*2*deltax)+f3*exp(2*1i*3*deltax)+f4*exp(2*1i*4*deltax)+f5*exp(2*1i*5*deltax)+f6*exp(2*1i*6*deltax)+f7*exp(2*1i*7*deltax)+f8*exp(2*1i*8*deltax)+f9*exp(2*1i*9*deltax)+f10*exp(2*1i*10*deltax)==cos(2*deltax)
U2 = 
U3=f1*exp(3*1i*deltax)+f2*exp(3*1i*2*deltax)+f3*exp(3*1i*3*deltax)+f4*exp(3*1i*4*deltax)+f5*exp(3*1i*5*deltax)+f6*exp(3*1i*6*deltax)+f7*exp(3*1i*7*deltax)+f8*exp(3*1i*8*deltax)+f9*exp(3*1i*9*deltax)+f10*exp(3*1i*10*deltax)==cos(3*deltax)
U3 = 
U4=f1*exp(4*1i*deltax)+f2*exp(4*1i*2*deltax)+f3*exp(4*1i*3*deltax)+f4*exp(4*1i*4*deltax)+f5*exp(4*1i*5*deltax)+f6*exp(4*1i*6*deltax)+f7*exp(4*1i*7*deltax)+f8*exp(4*1i*8*deltax)+f9*exp(4*1i*9*deltax)+f10*exp(4*1i*10*deltax)==cos(4*deltax)
U4 = 
U5=f1*exp(5*1i*deltax)+f2*exp(5*1i*2*deltax)+f3*exp(5*1i*3*deltax)+f4*exp(5*1i*4*deltax)+f5*exp(5*1i*5*deltax)+f6*exp(5*1i*6*deltax)+f7*exp(5*1i*7*deltax)+f8*exp(5*1i*8*deltax)+f9*exp(5*1i*9*deltax)+f10*exp(5*1i*10*deltax)==cos(5*deltax)
U5 = 
U6=f1*exp(6*1i*deltax)+f2*exp(6*1i*2*deltax)+f3*exp(6*1i*3*deltax)+f4*exp(6*1i*4*deltax)+f5*exp(6*1i*5*deltax)+f6*exp(6*1i*6*deltax)+f7*exp(6*1i*7*deltax)+f8*exp(6*1i*8*deltax)+f9*exp(6*1i*9*deltax)+f10*exp(6*1i*10*deltax)==cos(6*deltax)
U6 = 
U7=f1*exp(7*1i*deltax)+f2*exp(7*1i*2*deltax)+f3*exp(7*1i*3*deltax)+f4*exp(7*1i*4*deltax)+f5*exp(7*1i*5*deltax)+f6*exp(7*1i*6*deltax)+f7*exp(7*1i*7*deltax)+f8*exp(7*1i*8*deltax)+f9*exp(7*1i*9*deltax)+f10*exp(7*1i*10*deltax)==cos(7*deltax)
U7 = 
U8=f1*exp(8*1i*deltax)+f2*exp(8*1i*2*deltax)+f3*exp(8*1i*3*deltax)+f4*exp(8*1i*4*deltax)+f5*exp(8*1i*5*deltax)+f6*exp(8*1i*6*deltax)+f7*exp(8*1i*7*deltax)+f8*exp(8*1i*8*deltax)+f9*exp(8*1i*9*deltax)+f10*exp(8*1i*10*deltax)==cos(8*deltax)
U8 = 
U9=f1*exp(9*1i*deltax)+f2*exp(9*1i*2*deltax)+f3*exp(9*1i*3*deltax)+f4*exp(9*1i*4*deltax)+f5*exp(9*1i*5*deltax)+f6*exp(9*1i*6*deltax)+f7*exp(9*1i*7*deltax)+f8*exp(9*1i*8*deltax)+f9*exp(9*1i*9*deltax)+f10*exp(9*1i*10*deltax)==cos(9*deltax)
U9 = 
U10=f1*exp(10*1i*deltax)+f2*exp(10*1i*2*deltax)+f3*exp(10*1i*3*deltax)+f4*exp(10*1i*4*deltax)+f5*exp(10*1i*5*deltax)+f6*exp(10*1i*6*deltax)+f7*exp(10*1i*7*deltax)+f8*exp(10*1i*8*deltax)+f9*exp(10*1i*9*deltax)+f10*exp(10*1i*10*deltax)==cos(10*deltax)
U10 = 
linsolve(U1,U2,U3,U4,U5,U6,U7,U8,U9,U10,f1,f2,f3,f4,f5,f6,f7,f8,f9,f10)
Error using sym/linsolve
Too many input arguments.
im trying to solve this system of equations
any help would be great

回答 (1 件)

Walter Roberson
Walter Roberson 2023 年 2 月 6 日
M=10
M = 10
deltax=2*sym(pi)/M
deltax = 
syms f1 f2 f3 f4 f5 f6 f7 f8 f9 f10
U1=f1*exp(1i*deltax)+f2*exp(1i*2*deltax)+f3*exp(1i*3*deltax)+f4*exp(1i*4*deltax)+f5*exp(1i*5*deltax)+f6*exp(1i*6*deltax)+f7*exp(1i*7*deltax)+f8*exp(1i*8*deltax)+f9*exp(1i*9*deltax)+f10*exp(1i*10*deltax)==cos(deltax);
U2=f1*exp(2*1i*deltax)+f2*exp(2*1i*2*deltax)+f3*exp(2*1i*3*deltax)+f4*exp(2*1i*4*deltax)+f5*exp(2*1i*5*deltax)+f6*exp(2*1i*6*deltax)+f7*exp(2*1i*7*deltax)+f8*exp(2*1i*8*deltax)+f9*exp(2*1i*9*deltax)+f10*exp(2*1i*10*deltax)==cos(2*deltax);
U3=f1*exp(3*1i*deltax)+f2*exp(3*1i*2*deltax)+f3*exp(3*1i*3*deltax)+f4*exp(3*1i*4*deltax)+f5*exp(3*1i*5*deltax)+f6*exp(3*1i*6*deltax)+f7*exp(3*1i*7*deltax)+f8*exp(3*1i*8*deltax)+f9*exp(3*1i*9*deltax)+f10*exp(3*1i*10*deltax)==cos(3*deltax);
U4=f1*exp(4*1i*deltax)+f2*exp(4*1i*2*deltax)+f3*exp(4*1i*3*deltax)+f4*exp(4*1i*4*deltax)+f5*exp(4*1i*5*deltax)+f6*exp(4*1i*6*deltax)+f7*exp(4*1i*7*deltax)+f8*exp(4*1i*8*deltax)+f9*exp(4*1i*9*deltax)+f10*exp(4*1i*10*deltax)==cos(4*deltax);
U5=f1*exp(5*1i*deltax)+f2*exp(5*1i*2*deltax)+f3*exp(5*1i*3*deltax)+f4*exp(5*1i*4*deltax)+f5*exp(5*1i*5*deltax)+f6*exp(5*1i*6*deltax)+f7*exp(5*1i*7*deltax)+f8*exp(5*1i*8*deltax)+f9*exp(5*1i*9*deltax)+f10*exp(5*1i*10*deltax)==cos(5*deltax);
U6=f1*exp(6*1i*deltax)+f2*exp(6*1i*2*deltax)+f3*exp(6*1i*3*deltax)+f4*exp(6*1i*4*deltax)+f5*exp(6*1i*5*deltax)+f6*exp(6*1i*6*deltax)+f7*exp(6*1i*7*deltax)+f8*exp(6*1i*8*deltax)+f9*exp(6*1i*9*deltax)+f10*exp(6*1i*10*deltax)==cos(6*deltax);
U7=f1*exp(7*1i*deltax)+f2*exp(7*1i*2*deltax)+f3*exp(7*1i*3*deltax)+f4*exp(7*1i*4*deltax)+f5*exp(7*1i*5*deltax)+f6*exp(7*1i*6*deltax)+f7*exp(7*1i*7*deltax)+f8*exp(7*1i*8*deltax)+f9*exp(7*1i*9*deltax)+f10*exp(7*1i*10*deltax)==cos(7*deltax);
U8=f1*exp(8*1i*deltax)+f2*exp(8*1i*2*deltax)+f3*exp(8*1i*3*deltax)+f4*exp(8*1i*4*deltax)+f5*exp(8*1i*5*deltax)+f6*exp(8*1i*6*deltax)+f7*exp(8*1i*7*deltax)+f8*exp(8*1i*8*deltax)+f9*exp(8*1i*9*deltax)+f10*exp(8*1i*10*deltax)==cos(8*deltax);
U9=f1*exp(9*1i*deltax)+f2*exp(9*1i*2*deltax)+f3*exp(9*1i*3*deltax)+f4*exp(9*1i*4*deltax)+f5*exp(9*1i*5*deltax)+f6*exp(9*1i*6*deltax)+f7*exp(9*1i*7*deltax)+f8*exp(9*1i*8*deltax)+f9*exp(9*1i*9*deltax)+f10*exp(9*1i*10*deltax)==cos(9*deltax);
U10=f1*exp(10*1i*deltax)+f2*exp(10*1i*2*deltax)+f3*exp(10*1i*3*deltax)+f4*exp(10*1i*4*deltax)+f5*exp(10*1i*5*deltax)+f6*exp(10*1i*6*deltax)+f7*exp(10*1i*7*deltax)+f8*exp(10*1i*8*deltax)+f9*exp(10*1i*9*deltax)+f10*exp(10*1i*10*deltax)==cos(10*deltax);
eqns = [U1;U2;U3;U4;U5;U6;U7;U8;U9;U10];
vars = [f1;f2;f3;f4;f5;f6;f7;f8;f9;f10];
[A,b] = equationsToMatrix(eqns, vars);
sol = A\b;
vars(:) == sol
ans = 

カテゴリ

Help Center および File ExchangeSymbolic Math Toolbox についてさらに検索

タグ

製品


リリース

R2022b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by