sergsh_sergsh
Newbie | Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Цитата: Ну раз уж вы вроде как новичок то вот....
Код:
size=3;
A=magic(size)+0.5;
display(A);
fid=fopen('mytest.mat','wb');
fwrite(fid,size,'int32');
fwrite(fid,A,'float32');
fclose(fid);
Код:
program readnum
implicit none
integer*4 :: size
real*4,dimension(3,3) :: A
open(unit=100,file='mytest.mat',form='binary');
read(100) size
read(100) A
close(unit=100)
print *,'Size=',size
print *,'A=',A
end program readnum
|
Спасибо !!!!!!!!!!!!!
Жалко тут нельзя плюсы ставить, я бы вам поставил 100 штук !
в Гугле я искал долго, нашел несколько примеров и ни один у меня не заработал
Ваш фортран код заработал на Intel
а на gfortran написал ошибку
ошибка: FORM specifier in OPEN statement at (1) has invalid value 'binary'
но если сделать так :
open(unit=100,file='mytest.mat',form='unformatted', access='stream')
то и в GFORTRAN все читается
Intel Fortran у меня стоит, но Matlab его не видит.
Его переменные я запускаю "локально" : >source /opt/intel/bin/compilervars.sh intel64
и их Matlab видимо не находит
можно ли их как то указать явно (глобально) ? |
Всего записей: 31 | Зарегистр. 09-01-2012 | Отправлено: 22:35 22-06-2012 | Исправлено: sergsh_sergsh, 23:38 22-06-2012 |
|
Для практики в этом случае нужно как можно точнее выполнить численное дифференцирование. Классическая схема - по три точки на каждую переменную для первой-второй производных, плюс по точке на смешанную производную второго порядка, она при таком количестве параметров точно даст ужас, летящий на крыльях ночи, чисто из-за погрешностей такого примитивного дифференцирования. Если построить разумную схему с увеличенным количеством точек, возможно, и прокатит. 
