Ошибка: Subscripted assignment dimension mismatch. Почему ?

@khiggs2013 · Регистрация: 28.09.2016

Author24 — интернет-сервис помощи студентам

m-file1

Matlab M

function z=Bz(x0,y0,a,h,N)
psi=0:pi/50:2*pi*N;
F=f(psi,x0,y0,a,h);
z=trapz(psi,F);

m-file2

Matlab M

function F=f(psi,x0,y0,a,h)
Nx=length(x0);
Ny=length(y0);
for i=1:Nx
    s=0;
    for j=1:Ny      
        s=s+((y0(j)-a*sin(psi))*a.*sin(psi)+a*cos(psi).*(x0(i)-a*cos(psi)))./...
(x0(i)^2+y0(j)^2+a^2+(h/(2*pi)*psi).^2-...
2*(x0(i)*a*cos(psi)+a*y0(j)*sin(psi))).^(3/2);
     end;
       M(i,j)=s;
end;
F=M;

код программы:

Matlab M

N1=21;
i=1:N1+1;
j=1:N1;
x0min=0; y0min=0;
x0max=5; y0max=5;
x0(i)=x0min+(x0max-x0min)/N1*(i+1);
y0(j)=y0min+(y0max-y0min)/N1*(j+1);
a=1;
h=0.1;
N=1;
[x01,y01]=meshgrid(x0,y0);
mp=Bz(x0,y0,a,h,N);
contourf(x0,y0,mp,50);

Ошибки:
>> Test2
Subscripted assignment dimension mismatch.

Error in f (line 11)
M(i,j)=s;

Error in Bz (line 4)
F=f(psi,x0,y0,a,h);

Error in Test2 (line 13)
mp=Bz(x0,y0,a,h,N);

@Krasme · 20.06.2018, 11:01

Сообщение от khiggs2013

M(i,j)=s;

это должно быть внутри цикла

Добавлено через 49 секунд
на форуме есть теги matlab [MATLAB ] [ /MATLAB] - четвертый справа в первом ряду

@khiggs2013 · 21.06.2018, 12:35 **[ТС]**

Ребята, кто может подсказать что сделать с этим?

@VSI · 21.06.2018, 13:26

khiggs2013, а вот что...

Сообщение от khiggs2013

M(i,j)=s;

это должно быть внутри цикла

@khiggs2013 · 21.06.2018, 13:46 **[ТС]**

Дык я её внутрь цикла и помещал, то есть выше end , но результата то нет

@VSI · 21.06.2018, 13:52

khiggs2013, а ничего, что там есть еще один, вложенный цикл?

@khiggs2013 · 21.06.2018, 13:55 **[ТС]**

Дак что сделать то технически надо?

@SSC · 21.06.2018, 15:22

Сообщение от khiggs2013

Дык я её внутрь цикла и помещал, то есть выше end , но результата то нет

Конечно присвоение
M(i,j)=s;
должно быть выше еще одного end судя по индексу j, однако проблема в другом.
У Вас s получается массивом 1х101, из-за того что psi массив. И тут надо разбираться что Вы там и как интегрируете. Мне физика Вашей задачи неизвестна.
Из синтаксиса MATLABa возможно поставить присвоение

Matlab M

M{i,j}=s;

Однако перед интегрированием надо разбираться что с полученными данными делать.

@khiggs2013 · 22.06.2018, 00:19 **[ТС]**

Сообщение от SSC

У Вас s получается массивом 1х101,

А не могли бы вы пояснить как вы подсчитали массив что размерность его 1 на 101?

@SSC · 22.06.2018, 08:18

Сообщение от khiggs2013

как вы подсчитали массив что размерность его 1 на 101?

Вот эта строка

Сообщение от khiggs2013

Matlab M

1	psi=0:pi/50:2piN;

дает 101 элемент.
Но можно и не считать, а просто в отладчике поставить точку останова, и посмотреть размер массива

@khiggs2013 · 22.06.2018, 19:20 **[ТС]**

Сообщение от SSC

дает 101 элемент

Ага. понял. Ну скорее наверное 101*N так как N в моей задаче это число витков соленоида.
Дело в том, что все должно там вроде бы работать(где происходит интегрирование). Выражение

Сообщение от SSC

1
psi=0

i/50:2*pi*N;

определяет в функции

Сообщение от khiggs2013

z=trapz(psi,F)

пределы интегрирования, как я понял из: http://old.exponenta.ru/soft/m... /trapz.asp

Вообщем в моей задаче нужно взять например M(5, 6) (ну вот это:

Сообщение от khiggs2013

M(i,j)=s;

) и проинтегрировать его от 0 до 2*pi*N. Далее вынести на сетку по X=5, а по Y=6. И так перебрать каждый i и j.

@SSC · 22.06.2018, 21:38

Сущность проблемы мне не понятна, но возможно надо проводить интегрирование сразу перед присвоением
M(i,j)=s;
и в массиве хранить уже результаты интегрирования

@khiggs2013 · 23.06.2018, 17:02 **[ТС]**

Спасибо за помощь SSC!
Все таки надо было сначала интегрировать а потом уж перебирать по i , j.
Ответ получил в виде вот такого распределения. Это компонента по z магнитного поля однослойного соленоида рассчитанная по закону Био Савара в плоскости z-y.

@Krasme 6890 / 4940 / 2086 Регистрация: 02.02.2014 Сообщений: 13,126
	20.06.2018, 11:01	2
	Сообщение от khiggs2013 M(i,j)=s; это должно быть внутри цикла Добавлено через 49 секунд на форуме есть теги matlab [MATLAB ] [ /MATLAB] - четвертый справа в первом ряду 0

@khiggs2013 0 / 0 / 0 Регистрация: 28.09.2016 Сообщений: 10
	21.06.2018, 12:35 [ТС]	3
	Ребята, кто может подсказать что сделать с этим? 0

@VSI Модератор $Эксперт по математике/физике$ 5272 / 4054 / 1389 Регистрация: 30.07.2012 Сообщений: 12,429
	21.06.2018, 13:26	4
	khiggs2013, а вот что... Сообщение от khiggs2013 M(i,j)=s; это должно быть внутри цикла 0

@khiggs2013 0 / 0 / 0 Регистрация: 28.09.2016 Сообщений: 10
	21.06.2018, 13:46 [ТС]	5
	Дык я её внутрь цикла и помещал, то есть выше end , но результата то нет 0

@VSI Модератор $Эксперт по математике/физике$ 5272 / 4054 / 1389 Регистрация: 30.07.2012 Сообщений: 12,429
	21.06.2018, 13:52	6
	khiggs2013, а ничего, что там есть еще один, вложенный цикл? 0

@khiggs2013 0 / 0 / 0 Регистрация: 28.09.2016 Сообщений: 10
	21.06.2018, 13:55 [ТС]	7
	Дак что сделать то технически надо? 0

@khiggs2013 0 / 0 / 0 Регистрация: 28.09.2016 Сообщений: 10
	22.06.2018, 00:19 [ТС]	9
	Сообщение от SSC У Вас s получается массивом 1х101, А не могли бы вы пояснить как вы подсчитали массив что размерность его 1 на 101? 0

@khiggs2013 0 / 0 / 0 Регистрация: 28.09.2016 Сообщений: 10
	22.06.2018, 19:20 [ТС]	11
	Сообщение от SSC дает 101 элемент Ага. понял. Ну скорее наверное 101N так как N в моей задаче это число витков соленоида. Дело в том, что все должно там вроде бы работать(где происходит интегрирование). Выражение Сообщение от SSC* 1 psi=0i/50:2piN; определяет в функции Сообщение от khiggs2013 z=trapz(psi,F) пределы интегрирования, как я понял из: http://old.exponenta.ru/soft/m... /trapz.asp Вообщем в моей задаче нужно взять например M(5, 6) (ну вот это: Сообщение от khiggs2013 M(i,j)=s; ) и проинтегрировать его от 0 до 2piN. Далее вынести на сетку по X=5, а по Y=6. И так перебрать каждый i и j. 0

@SSC $Эксперт по математике/физике$ 3390 / 1913 / 571 Регистрация: 09.04.2015 Сообщений: 5,365
	22.06.2018, 21:38	12
	Сообщение было отмечено khiggs2013 как решение Решение Сущность проблемы мне не понятна, но возможно надо проводить интегрирование сразу перед присвоением M(i,j)=s; и в массиве хранить уже результаты интегрирования 1

@khiggs2013 0 / 0 / 0 Регистрация: 28.09.2016 Сообщений: 10
	23.06.2018, 17:02 [ТС]	13
	Спасибо за помощь SSC! Все таки надо было сначала интегрировать а потом уж перебирать по i , j. Ответ получил в виде вот такого распределения. Это компонента по z магнитного поля однослойного соленоида рассчитанная по закону Био Савара в плоскости z-y. Миниатюры 0

Ошибка: Subscripted assignment dimension mismatch. Почему ?

Решение